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PREFACE 


Engineering Aerodynamics has been written to meet a definite 
need for more practical information on aerodynamics, presented 
in forms suitable for direct application by aircraft designers 
and advanced students in aeronautical engineering. An endeavor 
has been made to include as much new material as possible 
while avoiding a needless repetition of well-known or easily 
obtained data. For this reason the conventional treatment of 
elementary aerodynamic considerations, such as airflow patterns 
and pressure distributions, is omitted, thus allowing more atten- 
tion to be concentrated on the problems encountered in design. 

A large amount of data has been analyzed and presented 
in the form of working diagrams and equations, while very little 
space has been devoted to the description of tests and test 
methods, or to the reproduction of undigested test data. Where 
given, derivations have been made as brief as practicable, but 
numerous references have been included for the benefit of those 
desiring more complete information. Many important calcula- 
tions are illustrated by examples. 

Most of the subject matter is original and much of this is 
now made available for the first time. A special effort has 
been made to present all data in the form best adapted to con- 
venient application in design problems. This has led to the 
use of numerous working diagrams, or graphs, from which 
required factors, constants, or ratios, may be readily obtained. 

The author deeply appreciates the encouragement received 
from many sources during the preparation of this volume. 

WALTER S. Drex#L. 


Washington, D. C., 
May 25, 1928. 
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ENGINEERING 
AERODYNAMICS 


CHAPTER I 
WING SECTION DATA 


The National Advisory Committee for Aeronautics has 
compiled a series of technical reports giving wind tunnel test 
data on more than 600 wing sections. These reports are issued 
under the title “Aerodynamic Characteristics of Airfoils,’’ and 
are numbered I to IV, as Report Nos. 93, 124, 182, and 244. 
Copies may be obtained directly from the Committee or from 
the Superintendent of Documents, Government Printing Office, 
Washington, D. C. 

The Committee also has published technical reports and 
technical notes on full scale airfoil characteristics, as measured 
in the variable density wind tunnel. These are of great value 
to the engineer: Technical Report No. 221, “Model Tests with 
a Systematic Series of 27 Wing Sections at Full Reynolds Num- 
ber,” by Max M. Munk and Elton W. Miller (1925); Technical 
Report No. 233, ‘‘Aerodynamic Characteristics of Seven Fre- 
quently Used Wing Sections at Full Reynolds Number,” by 
Max M. Munk and Elton W. Miller (1926). 

In view of the availability of wing section data, this chapter 
will give only such data as are necessary to illustrate general 
wing characteristics and the selection of a section best suited 
for a general design. 


Coefficients. Wing section characteristics are most conven- 
iently expressed in the form of “‘coefficients.” Practically all of 
the early work in this country was given in the form of coefficients 
having the dimensions of lbs./sq. ft. / (mi./hr.)?, using the 
symbols Ky for lift and K, for drag. These coefficients were 
based on air of standard density and defined by the relations 


I This expression frequently appears abbreviated as N.A.C.A. 
3 


- 
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Lien K,() SV? 


° 


Drag = D = K, ) mo 
where SS is the area in square feet and V the relative speed in 
miles per hour. 
About 1919 the National Advisory Committee for Aeronautics 
recommended the adoption of ‘‘absolute coefficients” having the 
same value in any consistent system of units. The coefficients 


tentatively adopted were defined by the relations 


ti Leo Ss 

D = Dy Q S We 
The coefficients L, and D, are the same as K; and Kp used hy 
the British. ‘, 


With the improvement in theoretical aerodynamics it became 
evident that there were many advantages in the use of the 
absolute coefficients employed by Prandtl. Therefore, in 1921, 
the National Advisory Committee for Aeronautics decided to 
recommend the general use of the coefficients Cz and Cp defined 
by the relations 


L= Cz, £SV?=CrqS andD=Cy)£ SV? =CpqS 


Where g is the dynamic pressure $e V*. These absolute coeffi- 
cients are now in extensive use. 
The relations between the three types of coefficients are: 
Gr = 2Le:= "201 1 Ky = .002558 Cr = .005116 Lc 
Cp =2Dc = 391K, Kz = .002558Cp = .005116 Dc 
Dynamic pressure. The dynamic or impact pressure 
qg=%0V” is the pressure developed in bringing a moving 
perfect fluid to rest. Since the standard density is eo = .002378, 


U 


Covt! Jo = -001189 V? 


where q» is in lbs./sq. ft., and V in ft./sec., or 
| Jo = -002558 V? 
where qo is in lbs./sq. ft. and V in mi./hr. 
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The value of g at any density other than the standard is 
obviously 
oe 2 
ane (=) 


In the metric system, V being in m./sec., and g in kg./sq. m. 
do = .0625 V? 


The reJation between velocity in standard air and a given 
dynamic pressure in lbs. /sq. ft., is 


V= 29.00 Vq_ ft./sec. 
V = 19.77 Vq_ mi./hr. 


In the metric system, g in kg./sq. m. 
V=4Vq m./sec. 


Wing section characteristics. Wing section characteristics 
are lift, drag, and moment coefficients as functions of angle of 
attack, Reynolds Number, absolute velocity ratio or any other 
non-dimensional ratio affecting these characteristics. Wing 
section characteristics as a function of angle of attack are usually 
determined by wind tunnel tests on an ‘‘airfoil”? model at a 
constant air speed. This model is normally of aspect ratio 6, 
and in this country, of the dimensions 3” x 18”, 5” x 30” or 6” x 36”, 
the smaller dimensions being the chord, or length in the direction 
of the relative wind. The variation of wing section character- 
istics with Reynolds Number cannot be fully determined in the 
average atmospheric wind tunnel because the combination of 
the largest model which can be tested with the highest test speed 
available or attainable gives a Reynolds Number which is of the 
order of 5% to 10% of the full scale value. The variable density 
wind tunnel is the most satisfactory method yet found for 
obtaining full-scale data. 

In comparing data from various wind tunnels, exact agree- 
“ment must not be expected. Good agreement should always be 
obtained however, if the tests are carefully made, particularly 
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with regard to corrections for the model attachment, and if the 
data are corrected for wall interference. (See Chapter II.) 
It has been found by experience that the test speed should be not 
less than 30 nor more than 120 miles per hour in order to avoid 
critical values of Reynolds Number and the effects of compress- 
ibility. Model chords of less than 3” are also likely to show the 
effects of critical flow conditions, due perhaps to unavoidable 
spindle or clip interference. 


Methods of plotting wing section data. Wing section char- 
acteristics in general may be plotted against angle of attack or 
lift coefficient. The former is convenient in presenting wind 
tunnel data, but the latter is more suitable for comparing various 
sections. 

The conventional method of plotting characteristics against 
angle of attack makes use of the “geometrical” angle of attack 
measured from an arbitrary chord line, which is usually the line 
of tangency to the lower camber of the section. The engineer 
should understand that this arbitrary reference line is adopted 
for convenience of definition and not because it has any con- 
nection with lift and drag characteristics. If the characteristics 
are plotted against the “absolute” angle of attack measured from 
the angle of zero lift, the difference between various sections is 
found to be much reduced. 

Moment coefficients should be plotted against lift coefficients, 


: Vole ee 
and center or pressure against speed range 75” Since in both cases 


the resulting curves will be approximately straight lines. The 
common plot of center of pressure against angle of attack is very 
difficult to read. 

Profile drags should be plotted against either lift coefficient 
or speed range. 


Selection of a wing section. The selection of the best wing 
section for any particular design depends upon a number of 
factors, the most important of which are: 
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Structural adaptability 

Stable airflow characteristics 

Low profile drag at all lift coefficients 

Low minimum drag 

High minimum lift 

Moderate moments, or a reasonable center of pressure 
travel 


By structural adaptability is meant the suitability of the 
section for the type construction under consideration. This 
usually sets a rather definite limit to the spar depths for an 
efficient structure. For example, in a biplane a much thinner 
section can be used than in a monoplane. The designer knows 
in advance the approximate spar depths desired and is able to 
set aside a group of wing sections for further investigation. 

Certain wing sections are very sensitive to slight mutilations 
or variations from the specified curvatures. Other sections can 
be modified considerably without appreciable change in char- 
acteristics. The difference is probably due primarily to some 
instability in the air flow around sections of the first type, and 
may be caused by irregularities in the variation of radius of 
curvature along the upper camber. In some cases the unstable 
flow is apparent from the conventional wind tunnel test and it is 
usually indicated by irregularities in pressure distribution curves. 

The requirement of either a low minimum drag or a high 
maximum lift serves to restrict the number of sections available 
for a given design. The final selection is then made on the basis 
of a low profile drag. A very convenient method of analyzing 


‘ ° Cro ) O ( Cr i 
section data is to plot (= Bs eeeinse Cr max! * 


(=) = 

CL max - W 

and hence is proportional to the total profile drag. This meth- 
od, therefore, compares the sections on the basis of total profile 


drag for constant gross load and stalling speed. The section 
selected will vary with the major requirement as follows: 
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Maximum speed. Section having least value of 


( =) at high speed ratios, is 25: 
Cz max S 


Maximum climb and ceiling. Section having least value 


Coo Via Vas 
of (es aa between ia = 7,10 and fee 1.50. 


Maximum endurance. Section having least value of 


Cos Vive 
Ges) at () = 1.10. 


General performance. Section having least average value 


of ( ae ) at all values of a. 


Lmax Vs 


A preliminary study of polar diagrams, moment curves and 
other data will ordinarily reduce to three or four of the number 
of sections to be analyzed for final selection. There is little to 
be gained by carrying out complete performance calculations in 
order to make a decision, but wind tunnel tests sometimes show 
a difference between two models which cannot be predicted from 
the wing section data. This difference is probably an inter- 
ference effect which properly belongs under the heading of stable 
airflow characteristics. For example, certain sections such as 
the U.S.A.-27 and G-387 do not show up as well in an airplane 
model test as do the G—398 and its derivatives N—1o and N-22, 
or the R.A.F.-15 and Clark Y. 


Example of analysis of section data. Values of 2 °_ for the 
Lmax 


U.S.A.-27, Clark Y, G-398, R.A.F.-15, and N-22 sections are 
plotted against Fr. on Figure 1. These data were obtained on 


5” x 30” models at 40 miles per hour in the 8’ x 8’ Washington 
Navy Yard Wind Tunnel and have been corrected for wall 
interference. This figure brings out the superiority of the G—308 
and N-z2 sections (the latter is 90% of the G-398 reduced on 
the chord line). For climb and ceiling the N-22 is exceptionally 
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Cp, , Ve 


from Washington Navy Yard Tests 
Cr max Vs e 1 


004 


CD. with 7 from N. A. C. A. Tests 


Figure 2. Variation of Full Scale 
Cr max 


good, while for all-around performance the G—398 shows up well. 
The poor showing of the R.A.F.—15 is interesting in view of its 
general excellence on the basis drag per unit area. 

This comparison must be made with full-scale data if the 
Cr, 
Cr max 


conclusions are to be reliable. Figure 2 is a plot of against 
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Io 


Figure 3. Effect of Reynolds Number on C Tear. N. Ay Cy Aeateste 


A. C. A. Tests 


Figure 4. Effect of Reynolds Number on Chang nk 
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i: for the U.S. A.-27, Clark Y, R.A. F.-15, and M-6 sections, 
Ss 


using data obtained in the N.A.C.A. variable density wind 
tunnel at approximately full-scale Reynolds Number. It appears 
that the first three sections are almost identical, but the M-6 is 
poor for climb and good for high speed. Unfortunately, similar 
data on the G-398 and N-22 are not now available so that the 
comparison can not be extended to these sections. However, 
preliminary tests indicate that the scale effect on the G—398 is 
similar to that on the Clark Y and flight tests indicate excep- 
tionally good full-scale characteristics. 
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ure 5. Full-Scale Profile Drag Coefficient. N. A.C. A. Tests 
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Variable density wind tunnel data. The principle of Munk’s 
variable density method for securing data at high Reynolds 


Oo 4 #48 I2° 16° 20° 24° 


Figure 6. Lift Coefficient at a= Ooas a Function of Trail Angle 8 


Number is well known. Comparison of flight test data with the 
results obtained in the N. A.C. A. variable density wind tunnel 
have shown such excellent agreement that there can be no doubt 
of the validity of this method. pis 
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Figure 3 gives the variation of Cz max and Figure 4 gives the 
variation of Cp min with Reynolds Number as obtained in the 
variable density tests on a series of popular wing sections. 
Figure 5 is a plot of profile drag Cp, against lift coefficient for 
various sections, using test data at Reynolds Numbers of 
approximately 3.6 X 10°, 


2.00 


8° 2°, PuGz aN 20”) 24° 


Figure 7. Maximum Lift Coefficient as a Function of Trail Angle B 


Trail angle. The acute angle which the tangent to the mean 
camber at the trailing edge of an airfoil section makes with the 
chord line of an airfoil is known as the trail angle. (See Figure 9.) 
This angle is a convenient substitute for the mean camber and 
certain characteristics may be referred to it for approximate 
estimates, etc. The effect of trail angle on lift at zero angle 
of attack (of the chord line) is given in Figure 6. The effect of 
trail on maximum lift coefficient is given in Figure 7, and the 
effect of trail angle on center of pressure travel is given in Figure 
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8. The center of pressure is given to a fair degree of approxima- 
tion by the empirical equation: 


Cp = .17 + .075 (~) + .00645 ()"o (1) 


These data apply with a fair degree of accuracy to all normal 
wing sections, including those with flaps on the trailing edge. 


90 
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Figure 8. Variation of Cp with Trail Angle and Speed Range 


Mean camber. The mean camber of an airfoil section lies 
midway between the upper and lower cambers. (See Figure 9.) 
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The mean camber has a greater effect on the aerodynamic 
characteristics of an airfoil than any other physical dimension. 
The general airfoil sections of the same mean camber have the 
same or very similar aerodynamic characteristics provided that 
the thickness of the section is not greater than 20% of the chord. 


UPPER CAMBER 


—_—_——nks oe 


LEADING \LoWwER CAMBER wees LINE TRAILING 


EDGE EDGE 
SHOWING “TRAIL ANGLE’ 6 


Figure 9. Illustrating Mean Camber and Trail Angle 


Wind tunnel tests at the Washington Navy Yard have shown that 
a thin section, like the well-known ‘‘Sloane,’’ may be increased 
60% in thickness in this manner without appreciably affecting 
the general aerodynamic characteristics. 

Wing flaps. Munk has shown’ that the effect of turning up 
or down a flap (or elevator) is 


Aa = K 28 Be 


where Aa, is the change in angle of attack of the entire surface 
which has the same effect as turning the flap (or elevator) to the 


angle 8 measured positive downward, “the ratio of the flap chord 
to the total chord, and K a constant which varies with the ratio 


! In the original paper, Munk gives the derivation of K and a 
C 


2 Max M. Munk “General Theory of Thin Wing Sections,” N. A. C. A. Technical 
Report No. 142 (1922).. 
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Effect of Flaps on Lift at C, = 1.13 


Figure 11. 
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plot of K against ~ Between -= to and 4= .60, which covers 
(s 


the whole practical range, K is given closely by 


s t 
K = 1.90 — 0.57 ~ (3) 
Equation (2) may be written in the more convenient form: 
eae 
AC, = ay K ; (2a) 


Analysis of a number of wind tunnel tests verifies Munk’s 
equation for the small values of @ on which it is based, but for 
large values of 8, K is considerably reduced. A plot of K 
against 8, as determined from the experimental data in R. & M.3 
Nos. 110, 152, and 3109, indicate that the variation of K with 6 
is approximately linear and according to 


° i ° 
Ke= K;— .0148° = 1.90 — hey (a .O14 8 (3a) 


where K, is the actual value and K, the theoretical value. 

More recently Glauert* has derived some theoretical relations 
for lift and moment of wings with flaps, which are somewhat 
complicated but appear to give excellent results. 


The engineer is interested primarily in flap hinge moments | 


and the changes in lift and drag due to flaps. For this purpose | 


a series of figures has been prepared from test data. 

The change in lift due to flaps of various chords set at various 
angles is shown on Figures 10 and 11. The former is fora lift 
coefficient Cz = 0.23 and the latter for a lift coefficient of 
Cz = 1.13, corresponding to angles of attack of o° and 12° on 
the R.A.F.-9 section. Figures 12 and 13 are companion curves 
showing the change in drag. These curves may be used to 
calculate rolling and yawing moments of ailerons, by assuming 
that only the wing area immediately forward of the aileron is 
affected together with the aileron. 

3 Reports and Memoranda, British Aeronautical Research Committee. 

4H. Glauert, “ Theoretical Relationships for an Aerofoil with Hinged Flap,” Br. A.R. C. 


R. & M. No. 1095 (1927). 
5 Taken from R. & M. Nos. 110 and 319. 
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Figure 13. Effect of Flaps on Drag at C TL = 1.13 
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Hinge moments are shown on Figures 14 and 15 for the same 
conditions applying to lift and drag. The hinge moment in 
Ibs.-ft. is given by 


H = Ku *acv (4) 


where a is the flap area (sq. ft.), c the flap chord (ft.), and V the 
air speed in miles per hour. 

These plots are of considerable interest in that they bring 
out the decrease in effectiveness of ailerons at large lift coeffi- 
cients and show the origin of the yaw due to ailerons. They 
also show how rapidly the hinge moment increases with increase 
of flap chord. 

' The increase in maximum lift which can be secured by the 
use of flaps is practically independent of the original wing 
section, although it varies slightly with the camber ratio and 
Ci mox- -It is approximately ACz me, = 0.70 ot AKy mox = 
.0o18 for a thin low-lift section, and ACz max = 0.60 or LI ae 
= .oo15 for a thick high-lift section. Flaps therefore show up 
best on thin low-lift sections. The only practical use of flap 
is to increase the maximum lift and decrease the stalling speed. 
It is very unlikely that climb and ceiling can ever be ap- 
preciably improved by the use of flaps of the conventional form 
on a good wing section. However, there appears to be some 
field for the future use of very narrow balanced flaps extending 
over the entire span and fitted with a gear which permits differ- 
ential operation as ailerons from any flap setting. ‘The purpose 
of such an arrangement is to reduce the flap forces and weight, 
and the only gain to be anticipated is a decrease in stalling speed. 

Slots. The Lachman-Handley Page slot is well known from 
the great amount of publicity that it has received as a means of 
increasing lift and preventing loss of control in a stall. The 
latest development is the automatic slot which opens without 
external control, by the reduction of pressure on the leading edge 
airfoil. This slot can be designed to open at an angle of attack 
just below the stall, and when fitted to the wing forward of the 
ailerons it gives remarkable lateral control at very high angles. 
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Figure 15. Flap Hinge Moment Coefficients Ky X 105 at Cy = 1.13 
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As a means of increasing lift and reducing stalling speed, the 
slot is very effective, although somewhat disappointing in that 
the full-scale net gain is much less than might be expected from 
wind tunnel tests on airfoil models with a slot across the entire 


200 
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1.20) 
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16° 24 32° 40° 
Figure 16. Effect of Slots and Flap on Lift Coefficient 


span. The loss may be ascribed to two causes. In most 
practical applications it is not feasible to extend the slot across 
the entire span or to keep the slot clear of obstructions, and any 
interruption to the slot or to the flow through it causes a very 
marked reduction in lift. The second cause for failure to obtain 
the full benefit of the slot is found in the very high angles of 
attack required for maximum lift. The slot merely extends the 
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Figure 17. Effect of Cutting off Trailing Edge 
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lift curve against angle of attack without appreciably changing 
the slope, and maximum lift may occur between 23° and 20°. 

Slots are most effective when used in conjunction with flaps 
as shown by Figure 16, which is based on wind tunnel tests on a 
model wing fitted with two slots and a flap. If the flap alone 
were used at an angle of attack of 16° there would actually be a 
decrease in lift coefficient from Cz =1.50to Cz =1.30 approx- 
imately, but if slots are used with the flaps, the lift coefficient 
increases from Cp =1.50toC =2.18. 


Wing tips; cut-outs. When the chord of a wing varies, the 
thickness of the section should vary in proportion, particularly 
at the tips. If the wing section is originally thick, an improve- 
ment is obtained by taper in camber as well as in chord, at the 
tips. 

For cut-outs at the center section, very little is to be gained 
by making the sections geometrically similar, and a fairing to 
the trailing edge is sufficient. Figure 17 shows how little the 
lift is affected by cutting off the trailing edge of a wing. The 
section drag at all lift coefficients for any particular amount 
removed is increased by approximately 1/10 of the reduction in 
C, based on the original area. 

The effect of center section cut-outs on wing characteristics 
has been investigated at the Washington Navy Yard°® and the 
conventional approximately semicircular form with rounded fair- 
ing into the normal trailing edge, was found quite satisfactory. 


6R. H. Smith, ‘‘Air Force and Moment for N-20 Wing with Certain Cut-Outs,” 
N. A.C. A. Technical Report No. 266 (1927). 


CHAPTER 
WING THEORY 


Hydrodynamical definitions. Before giving an outline of 
some of the outstanding developments in theoretical aerody- 
namics, it is desirable to give a brief review of the hydrodynamical 
conceptions on which the subject of modern aerodynamics is 
based. It is not possible in a work of this type to give more than 
a definition of the more common terms, but every aeronautical 
engineer should have a clear conception of what is meant by 
potential flow, circulation, vorticity, etc. 

For further information on theoretical hydrodynamics, 
reference should be made to Prandtl’s “Application of Modern 
Hydrodynamics to Aeronautics” (N. A. C. A. Technical Report 
No. 116), or to the following books: 


E. B. Wilson, “Aeronautics,”” John Wiley & Son (1920). 

.“The Mechanical Properties of Fluids.”” <A collective 
work. D. Van Nostrand (1924). 

H. Glauert, “The Elements of Aerofoil and Airscrew 
Theory.’’ Cambridge University Press. 

C. Eberhardt, ‘“‘Einfuhrung in die Theoretische Aero- 
dynamik.’’ R. Oldenbourg, Munchen (1927). 


H. Lamb, “Hydrodynamics.”” Cambridge University 
Press. 


Each of these books contains much of value to the student 
but the last three are particularly valuable. The beginner who 
can read German will find Eberhardt’s book very clear and 
instructive. Lamb’s Hydrodynamics is the classical general 
treatment of this subject. 


Bernoulli’s theorem. Bernoulli’s theorem states that the 
24 
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total energy of a fluid particle is constant at all points on its 
path in a steady flow. Bernoulli’s equation is 


ee 
ge tieVth=H (5) 


r: 
where a the “‘pressure head,” $e V? the “velocity head,” 


h the potential head, and H a constant. 
This equation is due to Daniel Bernoulli and dates from 1737. 
In the original form and as given it applies only to steady flow 
of an incompressible fluid but many of the more common problems 
of hydrodynamics are solved by its use, with the following 
restrictions: 
For a general fluid in general motion, H is never constant. 
For an incompressible fluid in general motion, H is con- 
stant for a given particle. 
For an incompressible fluid in steady motion, H is con- 
stant for all particles along a streamline. 
For an incompressible fluid in steady irrotational motion 
H is constant for all particles throughout the fluid. 


Stream function. Consider the flow of an incompressible 
fluid across an element dS of any curve in the plane of a two- 
dimensional flow. The flow across dS must be equal to the sum of 
the flows across dx and dy. The flow across dx is v dx and that 
across dy is —u dy, hence the flow across the element dS is 

dy = vdx — udy 
and the difference between the values of » at any two points 
along the curve is given by the integral | 
| S dy = f (vdx — udy) , 
v is called the “stream function.”’ 

Streamline. Mathematically a streamline is defined by 
 =C, or dy =o, since for this condition there is no flow across 
the curve so defined. 

Two-dimensional flow. A flow which is two-dimensional in 
the plane of X-Y will be exactly similar in any parallel plane. 
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An example of such a flow is that around the center of a very 
long cylinder or strut. ‘Two-dimensional flows can be completely 
investigated by considering the flow in a single plane. 

The mathematical conditions for two-dimensional flow are: 
(x) no velocity component along the z axis, and (2) no xory 
velocity gradients along the z axis. That is 


Ve a 
BVe _ 3Vs 
og oz 


Three-dimensional flow. The general motion of a fluid is 
three-dimensional, with acceleration and velocity components 
along all three axes. 


Velocity potential. If a fluid flow is such that the circulation 
about every closed curve vanishes, then the flow may be repre- 
sented by 

udx +udy = —d® 
where @ is a function of x and y, known as the velocity potential. 
Any fluid motion for which there is a velocity potential is said 
to be irrotational, or a potential flow. 

The velocity potential @ and stream function » occur in 


three-dimensional flow, in exactly analogous senses as in two- 
dimensional flow. 


Circulation. Circulation in a fluid is determined by the line 


integral of the tangential velocity taken around any closed circuit. 
That is 


T= Sf V;cos6-dS 


where V, is the resultant velocity making an angle 6 with the 
elements dS. Resolving V, into its components gives 


yy ee. dy 
V,cos 6 = (uagte a) 


hence 


d d 
ELEY e («sat @aa) dS = f (udx +vdy) 
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Many engineers instinctively associate circulation with lift. 
It should be understood that circulation in general is a type of 
fluid motion and that in any particular case it has the value given 
by the integral in the preceding paragraph. 

The lift on an element of unit length of a wing of infinite span 
varies directly as the circulation around it. That is 


L=elVv 


where ¢ is the density and V the relative velocity measured at 
‘a great distance from the element. This is the well-known 
Kutta-Joukowski equation. 

For a wing of finite span, the lift and from it the induced 
drag can be calculated only when the distribution of circulation 
along the span is known. If the circulation is constant along the 
span, then $l is zero except at the tips where it equals I. This 
would correspond to a constant lift along the span, and there 
would be a vortex at each tip only. Actually there is a vortex 
at each tip and a gradient of I’ along the span from a maximum 
value at the center to zero at the tips. The vortices which peel 
Bt 
db 
variation in vorticity along the span is from a maximum positive 
value at one tip to an equal negative value at the other tip, 
passing through zero at the center. 


off of the trailing edge vary in strength with —. Hence the 


Vorticity-vortices. The circulation around an infinitesimal 
area is proportional to the area, and the ratio of circulation to 
the area is the measure of the intensity of rotation or vorticity. 

The so-called ‘‘ Vortex Theory” of Prandtl, by means of which 
solutions have been made for many problems involving wing lift 
and drag, is perhaps better considered as a group of theorems 
explaining the action of vortices and the method of determining 
the resultant velocity at any point from a given distribution of 
vortices. 

In this connection a vortex is a mathematical rather than a 
physical conception. The violent whirling motion in water, 
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which is called a vortex, certainly has some rotation although it 
may be irrotational. It would be irrotational if the velocity 
varied inversely as the radius, and it is rotational to the extent that 
the velocities deviate from the inverse law. The simplest case 
of rotation is a mass of fluid rotating as a rigid body, and for this 
case the rotation or vorticity is twice the angular velocity. 

The mathematical vortex is a convenient method and not a 
necessity as Munk? has shown. The same results are obtained 
by kinetic energy relations. 


Induced drag. In 1911, Dr. Prandtl and his assistants at 
Gottingen derived a relation between the circulation and the 
vertical or downwash velocity component w, due to a lifting 
wing. At the same time it was proved that half of the final 
downwash velocity was acquired forward of the center of pressure, 
or in other words, a downward acceleration was imparted to the 
air at some distance forward of a lifting wing which therefore 
operated in a downwardly inclined airstream. The average 
downwash velocity at the center of pressure was found to be 
greater, the greater the lift, and the nearer to the center of the 
wing the main production of vortices. 

As a result of the virtual inclination of the airstream through 
the angle 9=tan™ a , the wing ‘“‘goes up a hill” having the 
slope 9. Consequently the lift, which is vertical to the re- 
lative wind, now has a rearward or drag component. With- 
out going into the details of the derivation? it was shown 
that the inclination of the lift vector, and hence the drag, 
was a minimum when the downwash was constant along the 
span, and that this condition corresponded to a lift distri- 
bution proportional to the ordinates of an ellipse having the span 
as a diameter. The constant value of the downwash velocity 


* Max M. Munk, ‘‘Elements of the Wing Section Theory and of the Wing Theory,” 
N. A.C. A. Technical Report No. ror (1924). 


? See L. Prandtl, ‘‘Applications of Modern Hydrodynamics to Aeronautics,” N. A. C. A. 
Technical Report No. 116. 


3 The mathematical proof was first given by Munk in his Géttingen Dissertation which 
was subsequently translated and published as N.A.C.A. Technical Report No. 121 
“The Minimum Induced Drag of Aerofoils.”’ 
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resulting from the elliptical lift distribution was shown to be 
= Po 

2b 
where I’, is the circulation at the center of the span b. It may 
be shown that the value of T, is 
wna 

ToVb Take 

where L is the lift and pe the density. Hence the angle of 
downwash is 


DP, 


ae 2L 
~ to Vb? 
Since the downwash is constant along the span the drag D is 
w aie iis 
ng aie age (6) 


It was found in 1913 that the actual measured drags were 
greater than this theoretical minimum but the wing sections in- 
vestigated were very poor. Two years later, investigations on 
much better wing sections showed close agreement between the 
theoretical and the measured drag. ‘The investigation was then 
continued on wings of various aspect ratios and a very important 
discovery made. At the same lift coefficient with the same 
wing section the difference between the theoretical and measured 
drag coefficients was always the same, within the experimental 
error of course, for any aspect ratio. The part of the total drag 
which does not change with aspect ratio is due to the shape of 
the wing section and hence it was given the designation “Profile 
Drag”’ or “Section Drag.”’ The theoretical drag, varying only 
with lift and span, was given the designation “Induced Drag” 
from the analogy to electrical induction phenomena. This 
constituted one of the most important advances ever made in 
theoretical aerodynamics and opened an immense field of 
practical application. 


Substituting Cas SV? for L and Cp ESV? for D in equation 
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(6) gives the coefficient of induced drag 
Coi= cee (7) 
TD 


The section drag is Cp>=Cp—Cp; and this is constant for any 
given section and lift coefficient. Hence at a constant value of 
C_ the relation between the drags for two aspect ratios is 


Cr AJ acs Crs; 
Chas = mb,” Co, Tb”? 
or 
Coz — Cpr + cl Ss Ss =) (8) 


which enables the drag to be calculated for any aspect ratio 
when the drag is known for one aspect ratio. 


Induced angle of attack. Since the wing is operating in a 
vertical downwash velocity of 
_ 2h, 
mp Vb? 
and a horizontal velocity of V, the relative wind direction is 
-inclined downward to the rear of the wing by the angle 9 having 


the value 
=o (2) = tn (2) 
9 = tan e tan m0 V0? 
? is always small so that tan 9 = 9, hence 


oe 2L fey C DS 
fei moV?b? xb? (9) 


9 


g is the “induced angle of attack.’ It increases as the aspect 
ratio decreases. The physical significance is that as the aspect 
ratio is decreased the downwash increases and the wing must be 
turned up to a higher apparent or geometrical angle of attack in 
order to obtain a given lift coefficient. 


Munk’s span factor. Equation (7) was completed by Munk, 
who showed that in order to apply to biplanes (or multiplanes) 
the maximum span 6 must be replaced by kd, which is the span 


4 
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of the monoplane having the same area and induced drag as the 
biplane (or multiplane). For a monoplane k = 1.00, but for a 
biplane & varies with the ratio of gap to span, the ratio of the 
spans and the proportional area in the two wings as will be shown 
later. 

The introduction of the span factor was of great practical 
importance. With this factor equation (7) becomes 


go OS 
Coi = m (kb)? (7a) 
and equation (8) becomes 
= Ce AYA a Se 8 
Co; — Cr: = [ ahs at | ( a) 
In a similar manner, the induced angle of attack in radians is 
A 1Gps 
OL n(kb)? (9a) 


and in degrees 
o 57-3CL 1S, Ss: | 
G2 _— Gy + . la aye (k,b,)? (10) 


Equation (10) is only approximately correct. Munk com- 
pleted it later+ by dividing into three parts the angle of attack 
necessary to produce a given lift coefficient. These parts are as 
follows: (a) the intrinsic angle of attack for the given wing 
section and lift coefficient, (6) the additional induced angle of 
attack, and (c) the additional interference angle of attack. 
With this modification equation (10) becomes 


= A+ s7Ce| (5, + 1.) ~ (ais: + 1) | (10a) 


where J is the interference factor. J varies slightly with stagger 
and with wing section, and is less for a lift produced by curvature 
than for lift produced by angle of attack. However, I is approx- 
imately a function of gap chord ratio only, with the following 
average values according to Munk: 


4 Max M. Munk, “General Biplane Theory,’’ N. A.C. A. Technical Report No. 151 
(t922). 


Gap_ =~ | 2.87 | 2.02 | 1.46 | 1.11 98 79 64 
Chord 
I = fo) O12 024 030 055 060 082 104 


These values are plotted in Figure 18. 
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Figure 18. Munk’s Interference Factor for Induced Angle of Attack 


ONNTERFERENCE FACTOR-L 


Prandtl’s interference factors. Prandtl has shownS that the 


drag of one wing of a biplane in the presence of the other may be 
expressed in the form 


Lk; 


D=s en a (11) 


where D,, is the drag of wing 1 (having lift L, and span b,) in 
the presence of wing 2 (having lift L, and span 0,) and o is an 
“interference factor’’ which varies with the ratios of gap to span 
and shorter span to longer span. 


5‘*Techniche Berichte,’ Vol. II No. 6, (N.A.C.A. Technical Note No. 182) and 
N.A.C.A. Technical Report No. 116. 
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Prandtl gives values of ¢ in terms of b. and ——8°P 
b, average span 


2G ees 
or ea) as shown in Figure 19. These data have been con- 
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Figure 19. Prandtl’s Drag Interference Factor for Biplanes in Terms of 
Average Span 
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verted and replotted in Figure 20, using the ratio of gap to 
maximum span, instead of the ratio, gap to average span. 


Induced drag of biplanes. The interference factor ¢ may 
be used to calculate the value of the span factor k for any biplane, 
as shown by Prandtl. The method is as follows: 
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Let b, = span of longer wing 
b, = span of shorter wing 
y = — = ratio shorter span to longer 


b 

L, = lift on longer wing 
L, = lift on shorter wing 
L = L,;+ £,= total litt 


Prandtl’s equation for total induced drag is, 


2 


I 1 NES 
Dante ate oped (12) 


Maximum Span 
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This is a minimum when 


L, vrs 
= fas 
——o6 
p. 
and has the value 
oes yi I — o? 
Minimum Di xqb, Leroeand (13) 


Assuming the lift proportional to the area and setting 
S,=rS, it follows that S,=(1-7)S, Ly=rL, and L,=(1-r)L. The 
factor 7 is obviously the ratio of area (or lift) of the longer wing 
to the total area (or lift). Substituting these relations into 
equation (12) gives 


Ds = =| 2 Pee e meatal 


from which the span factor is seen to be 


ae yw? z 

r? (u? — 2uo-+1) +2r(yo—1) +1 (14) 

Figures 21, 22, 23, 24, and 25 give the values of & against u 

and r for = ,05, .10, .1§, .20 and .25. This covers the 

extreme range ordinarily used. The variation of k is substantially 

linear with . between any two adjacent values of ¢ and 

I Ir 

hence k may be obtained by interpolation. However, a great 

majority of biplane designs have wings which approximate either 

equal chords or equal aspect ratio. For this purpose Figures 

26 and 27 have been prepared. These give the value of k 
directly in terms of b. and & 
b iy be 

Proportions of the most efficient biplane. Figure 28 gives 

the proportions of the most efficient biplane as determined from 

Figures 21 to 25. From this diagram the best value of any one 

variable, Fe a or zo is determined when the other two are 

assumed or known. 
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Induced drag charts. The importance of obtaining a high 
Co: 
Ce 
effective aspect ratio as in Figure 29. The change in induced 
drag in passing from aspect ratio 1, to m, is 
DX Cpi = (Ke as K,) Cr = AW OF he 

The value of ACp; is readily found from the chart, Figure 30, 
in which AK is plotted against n, and n,. Note that , and 


ry kb, = k bs i 
n, are effective aspect ratios, 7.€., #1, = tees and 2, = oe . 


aspect ratio is clearly brought out if K = be plotted against 


Kb)* 
$s 


— 
2 
re 
< 
c 
re 
oO 
Ww 
2 
ul 
2 
F 
oO 
Ww 
L 
L 
uj 


Figure 29. Induced Drag Factor K = ©2# 


Cr 
Practical applications of induced drag. The division of wing 
drag into induced and section drags enables the value of a 
particular wing section to be appraised at a glance. This is 


particularly true if the section drag is plotted against lift 
coefficient. 
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Performance calculations are immensely simplified by the 
use of induced drag, as explained in Chapter VII. The effect of 
aspect ratio on performance is easily obtained. If the perform- 
ance of an airplane is known with one wing section, the per- 
formance with another wing section is obtained by correcting 
for the difference in profile drag. 
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Figure 30. Change in Induced Drag Factor with Change in Aspect Ratio 


Center of pressure. Munk has shown in a remarkable 
paper® that the classical treatment of wing theory by means of 
vortices may be replaced by kinetic energy considerations. 
In this paper he shows that the lift due to the curvature of a wing 
acts at 50% of the chord while lift due to angle of attack acts 
at 25% of the chord. Consequently the moment about a point 
on the chord 25% aft of the leading edge is constant. This 


6 Max M. Munk, ‘Elements of the Wing Section Theory and of the Wing Theory,” 
N.A.C.A. Technical Report No. 191 (1924). 
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relation has been amply verified by wind tunnel tests. It is of 
considerable value in the solution of a number of practical 
problems, enabling the center of pressure to be estimated with 
considerable accuracy. The center of pressure is readily found 
from the moment coefficient about the 25% point, by the 
relation 


= pee 
Cp = 0.25 Cy 
where Cy is the normal force coefficient (Cy = Cz cos @ + 


Cp sina). For all but very small lift coefficients and very 
large angles of attack C; may be used for Cy. 


Zero lift—zero moment. Munk has shown? that the angle 
of attack for zero lift is given approximately by the line drawn 
through the trailing edge and a point located on the mean 
camber at 50% of the chord. The exact angle of zero lift is 
determined as follows: pass a straight line AB through the 
trailing edge A and a point on the mean camber 11% of the chord 
aft of the leading edge, pass a second straight line AC through 
the trailing edge and a point on the mean camber 89% of the 
chord aft of the leading edge. The line AD bisecting the angle 
BAC is the direction of zero lift. These constructions are shown 
on Figure 31. 

The angle for zero moment about the leading edge is found 
by passing a straight line through two points on the mean 
camber. The first point is at 444% of the chord and the second 
is at 9534% of the chord measured from the leading edge. This 
construction is also shown on Figure 31. 

If the wing section has a large leading edge radius 7, the 
mean camber curve should pass through the center of this 
arc and be prolonged to a point P, which is one-half 7 distant 
from the center. The chord length should then be measured 
from the point P instead of the actual leading edge. 


7 Max M. Munk, ‘“‘The Determination of the Angles of Attack of Z Lift 
Moment, Based on Munk’s Integrals,’” N. A. C. A. Technical Note 122 (Kesar aus Zee 
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APPROXIMATE ZERO LIFT LINE 
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Figure 31. Munk’s Methods for Finding Zero Lift and Zero Moment Lines 


Wing with a fore-and-aft slot. Prandtl gives the theoretical 
solution® for the effect of a fore-and-aft slot in a wing, such as 
that existing at a panel joint. This solution is of great practical 
value in that it indicates the necessity for avoiding any kind 
of ‘‘leakage”’ joints in wing construction. 


8 Due to Grammel and Pohlhausen. See N.A.C.A. Technical Report No. 116. 
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Figure 32. Munk’s Span Factor for Wing with Fore-and-Aft Slot 
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Figure 33. Observed Increase in Drag Due to Fore-and-Aft Slot 
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For a monoplane with a fore-and-aft slot, Munk’s span 
factor k is no longer unity, but much less than unity, even for 
a quite narrow slot. The variation of k with the width of the 
slot is given on Figure 32. The curve marked A on this figure 
is an enlargement for small slot widths. A slot width of .oor 
X span reduces k& to 0.87, which is equivalent to a reduction of 
about 24% in aspect ratio. This effect has long been known 
from Munk and Cario’s wind tunnel tests.? 

The average increase in drag due to various slot widths, as 
found in these tests, is given in Figure 33. A slot width of 
1% of the chord on a wing of aspect ratio 6 increases the drag 
about 6%. 


Correction for wall interference. Prandtl has shown’® that 
the boundary of any finite wind stream, in either a closed or an 
open-working section wind tunnel, restricts the flow past any 
object under test. The boundary conditions of constant pres- 
sure for the free jet and zero normal velocity for the closed 
tube are satisfied analytically by the introduction of a series 
of images which give an induced velocity at the model under 
test. This induced velocity is equivalent to a downwash for 
an open section or an upwash for a closed section. ‘This causes 
the angle of attack and the drag, as measured, to be higher in a 
free jet and lower in a closed tube, than in an infinite air stream 
by the amounts 


57-3.Cr) 7.16 Crs 


(Aa)° me 2nD? A (15) 
PEG... 
De ay BA (16) 


Where SS is the model wing area, D the jet diameter, and A 
the jet area. The corrections must be added to the readings 
obtained in wind tunnels having an enclosed test section. The 
corrections are subtracted when the jet is of the open type. 


9 Technische Berichte I-6. 
10 G6éttingen Nachrichten 1919, also N. A. C. A. Technical Report No. 116 (1921). 
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Glauert has shown" that the corrections for a square jet 
are about 10% greater than for a circular jet. 


The general drag equation. The chief variables which affect 
the air force on any aerodynamic object, such as a wing, are: 
(x) angle of attack a, (2) relative wind velocity V, (3) scale of 
the object ZL, (4) density of the fluid », (5) compressibility of 
the fluid, (6) gravity g, (7) surface roughness, and (8) texture 
of the air flow. From the theory of dimensions it may be 
shown that the general drag equation, for example, is 


Eigen ePVLV VI » 
D = eViir| &  @EEYV | 


(17) 


We are accustomed to the use of a drag coefficient Cp instead 
of the function in the brackets, but in so doing, sight must not 
be lost of the fact that Cp may and does vary with a number 
of non-dimensional terms, and that this variation must in 
general be determined experimentally. The purpose of the 
ordinary wind tunnel test is to determine 9, («) or the variation 
of Cp with angle of attack. The purpose of testing at various 


oVL 


speeds is to determine 9, (ee over a part of its range. Propeller 


designers understand the effect of ,(*) on propeller character- 
a 


istics in general, if not in particular, when they try to hold the 
tip speeds well below the velocity of sound. Model tank ex- 


periments are usually made on the basis of a because this 
term brings in the gravity effects associated with wave making. 
If tests on two models do not agree, the surface roughness £ may 


not bear the same relation in the two cases, and this is often 
observed in skin friction tests. Finally, if two wind tunnels 
do not agree in tests on the same model, it may be due to a 


Irs The atc of Wind Ch 1 Wall; 
Be Ane oie es "i he spears alls on the Aerodynamic Characteristics of 
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difference in turbulence represented by the ratio of the average 
lateral turbulence velocity to the measured axial velocity. 
These terms do not have the same weight. In aeronauti- 
cal engineering, angle of attack is considerably more important 
than turbulence or surface effect. Next to angle of attack in 


importance is either the Reynolds Number Le or the com- 


pressibility effect Yr So long as the relative speed is low the 


compressibility effects are negligible, but as the speed increases 
they become increasingly more important. It is for this reason 
that wind tunnel tests at very high speed do not agree with 
tests at the same Reynolds Number at a lower speed. Such 


tests may include both 9, (=) and 93 (*). 


Determination of Reynolds Number. ae will have the 


same value in any consistent system of units. If the units for 


oe, V, and L are respectively Ibs. mass per cubic foot (e = 2) 


ft./sec., and ft., then » must have the dimensions of In 
the c. g.s. system the value of u is given for air by Holman’s 
formula 
107 @ = 1,715.5 (1 + .00275¢ — .00000034 #”) (18) 
In ft.-lb.-sec. units , is 
10?° wp = 3,582.9 (1 + .00275 # — .00000034 2”) (18a) 
t being in °C. 
The deviation of » from a straight-line function of ¢ is less 
than 1 part in 3,000 over the usual working range and one may 


write 
107°. = 3,408. + 5.483¢ (°F.) (19) 
= 3,583. + 9.870¢ (°C.) (19a) 


for use in determining the value of Bre. The value of y at 
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15° C. or 59°F. by the approximate formula is 3,731. X Io. 
while Holman’s formula gives 3,730. X Io. ne 34 has 


been prepared for use in reading the value of & ae or 


ile AIR TEMPERATURE %-°F SS 
30 40 50 60 70° 80 90.100 HON @ieG 
Figure 34. Reynolds Number Coefficient 


for checking calculations. Reynolds Number is obtained by 
multiplying the value of : from this figure by V in ft./sec. and 


Lin ft. The engineer will find it convenient to remember that 
in standard air a wing of 1.0 ft. chord at 100 mi./hr. gives a 
Reynolds Number of 925,000, or approximately one million. 


CHAPTER III 
AIRPLANE MODEL TESTS 


The purpose of a routine wind tunnel test on an airplane 
model is to determine performance, stability, and control, with 
reasonable engineering accuracy, and the advantage of a wind 
tunnel test over any other method lies chiefly in its quick and 
definite indications when reasonable deviations from mathe- 
matical exactitude are allowed. Routine testing of airplane 
models cannot well be justified on any other basis. 

It is an easily demonstrated, but not a widely recognized 
fact, that any model test is a compromise. There are always 
some conditions of similitude which conflict or cannot be met, 
but these conditions do not all have the same weight in determin- 
ing air forces. Since there is no known theoretical method of 
assigning the proper relative values to the various conditions, a 
practical solution based, for example, on experience or on trial 
and error must serve instead. The principle of dimensional 
homogeneity leads to an equation for air force in the form 


oVE Vive th oe a 


u. ? a gL’ fo V’ c (17a) 


F = pV?L?-9 | « 


where @ is the angle of attack, ue the well-known Reynolds 


2 
Number, Le the ‘‘wave making” or gravitational term, e the 
§ 


ratio of the relative wind to the velocity of sound, ra measure 


C 

the aspect ratio. These are by no means all of the dimensionless 

combinations which may be written, but they comprise what are 
53 


v : : b 
of surface roughness, 7? measure of turbulence in the air, and — 
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usually considered the most important ones. It is desired to 
point out three facts: (1) theoretically no model test can com- 
pletely represent a full-scale condition unless all of the dimen- 
sionless ratios are held constant; (2) it is impossible to hold all 
of these ratios constant at the full-scale value in a model test, 
for some of them are contradictory; and (3) the dimensionless 
ratios may not have equal weight but theory does not show it. 
Experience tells us, however, that the angle of attack @ is the 
most important of all the terms given, and that Reynolds 
Number probably comes second in importance. The “wave 


making” term 5 which is so important in testing ship models, 


is negligible in air work. The term r does not have any great 


importance at speeds less than about 300 ft./sec., where com- 
pressibility effects begin to appear. The surface roughness 
has a very great effect under many conditions and cannot be 
neglected. Turbulence is known to have a fairly large influence, 
particularly on critical flow conditions. The large effect of 
aspect ratio is well known, and so on through the list of minor 
unconsidered terms. 

The only major term which cannot be satisfactorily met is 
the Reynolds Number ee. If the value of the function was 
the same for all airplanes, or even for all parts of an airplane, the 
problem would be capable of a satisfactory solution. Actually 
the effect of Reynolds Number is different on each item; wings, 
fuselage, struts, wires, etc. Furthermore, the effect is not the 
same on all wings or on all struts but varies widely in each 
group. 

As a concrete example, the drag coefficient on streamline 
struts may be four times as great on a 1/24 scale model as on the 
full-size airplane; streamline wire may show even greater scale 
effects. This is very much greater than any scale effect on 
wings or fuselage. The conclusion is inevitable, that reasonable 
engineering accuracy is best secured by constructing a model 
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airplane with no streamline struts, wires, fittings, or minor de- 
tails, and correcting the measured model drag for the scaled-down, 
calculated full-scale drag of the omitted parts. The wings of 
the simplified model can be held together and the landing gear 
attached by means of a simple system of round brass wire 
struts (about 3/32” diameter), threaded right-hand on one end 
and left-hand on the other; the minimum length of wire necessary 
to insure rigidity being used. This method has three very 
important advantages. 

The drag of the round wire struts is quite definite and can 
be determined very accurately from runs with an exact dupli- 
cate set of struts attached symmetrically and in such a manner 
as to avoid interference while duplicating eachitem. If ADbe 
the drag of one set of struts, it is obvious that 


D, = (D. + AD) — [(D.+2AD) — (Do + AD)| (20) 


where D, is the drag of the model without struts, wires, fittings, 
or miscellaneous omitted parts. The full-scale drag of these 
parts can be calculated with considerable accuracy. Assume 
that the value so calculated at a speed V, be d. Then assuming 
d to vary as (VL)? 
 aS2 ay 

dm = dS ( a (21) 
where § is the (fractional) model scale ratio (for example, 1/20), 
V the model test speed, and d» the correction to be applied to 
_ the model test. The correct model drag will be 


D= Ds ie dm (22) 


Lift and moment corrections are also obtained from the 
“duplicate strut” runs. An additional moment correction due 
to the omitted parts is obtained by calculating the line of action 
of d (and therefore d,,) so that the distance of this line of action 
from the moment axis is known. The table shown in Figure 35 
is a very convenient form for calculating simultaneously the 
drag correction and its line of action. A typical lift wire calcula- 
tion is given to illustrate the steps followed. 
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The other advantages of the method are that it allows very 
accurate alignment of the model, and reduces greatly the cost 
of model construction. The necessity for very accurate model 
alignment is not generally appreciated, but it has been shown 
at the Washington Navy Yard that any difficulty in checking 
test data after a lapse of time is almost entirely due to warping 
or change in alignment of the model. Of course, this difficulty; 
is partially obviated when metal wings are used. 
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Figure 35. Tabulation of Drag Correction for Airplane Model Test 


The simplified model cannot allow for scale effect on wings or 
fuselage, but in general it does confine the total scale effect to 
these two items. The full-scale values of profile drag for various 
wing sections are given in Chapters I and IV and an estimate 
can be made of the scale effect on the wings if greater accuracy 


is required. The scale effect on fuselage drag is not normally 
large. 


Ch. 3] AIRPLANE MODEL TESTS 57 


To sum up briefly, it may be stated that in order to obtain 
good results in an airplane model test it is necessary to omit 
all minor parts which show excessive scale effects, such as 
struts, wires, fittings, small attachments, etc. A drag correc- 
tion which includes the resistance of all omitted parts must be 
calculated and added to the basic model drag. 

The objections to this method are usually based on the two 
arguments: (1) that the principle of dynamic similarity is 
violated, and (2) that interference effects are neglected. These 
objections cannot be maintained since it has been shown that 
the principle of dynamic similarity does not apply for an exact 
scale model and since the interference effects of a strut, for 
example, certainly cannot be more than a small fraction of the 
four-fold resistance which this method corrects. On the basis 
of results, the method is quite satisfactory. The high speed 
predicted by a careful wind tunnel test on a simplified model 
normally differs less than about 3% from that obtained in 
flight tests, and it is not unusual to obtain less than 1% difference. 

For additional detail on model construction and test correc- 
tions see N. A.C. A. Technical Notes No. 82, ‘‘Notes on the 
Construction and Testing of Model Airplanes,” W. S. Diehl 
(1922), and No. 254, ‘Method of Correcting Wind Tunnel 
Data for Omitted Parts of Airplane Models,” R. H. Smith (1927). 

Pitching moments. Information concerning the longitudinal 
static stability and control is given by pitching moment curves 
for three or more elevator settings at a selected stabilizer setting. 
This part of the tests and the lift and drag data should be 
considered of equal importance. 

In general, the pitching moments will be measured about 
an axis which is parallel to, but considerably displaced from, 
the gravity axis. Since the final moments must be referred to 
the c.c., either a graphical or an analytical method must be 
used in the conversion. ‘The graphical method is simple, quick, 
and accurate if care is used in drawing the vector diagram. 
This diagram once drawn is always available for rapid calcula- 
tions of the pitching moment about any new c.c. While an 
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experienced engineer can obtain an excellent idea of stability 
characteristics directly from the vector diagram, the slope of 
the moment curve is the only reliable criterion on which to 
base definite conclusions regarding static stability. 

All airplane model tests should include pitching moment 
curves about the center of gravity for at least two different 
stabilizer settings with elevators at o°. The purpose of this © 
being to determine the moment change due to a definite incre- 
ment in stabilizer angle, so that settings required for balance, 
or moments for new conditions, may be accurately calculated 
from the test data. It is easily shown that changing the stabil- 
izer setting or élevator angle does not appreciably affect the 
slope of the moment curve, and merely shifts it up or down by a 
definite increment. 

The elevator angles usually investigated are 0°, 5°, +10°, 
—15°, —20°, and —30° where the + sign refers to “down” 
elevator. With the exception of the o° and +10° settings, 
it is unnecessary to cover a greater angular range than that 
required to give the zero moment, or balance, condition. In 
general, the moment curve is not satisfactory unless it has 
a negative slope over the entire range of flight angles. Slight 
waves in the curve giving zero or slight positive slopes are 
not actually dangerous, but they are indicative of some very 
undesirable condition. The desirable slope of the moment curve 
varies with the type and purpose of the airplane. The slope 
of the full-scale moment curve should be either 


dM 
or 
dM 
Pere ao (24) 


dM . : 
where aa 8 the slope (in Ibs.-ft. per degree), g the dynamic 


pressure 3 ¢V*, W the gross weight (Ibs.), C the wing chord (ft.), 
and S the total wing area (sq. ft.). The value of K should lie 
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between .coo6 and .ooro according to the stability desired. 
The corresponding values for K, are .co6 and .o1o, but K, is 
not necessarily equal to 10K. Some additional data on slope 
of the pitching moment curve may be found in Chapter V. 

Whether or not the elevators give adequate control is usually 
quite obvious from an inspection of the moment curves or 
vector diagrams. It is a common rule to assume that 20° up 
elevator should balance the airplane at the angle of attack 
where maximum lift is obtained, assuming that balance with 
neutral elevators is in the normal range between o° to 6° angle 
of attack. 


Rolling moments. Rolling moments need not be measured 
unless there is considerable doubt in regard to lateral control. 
Under these conditions it is advisable to measure both rolling 
moments and yawing moments due to the ailerons. The rolling 
moment due to yaw is not often measured except in. stability 
investigations. 

Experience has shown that satisfactory lateral control is 
assured when the full-scale aileron rolling moment is equal 
to or greater than 


L = .0060 qWb (25) 


for 10° setting at o° angle of attack. W é is the gross weight, q¢ 
the dynamic pressure, and 6 the maximum span. This is 
based on satisfactory airplanes, and will assure adequate rolling 
moments up to stalling angles for a normal design. 


Yawing moments. The conventional tests to determine 
directional stability and control consist of three parts as follows: 


Yawing moments are measured with fixed rudder setting 
and angle of attack. Usual range +20° in yaw. 

Yawing moments due to various rudder settings are 
measured with body at 0° yaw and fixed pitch. Usual 
range of rudder angles 0° to 20°. 

Angles of yaw held by various rudder settings up to 20°. 
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While it is desirable that these tests be made at several 
angles of attack, the time required is often an important factor. 
If the stability and control shown by the first tests are fully 
satisfactory, little is to be gained by tests at other angles of 
pitch unless the lateral stability is being investigated completely. 
The chief criterions for a satisfactory condition in yaw are: 


Definite slopes to the yawing moment curves. Absence 
of either irregularities in the curves or unsteadiness in 
test readings on balance, particularly at small angles of — 
yaw. 

The slope of the curve of angle of static yaw held by the 
rudder when plotted against rudder angle must be 
definite, and the curve must be free from irregularities. 
The ratio of static yaw angle to rudder angle should be © 
greater than 0.5 and less than 1.0. A value less than 
0.5 may indicate either too much directional (weather- 
cock) stability or deficient rudder control. A value 
greater than 1.0 may indicate either too little directional 
stability or too much rudder control. 


s 


The slope to the curve of yawing moments with rudder 
neutral may be written in the form 
dN 


ay = KaWe ess 


in which i is the slope (Ibs.-ft. per degree), W the gross oat 


(Ibs.), and 6 the span (ft.). An analysis of 61 airplane model 
tests in comparison with reported full scale characteristics shows 
that K should be between .cooo40 and .ooco80 for a normal 
airplane. Values slightly lower or higher are not entirely 
objectionable, but the most favorable comments have been 
made on airplanes for which K is within the limits given. 


Calculation of performance from model tests. Assuming 
that the air forces vary as (V1), the relations between model and 
full scale are readily obtained. Designating the model values by: 
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the subscript m and letting the scale ratio pe EN =e bm the 
Full scale 1 


drag relations at a given angle of attack are 
Model drag n= (Co) m§ (Vin lin)? (27) 
Full-scale drag D= Cot (V1)? (28) 


Since the air forces are assumed to vary as (V/)?, the drag co- 
efficients will be equal Cp = (Cp)m. Hence 


D De 
(VI)? . (Vin bm)? 
EVA (VN? 
or D= Dn (*) (7) (29) 
_ Similarly L=Lnm cae es) (30) 
m! \Vin 
In steady horizontal flight L = W. Hence 
ie E Ww 
V-= Va X ( ; a (31) 


The minimum, or stalling, speed is given by substituting the 
maximum model lift Ly» mar in equation (31) 


Fi In) eae 
Vee Vik @ \ ae (32) 


-It*‘is often more convenient to work with speed ratios, rather 
than to calculate each speed separately, 


V ae. \ os max = “|= max 
ig i Te or V Se Vs 5 (33) 


but equation (31) can be simplified for any given model, to 


fins AS oe | Ga) 


where K = W(Vm) : x (‘") Equation (34) can be solved very 
| : : if 


ra , 
Ze © (/ fUc el 
’ 


a 


X 
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rapidly on a slide rule by setting the runner to K on the A scale 
and moving the slide so that L,, on the B scale is in line with K. 
Then the index on the C scale will rest at the desired value of V 
on the D scale. 

Full-scale drag values may be calculated by equation (29) 
but time can be saved by using the values of L/D from the 
model test. Where curves of power required are to be calculated 
for two or more weights the values of L/D can be plotted against 


7 This curve can then be used in conjunction with equation 


Ss 


(33) since D = OB when L = W. 

In order to secure accurate results, particularly at the lower 
speeds corresponding to large angles of attack, it is advisable 
to correct both lift and drag for the elevator angle necessary to 
maintain balance. If lift, drag, and pitching moments are 
measured for several elevator settings, this may be done with 
sufficient accuracy by simple interpolation. Elevator angles of 
— 10° (down), o° (neutral), —5° (up), — 10°, —15°, and —20° 
are usually required, but the —5°, —15°, and — 20° settings 
need be investigated over only a limited angular range at high 
angles. 

The methods employed in a typical performance calculation 
will now be illustrated by an example. Table 1 contains the 
wind tunnel data obtained in the Washington Navy Yard 
tests on a 1:16 scale model of the DH-4B, the lift and drag 
values, L and D, having been corrected for elevator setting. 

The maximum lift is 7.571 Ibs. (at 20°) at a speed of 4o 
mi./hr. The gross weight is 3,876 lbs. Hence, according 
to equation (32), the stalling speed is 


_ 40 [3,876 _ r 
Te Nosh 56.6 mi,/hr. 
For any other weight, the stalling speed is 
40 W 
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A 
1 = = 
also Vs = Vs, W, 
or, the speed at any lift is 
Le GMS 
‘es 16 ibe 


Table 2 contains the calculations for power required in 
horizontal flight at sea-level with two gross weights. The 


So 60 70 80 90 100 ito 120 130 
Figure 36. DDH-4B Power Curves from Wind Tunnel Tests 
corresponding powers required at any desired altitude may be 
obtained by multiplying each value of V and 7.Hp., by the 


value of \ = for the altitude under consideration, as explained 


in Chapter VII. Figure 36, is a plot of the 7.Hp. available and 
T.Hp. required against air speed at sea-level. The performance 
indicated is in excellent agreement with the flight test results on 
this airplane in the condition represented by the model. 
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TABLE 1. WinpD TUNNEL Test Data on DH-4B AtrpLANE MODEL 
ScALE Ratio 1/16. TEST SPEED 40 mi./hr. (standard air) 


Angle Elevator ° 
a” | vanale | alto leaceeeeg | 7 aire ae 

Attack for ro ie D Vs 
a Balance iy a 
I a7 1.214 .392 3.10 2.500 
2 -4 1.853 . 404 4.58 2.020 
3 fe) 2.475 -420 5.90 1.746 
4 — .4 3.040 -458 6.63 1.576 
6 —1.2 4.053 -527 7.68 1.366 
8 2.0 5.045 .623 8.10 T2223 
10 3.2 5.858 -747 F303 1.136 
12 a7, 6.574 .903 GP az 1.073 
16 10.0 7.498 1.745 4.29 1.004 
20 18.5 757i 2.641 310 1.000 


TABLE 2. CALCULATION FOR T.Hp. REQUIRED IN HORIZONTAL FLIGHT 
AT SEA-LEVEL. DH-4B ArrpLane. 


(Based on Wind Tunnel Data in Table 1) 


Wind Tunnel 


Pata W = 3876 W = 4300 
Angle 
* of ; L re ot Drag Air Drag 
ttac a atio peed Speed 
a D V V He T.Hp., | V D Ded BG D>. 
Vs ||mi-/or.| “PS mi./hr.| Ibs 
I 3.10 | 2.500 |} 141.4 | 1,250 | 472.0 || 149.0 | 1,385 | 550.0 
2 4.58 | 2.020 |] 114.3 847 | 258.0 || 120.3 938 | 301.0 
3 5.90 | 1.746 98.8 657 | 173.3 |} 104.0 728 | 202.0 
4 6.63 | 1.576 89.2 585 | 139.0 || 93.8 648 | 162.2 
6 7.68 | 1.366 ee 505 | 104.0 81.3 560 | 121.3 
8 8.10 |) 55223 69.2 478 88.2 72.8 530 | 103.0 
10 7103 SO! 64.2 495 85.0 67.7 548 98.8 
12 Tee Th O73 60.7 533 86.3 64.0 591 | 100.7 
16 4.29 | 1.004 56.8 903 | 136.8 59.8 | 1,002 | 160.0 
20 3.10: | Ls000 56.6 | 1,250 | 189.0 59.6 | 1,385 | 220.0 
T.Hp.r = VD 


CHAPTER IV 
PARASITE DRAG DATA 


The engineering drag coefficient K is used frequently in this 
chapter. K is defined by the drag equation - 


Drag = K (2) AV? 


where A is the cross-sectional area in sq. ft., and V the air speed 
in mi./hr. K is therefore the drag per sq. ft. at 1 mi./hr., and 
is found from 


K = Cp® (Fy)? 


where (0 = .002378 and Fy is the conversion factor from mi./hr. 


to ft./sec. Fy = = (Fy)? = 2.151. Therefore K = .002558Cp 


and Cp = 391 K. 


Square flat plates. Eiffel investigated the drag of square, flat 
plates over a limited range in Reynolds Number * and the author 
is unable to find any more extensive series of tests. Eiffel’s 
results were originally reported as the effect of area on drag 
square plates as shown on Figure 37. These data have been 
converted and for a test speed of 10 m./sec. are as follows: 


Length of Side Reynolds 
Cm. Cp K Number 
10.0 1.040 . 00266 68,000 
15.0 1.055 . 00270 102,000 
25.0 1.072 .00274 170,000 
B75 1.140 . 00292 255,000 
50.0 1.193 . 00305 340,000 
70.7 1.234 . 00316 482,000 
100.0 1.263 - 00323 680,000 


I‘ The Resistance of the Air and Aviation,’’ Chapter IT. 
65 
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Eiffel states that for a very large plate Cp approaches the 
value of 1.28, or K = .00328. This value is generally used in 
calculating the drag at very high Reynolds Numbers, and in 
calculating the equivalent flat plate area of parasite drag. 


Rectangular flat plates normal to wind. The drag coefficient 
of rectangular flat plates increases as the aspect ratio is increased. 
Eiffel’s tests? are well known: Foppl made some similar tests 


a ee 
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Pakaeese: 
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Zoe 
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Figure 37. Effect of Area on Drag Coefficient for Square Flat Plates 


at Gottingen.* Eiffel’s and Foppl’s data do not agree very well, 
probably on account of the differences in Reynolds Numbers. 
Wieselsberger* has recently made some tests which appear more 
consistent; his data are plotted on Figure 38. He also tested 
a number of annular discs, the data for which are plotted on 
the same figure. It is interesting to note that the drag coefficient 
for infinite aspect ratio is Cp = 2.00, or K = .00512. 


Circular discs. The resistance of circular discs has been 
measured in various wind tunnels, and the results are in excel- 
lent agreement. Eiffel found Cp = 1.055 and 1.08 for 15 cm. 


2 Chapter II, ‘‘ Resistance of the Air and Aviation.”’ 


_3 Reported in the ‘Jahrbuch der Motorluftschiff—Studiengesellschaft,” 1910-11 
Julius Springer, Berlin. i 


4 Gottingen Ergebnisse II (1923), R. Oldenbourg, Miinchen. 
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and 30 cm. discs respectively, when tested at a speed of 10 m. /sec. 
Wieselsberger® tested four discs 3, 6, 15.1, and 4o cm. in diam- 
eter over a wide range in Reynolds Number obtaining values of 
Cp varying from 1.07 to 1.13, with an average value of about 1.11. 
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ABSOLUTE DRAG COEFFICIENT-CD 


BREAOTH any INNER  OIAMETER 
LENGTH OUTER DIAMETER 
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Figure 38. Drag Coefficients for Rectangles and Annular Discs 
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Figure 39. Drag Coefficients for Spheroids and Ellipsoids 


A 15 cm. disc tested in the N. A. C. A. variable density wind 
tunnel® at Reynolds Numbers ranging from 210,000 to 4,440,000 


5 See “Ergebnisse der Aerodynamischen Versuch-Anstalt Zu Géttingen” II (1923), 
Pp. sg R. Oldenbourg, Miinchen. 
6 James M. Shoemaker, ‘Resistance of a Fifteen-Centimeter Disc,” N. A.C. A. 
Technical Note No. 252 (1926). 
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gave values of Cp between 1.077 to 1.139 with an average value 
of about 1.11. 

In both the Gottingen and the N. A. C. A. tests there was no 
definite change of Cp with Reynolds Number, although there 
is a little evidence of the existence of the critical Reynolds 


0.8 
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04 

0.3 ae a 
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0.1 


100" 
Figure 40. Variation of Cp with ai Notte eden ye Spheres 


Number which is found for cylinders and spheres. For any 
ordinary size of disc or wind speed the average value Cp = 1.11 


may be used. For Reynolds Numbers above 6,000,000 the coeffi- 
cient is unknown. 


Spheres. The drag of spheres has been measured over an 
extremely wide range in Reynolds Number. The variation of 
Cp with Reynolds Number given on Figure 4o is based on a 
compilation of test data from various sources made by Dr. A. F. 
Zahm? and Mr. F. A. Louden of the Washington Navy Yard. 


Spheroids and ellipsoids. Drag coefficients of two ellipsoids 
are compared with those for a sphere® in Figure 39. 


iA. F. Zahm, ‘‘Flow and Drag Formulas for Simple Quadrics,’’ N. A. C. A. Technical 
Report No. 253 (1926 ys 


8 Based on data in ‘“‘Géttingen Ergebnisse” IT. 
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Hemispheres—parachutes. LEiffel tested a hemispherical cup 
and found Cp = 0.33 when the convex side was presented to 
the wind, and Cp = 1.33 when the concave side was presented 
to the wind. The latter value is used for parachutes, on the 
basis of projected area when inflated. The form of an inflated 
parachute is approximately a part of the surface of an ellipsoid 
and the projected area is about 55% of the surface area. The 
limiting velocity of a parachute is therefore given by 


2Ww 


ieee DES 
5 1.33e (0.55 9) 


(35) 


where W is the total weight (lbs.) of parachute plus its load, 
S the total surface area (sq. ft.), and o the air density. For 
standard density the limiting velocity in ft./sec. is 


Vi= 34 Ve (36) 


Cylinder with axis normal to wind. The drag of cylinders 
normal to the wind has been thoroughly investigated. Figure 4o 
gives the variation of Cp with Reynolds Number both for a 
section of an infinite cylinder and for a section 5 diameters long. 
These curves are based on test data obtained at the Washington 
Navy Yard and other laboratories, as compiled by Dr. A. F. 
Zahm? and Mr. F. A. Louden. 


Cylinder with axis parallel to wind. Eiffel tested two series 
of cylinders with the axes parallel to the wind. In one series 
the diameter was 30 cm. and in the other 15 cm. Lengths from 
a thin disc up to 7.0 diameters were tested at a wind speed of 
to m./sec. The ratio of the drag of the cylinders to the drag of 
the discs are plotted against length/diameter ratio in Figure 41. 

Subsequent tests on the cylinder with length/diameter = 7.0 
showed the value of Cp to be reduced from 0.94 to 0.19 by fitting 
hemispherical ends. 


9See N. A.C. A. Technical Report No. 253. 
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Figure 41. Drag Coefficients for Cylinders Parallel to Wind 


Skin friction. The skin friction resistance of a very smooth 
surface such as glass, varnished wood, or doped fabric, is given by 


Boece SV (37) 


where ¢ is the air density, S the fotal area or “‘wetted surface,” 
V the air speed, and Cpr the frictional drag coefficient given on 
Figure 42. Cp, is also given by the relation 


—O0l5 
= .0375 (R N)-™'s (38) 
hence 
Dr = .0375 (RN)-o1s* £. Sy (39) 


The frictional drag coefficient of a surface which is not 
very smooth, such as unvarnished wood, undoped cloth, etc., 
varies with the roughness of the surface between .oo5 and .o10 
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t 5 2 vg 4 5 6 8 40 
Figure 42. Skin Friction Drag Coefficient Cor 


and is substantially constant for all Reynolds Numbers over 
a wide range including the common full-scale values. 

Some doubt exists as to the correct form of equation (39) 
on the basis that breadth should enter directly, so that the 
frictional coefficient varies as a function of length and Reynolds 
Number. There is also some doubt as to the validity of the 
equation for short lengths, but this doubt probably arises from 
the difficulty which some experimenters have had in separating 
form resistance from frictional resistance. The total “profile 
drag”? of a thin section may be less than the skin friction. 
Equation (39) is certainly sufficiently accurate for most engineer- 
ing purposes. 


Wing section drag. Wing section, or profile drag is largely a 
function of camber and thickness. If the full-scale values of 
minimum drag coefficients obtained in the N. A. C. A. variable 
_ density wind tunnel be plotted against thickness, as in Figure 43, 
it is found that the variation is linear over the range investigated 
and the straight line extended to zero thickness intercepts the 
drag axis at Cp, = .00375, which is just half of the skin friction 
coefficient (for the double surface). Actually the minimum 
section drag coefficients may be less than the skin friction 
coefficients for moderate cambers, but it is unlikely that any 
observed value would ever go below, say Cp = .oo6o. 
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Struts. Figure 44 shows the drag coefficients. obtained at 
various speeds, on a series of four geometrically similar Navy 
No. 1 struts of fineness ratio 3.0, tested at the Washington Navy 
Yard. This strut*® has the lowest drag coefficients yet recorded 
for a strut section. The offsets are given in Table 3. 

Drag, in lbs./ft. at 100 mi./hr. in standard air for a strut of 
fineness ratio 3, is plotted against strut thickness in Figure 45. 
These values should probably be increased about 20% for the 
commercial shapes. Rounding the trailing edge does not appear 
likely to give large increase in drag, providing the radius is small. 
Cutting off about 15% from the rear of a strut does not affect 
the drag at all, according to Munk’s tests reported in the Tech- 
nische Berichte. 

The change in drag with fineness ratio is shown on Figure 46. 
This is a composite curve based on British, German and American 
tests, which differ very little. 


10 A. F. Zahm, R. H. Smith, and G. C. Hill, ‘Point Drag and Total Drag of Navy Struts 
No. 1 Modified,’ ’'N. A. C. A. Technical Report No. 137 (1922). 
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Figure 44. Drag Coefficient Cp as a Function of Vd for Navy Struts 


TABLE 3. 


OFFSETS FOR NAvy STRUT AND C-CLAss AIRSHIP 
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Navy C-Class 
No. I Air- 
Strut ship 
100.0 99.9 
99.5 99.0 
95.0 95.0 
86.1 88.5 
73-2 79-9 
56.2 66.5 
33.8 49.3 
19.0 36.2 
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Figure 45. Drag per Foot at yoo mi./hr. Navy Struts 
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Figure 46. Variation of Drag Coefficient with Fineness Ratio for Struts and 
Streamline Bodies 
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Aircraft cable. The variation of Cp with Reynolds Number 
for round wire is given on Figure 47. The range in Cp actually 
encountered is not very large and it is customary to take Cp = 
1.25 for all normal sizes and speeds. This applies only to full 
scale, and gives 


, 
D lb/ft. = 2.56d (+) (is) 


_ where d is the diameter in inches and V the air speed in mi./hr. 
at standard density. The drag of various sizes of wires and 
turnbuckles may be found with other data in Table 4. 


TABLE 4. DATA ON AIRCRAFT CABLE 


Diameter Resistance at 100 mi./hr. 
Approxi- Breakin 
g 
Nomi mate ae Strength Wire Turn- Turn- 
i : oe min. Alone buckles buckles 
ee a ibs Ibs. per ft. | Plus Eyes} Plus One 


run Complete | Eye Only 


a .062 . 81 480 .16 .42 .25 
er .078 295 550 .20 .48 .29 
Pe - 094 1.45 920 24 -§2 - 33 
i .125 2.45 1,350 132 .65 .43 
as . 156 4.67 2,600 .40 .82 -53 
3 . 187 5.80 3,200 .48 1.01 W563 
a5 .218 8.30 4,600 .56 1.24 TA: 
4 250 10.50 5,800 64 1.50 . 86 
i 312 16.70 9,200 80 2.05 1.08 
3 375 5 96 2.60 133 


Figure 48 is a plot of total drag divided by total length, that is, 
drag per foot, at 100 mi./hr. for various sizes and lengths of cable, 
complete with terminals. 

The effect of fore-and-aft spacing is shown on Figure 409. 

The effect of inclination is given on Figure 50. The curves ; 
marked A and B are two sets of test data for which the drag ratio 
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Figure 48. Drag per Foot at 100 mi./hr. for Cables with Turnbuckles and Eyes 
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TWO WIRES 


RESISTANCE OF SINGLE WIRE 


Ri 
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1.20 


Figure 50. Effect of Inclination on the Drag of Wire and Cable. Curves A 
and B based on Unit Length, Curve C on Projected Length 
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refers to unit length. The curve marked C is the average drag 
ratio per unit of projected length. 

Lift and drag coefficients as a function of angle of attack are 
given for both smooth wire and cable on Figure 51. These are 
the usual absolute coefficients, for example, 


Litt ESV? 


S being the product of the wire length in feet by its diameter in 
Peet. 


Q° 20° 40° 60° so° tao° 
Figure 51. Lift and Drag Coefficients for Round Wire and Cable 
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Streamline wire. The drag of streamline wire has been 
measured by N. P. L.%" and N. A. C. A. and the variation of 
Cp with Reynolds Number is given on Figure 47. The two sets 
of data fall on a single curve. 

Figure 52 gives the drag in lbs./ft. at 100 mi./hr. and also the 
drag of two standard terminals at the same speed. It is custom- 
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Figure 52. Drag of Streamline Wire and Standard Terminals 


calculating the drag. This is an approximation which is 
sufficiently exact if the angle is not less than about 40°. 

Figure 53 gives the total drag divided by total projected 
length for various sizes and lengths of streamline wires with 
standard terminals. 

The effect of fore-and-aft spacing on double streamline wires 
is shown on Figure 49, and the effect of yaw on the drag of a 
single wire is shown on Figure 54. 

Table 5 contains data on the standard sizes of streamline wire. 


IR. & M. 256. N.P.L. is abbreviation for National Physical Laboratory (British). 
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Figure 53. Drag per Foot at 100 mi./hr. for Streamline Wire with Standard 
Terminals 
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Figure 54. Effect of Yaw on Drag of Streamline Wire 
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; Figure 55a. Fuselage Models 
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Figure 55b. 


Fuselage Models 


Figure 55c. Fuselage Models 


Figure 55d. Fuselage Models 


Figure 56a. Float Models 


Figure 56b. Float Models 


Figure 56c. Float Models 


Figure 56d. Float Models 


Figure 56e. Float Models 
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Fittings. The average drag of a fitting mounted on a wing, 
fuselage, float, or tail surface is given approximately by consider- 
ing the fitting as a flat plate of double its projected area. This 
allowance is necessary in order to account for the “‘interference”’ 
drag due to the fitting. The drag of a fitting is surprisingly high. 
The drag of an average strut end fitting with wire attachments on 
a 4,000-lb. airplane will probably be about 3.0 lbs. at 100 mi./hr. 
Most of this is in the wire terminals and can be included in the 
wire drag where it belongs. 

Miscellaneous projecting parts, not streamlined, should be 
treated as fittings. 


Control horns. The value of K for a well-streamlined control 
horn is about .coo8o, but if bolts and nuts are exposed K = .oo15 
is probably a fair average. 


Fuselage drag. Drag data on fuselages and similar airplane 
parts have been compiled by the Bureau of Aeronautics, Navy 
Department, and published as N. A. C. A. Technical Report 
No. 236.? Since the publication of this report additional tests 
have been made at the Washington Navy Yard on the models 
illustrated in Figure 55. ‘The data on these models are given in 
Table 6. 

Models Nos. 3, 4, and 5 are practically pure streamlines, 
which were tested in order to supply a basis for comparison. 
Models 7 and 9g were tested in order to find some form of simple 
model construction that would have the same drag as an air- 
cooled engine; No. 9 gives a close approximation. Several 
additional models not shown were tested in this investigation; 
the conclusion was reached that the model cylinders should be 
solid to the outer edge of the fins and up to the top of the cylinder. 
The remaining models closely represent well-known airplane 
fuselages. 

If all the engine, radiators, windshields, fittings, and other 
details are considered separately, the average value of K is 


12 W.S. Diehl, ‘‘ Tests on Airplane Fuselages, Floats, and Hulls,’’ N. A. C. A. Technical 
Report No. 236 (1926). 
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approximately .ooo4o. If all of the minor details are not con- 
sidered, the average value of K will be about .oooso. The latter 
figure takes care of the minor items only, such as fittings, and 
irregularities of construction. 

The drag coefficients for nacelles are not appreciably different 
from those for fuselages. 


TABLE 6. DATA ON AIRPLANE FUSELAGE MODELS SHOWN IN FIGURE 55 


Dimensions 
M Cross Measured Absolute 
is Gece Section Drag at K Drag 
engt Atea 40 mi./hr. Coefficient 
si lb 
A . 
sq. ft 

I 1.670 0457 . 066 00091 353 

2 1.668 0459 .050 00068 266 

3 1.669 0472 O12 Ooo16 062 

4 1.668 0548 .O16 00018 o71 

5 1.673 0621 .O16 00016 063 

6 1.673 0463 .022 00030 116 

Uf 1.673 0708 058 00051 200 

8 1.663 0550 021 00024 094. 

9 1.662 0662 068 00064. 251 
10 1.667 . 0862 .036 .00026 . 102 
II 1.667 .0481 .026 - 00034 . 132 
12 1.667 .0397 .027 . 00042 . 166 


Effect of pitch on fuselage drag. The drag of a fuselage with 
well-rounded section does not change greatly with angle of pitch 
while the drag of a fuselage with square or rectangular sections 
increases rapidly with increasing pitch. Consider models Nos. 2 
and 3 in Figure 55. 

The measured model drag in lbs. at 40 mi./hr. varies with 
pitch as follows: 


Angle of Pitch.... °° ce 10° 15° 20° 


Model No. 2...... .050 .057 .070 .083 .129 
Model No. 3...... .O12 .O12 .O10 .O14 .022 
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Rounding the sections of a fuselage reduces the drag at 0°, 
and gives an even greater improvement at cruising and climbing 
speeds. 


Floats. Table 7 gives data for a number of seaplane 
floats shown on Figure 56. These models, representing 
about an equal number of actual floats and of design studies, 
were tested at the Washington Navy Yard. The actual 
floats represented by models Nos. 1-10, and No. 19, show 
fairly uniform drag coefficients. For the average float K= 
.0005 5. 

Float No. 12 is practically a pure streamline form tested for 
comparison. Models Nos. 11, 12, 13, and 14 are included for 
comparison purposes only. They cannot be used as actual 
floats on account of undesirable water performance. 

The drag of a float may be calculated from the “shape 
coefficient’”’ Cs, the gross weight of the seaplane, and the excess 
buoyancy. Since the weight of sea water is 64 lbs./cu. ft., the 


cubic displacement for a weight W is (=). Letting the ratio 


of submerged displacement to load displacement equal (1 + e), 


the total volume of the float is sees and the drag 


ye? e+) Vate il 
D=csvr{ eA | @D 


The factor e is the decimal value of the usual percentage excess 
displacement. For example, with 80% excess, the value of ¢ is 
0.80. 

The use of a pronounced hollow bottom of the Short-Macchi 
type increases the drag 10% to 15% over that for the conven- 
tional ‘‘Vee” bottom, according to wind tunnel test data 
now available. However, a moderate amount of hollowness 
does not appear to affect the drag appreciably. 

Wind tunnel tests have been made to determine the effect 
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of astep on float drag. In general there is no measurable change 
in drag when a fairing is fitted behind the step. Test data 
indicate that a good streamline body fitted with hydrovanes 
has more drag than a conventional float even when due allow- 
ance is made for scale effect. 

The drag of an average float is comparatively low, and the 
drag of two floats may be slightly less than the drag of the two 
wheels in the corresponding land-type landing gear. The float 
bracing has more drag, however, than the land-type bracing. 
For an airplane of 3,000 lbs. gross weight, the parasite area for 
the seaplane arrangement is normally 1.0 to 1.5 sq. ft. greater 
than for the landplane arrangement. . 


TABLE 7. DATA ON SEAPLANE FLOAT MODELS SHOWN IN FIGURE 56 


é : | Absolute 
Dimensions Coefficients 
oe Drag :. 
Mode Cees at 40 odel 
No. Length Section eee x 
ft. Area A = Cs 
sq. ft. 

I 1.670 .0272 .0252 .00058 . 226 .060 

2 1.664 .0258 .0253 . OOO6I 2239 . 063 

3 1.667 .0276 . 0236 .00053 - 209 .058 

4 11670 . 0304 .0268 .00055 .215 - O61 

5 1.670 .0274 .0239 .00055 7213 .059 

6 1.667 .0269 .0228 . 00053 1207. .056 

7 1.667 .0340 .0286 .00053 .205 .063 

8 1.667 .0300 . 0339 . 00071 .276 .076 

9 1.659 .0292 .0243 . 00052 2203 .054 
10 1.655 .0305 .0178 . 00036 Sig) . O41 
II 1.665 .0365 .0139 . 00024 -093 .029 
12 1.667 . 0369 .0085 . 00014 .056 .020 
13 1.666 .0368 .O118 . 00020 .079 .026 
14 1.662 .0370 .O105 . 00018 .070 .025 
15 1.665 . 0339 0184 . 00034. 132 -O41 
16 1.665 .0253 .0247 . OOO6I .238 . 062 
17 1.654 .0405 .0313 .00048 . 189 .058 
18 E733 . 0319 . 0163 . 00032 .124 .034 
19 1.605 .0305 .O165 . 00034 . 132 .038 
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Wing tip floats. The drag of wing tip floats varies between 
_ .00040 AV? and .oco80 AV?. In the absence of specific data, the 
average value of K = .ooo60 may be used. 

Flying boat hulls. The drag of the F5L hull is approxi- 
mately .coo52 AV? and the drag of the NC hull approximately 
.00040 AV?. The NC hull is exceptionally clean, and the value 
of K is probably about .oooso for the average flying boat hull. 


LENGTH 
DIAMETER 
XO 2 4 6 8 
Figure 57. Effect of Fineness Ratio on K and min Equation for Drag of 
C-Class Airship 


Streamline bodies. The drag of streamline bodies at 0° 
pitch and yaw is not only very difficult to measure accurately but 
it is also greatly affected by surface irregularities, disturbance, 
or turbulence in the air flow, and by change in Reynolds Num- 
ber. The form having the least drag per unit of cross-section 
area does not necessarily have the least drag per unit of 
volume, and the form which shows up best at one Reynolds 
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Number may be comparatively poor at some other Reynolds 
Number. 

The C-Class airship is a form which has a very low drag under 
all conditions. This shape has been thoroughly tested at the 
Washington Navy Yard for the effect of inserting parallel middle 
body? and for the effect of varying the spacing between stations. 
The drag may be written in the form 


p= x (2) (¥8) “ 


where D is the drag in lbs., V the air speed in ft./sec., and d the 
diameter in ft. Both K and 2 vary with fineness ratio as shown 
on Figure 57. Allowing for the fact that the value of is very 
difficult to determine, the test data are quite consistent. 

Offsets for the C-Class are given in Table 3. 


Cellular radiators. Tests at the Washington Navy Yard*"4 
on cellular radiators gave values of Cp between 0.57 and 0.78, 
with an average value of about 0.70. This coefficient is for 
the core alone and it decreases slightly as the core thickness 
is increased. The average Cp corresponds to about 18.0 lbs./sq. 
ft. of core at 100 mi./hr. 

When a cellular radiator is covered by closed shutters, the 
drag coefficient is substantially that of a flat plate, i.e., Cp = 1.10 
approximately. Allowance must be made for the additional 
resistance of the internal flow over engine, etc., when calculating 
the drag of nose radiators. The drag of headers on free-air 
radiators should be considered as that of a flat plate of the same 
area. 


Lamblin radiators. Tests at Issy-les-Moulineaux'5 on a 
Lamblin radiator 10” in diameter by 23.6” in length, gave D = 
.0038 V ? in kg.-m.-sec. units which is equivalent to 


73 A. F. Zahm, R. H. Smith, and G. C. Hill, oe Drag of the C-Class A 
hee Varying Lengths of Cylindric Midships,’”” N. A eu Technical Boe et eer 
1922). 


4 R. H. Smith, ‘‘ Resistance and Cooling: Power of 
Tocheeeoes aes peak ing Power of Various Radiators,” N. A. C. A, 


tS “Bulletin Technique” No. 28, August, 1925. 
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D = .001675V? lbs. 


for V in mi./hr. at standard density. Tests at the Washington 
Navy Yard’® on a larger size Lamblin radiator gave 


D = .0023 V? lbs. 


Drag of Sperry Messenger airplane parts. The drag of var- 
ious parts of a full-scale Sperry Messenger airplane has recently 
been measured in the N. A.C. A. 20-ft. wind tunnel at Langley 
Field.t? The results were as follows: 


Ibs. Drag 
Part at 80 mi./hr. 
BSMOMIOSER TISCIARE 0 ord v stds ie es Bee es es Oe 15.95 
Engine (3-cylinder air-cooled)................... II.30 
(CXS GUISE AR, ARRten oF 6 Os DENIS eare ne RE rei meeres ee 2.00 
Windshield and upholstering................:0. 2.00 
Pan sumaces. 3.555. 62142 Ub Miers Oa apiceeln pape 8.00 
MATIN SOG CAl fas rorieeisca setele cles tec Oa 26.00 


The drag coefficient of the bare fuselage with faired-over 
cockpit is approximately K = .c0o42. The addition of the wind- 
shield and upholstering and the opening of the cockpit raises the 
value of K to approximately .0048. 

The drag of the 3-cylinder air-cooled engine is equivalent to 
D = .00176 V? for V in mi./hr. On this basis the drag of a simi- 
lar 9-cylinder air-cooled engine, without additional interference 
allowance, would be D = .0053 V?._ This is comparable with the 
average value D = .oo60 V? obtained from flight test data on 
the Wright ‘‘ Whirlwind.” 

The total area of the tail surfaces is about 30 sq. ft. The 
measured drag indicates a value of Cp = .o11 (K = .000042). 
This is equivalent to about 0.40 lb. per sq. ft. at 100 mi./hr. 

Air-cooled engines. The drag of an air-cooled engine de- 
pends greatly on the cowling. With average cowling the drag of 
the Wright “ Whirlwind ” is about .oo60V ? lbs., V being in mi./hr. 

16See N.A.C.A. Technical Report No. 261. 


17 Fred E. Weick, ‘‘Full-Scale Drag Tests on Various Parts of Sperry Messenger Air- 
plane,’ N. A. C. A. Technical Note No. 271 (1928). 
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Flight test data indicate values for this engine from approxi- 
mately .cosoV? to .o1oV?. The drag of the Pratt and Whitney 
“Wasp” has been calculated for normal cowling as D = .oo80V?. 
The drag of the Pratt and Whitney “Hornet” or the Wright 
“Cyclone,” on the same basis, would be D = .cogoV’. 

Data now available show a consistent indication that the drag 
of a 9-cylinder air-cooled engine is equivalent to an increase 
between 1 and 2 sq. ft. of parasite flat-plate area over the value 
for the corresponding water-cooled installation with a good free- 
air or tunnel type radiator. Most of this increase is unneces- 
sary and may be avoided by proper cowling. 

Wheels—landing gear. British tests'® gave the drag of 
wheels as follows: 


(a) Wheel without fairing, Cp = 0.70 
(D = 18.0 lbs./sq. ft. at 100 mi./hr.) 
(b) With the usual fairing from hub to rim Cp = 0.43 
(D = 11.0 lbs./sq. ft. at 100 mi./hr.) 
(c) With complete fairing from hub to edge of tread, 
Cp = 0.23 
(D = 6.0 lbs./sq. ft. at 100 mi./hr.) 


French tests (Bulletin Technique No. 28, August, 1925) on 
the complete landing gears for the Spad 33 and Spad 40, in- 
dicated a total drag of 55 and 61 lbs. respectively at 100 mi./hr. 
In both cases about 31 lbs. of this was due to the two wheels and 
the remainder to the four se and faired axle. The wheels 
were 800 x 150 mm. or about 32” x6”. The drag of one 32 
6” wheel is, therefore, about 15.5 Ibs. at 100 mi./hr. 

Tests on an 800 x 150 mm. Breguet shock-absorbing wheel 
fitted with side plates were reported in the same bulletin. This 
wheel had a drag of 16 lbs. at 100 mi./hr. These values agree 
very closely with flight test data. 

The drag in lbs. per wheel at 100 mi./hr. in standard air is 
approximately: 


Br. A.C.A. R.& M. No. 207. 
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Ibs. Drag ie 

Size of Wheel at 100 mi./hr. 
20x 3 7-5 
28x4 10.0 
30X 5 T25 
32x6 16.0 


Machine guns. Flight tests indicate that the drag of a 
machine gun, lined fore-and-aft, with a scarff ring mount, is 
about .co20V?. This value is increased about 50% by putting 
the gun broadside. 


Searchlights. A searchlight 21 cm. (8.27’’) maximum dia- 
meter by 42 cm. (16.55’’) in length, complete with wing-mounting 
brackets was tested at St. Cyr*® and the drag found to be 2.9V? 
in gr-m-sec. units which is equivalent to 


D = .00132-¥7 Ibs. 
for V in mi./hr. at standard density. 


Interference. The ‘‘interference” effect between two ob- 
jects varies greatly, sometimes decreasing the total resistance 
and sometimes increasing it. The reduction of drag when two 
wires are lined up fore-and-aft has been given. When two 
wires are side by side the total resistance, for spacings less than 
about 6 diameters, is greater than twice that of a single wire, 
but very little actual test data are available. It is usually 
assumed that no interference exists for lateral spacings of 10 
diameters or greater. This has been verified at the Washington 
Navy Yard, but no record was made of the tests. 

The ‘“‘upstream wake” of a cylinder, which is very aptly 
designated as an “‘air prow” by Dr. Zahm, extends to a distance 
of about 5 diameters and the lateral disturbance extends to a 
similar distance before the velocity is normal. The theoretical 
velocity distributions about a number of simple shapes are 


19 ‘Bulletin Technique” No. 28, August, 1925. 
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shown in N. A. C. A. Technical Report No. 253.7? From these 
one may secure an idea of the extent of the disturbance created 
by a moving object. 

The drag of two adjacent struts was partially investigated by 
Nayler and Jones.?! Two 1’’ x 3” streamline struts were tested 
at several spacings and the average drags found to be as follows: 


Spacing _ © 4.90 3.970 2.890 
Thickness 
Total Drag, lbs. 0.048 ©.050 0.052 0.053 


These data indicate that two struts spaced laterally 3 ‘‘dia- 
meters” apart on their center lines have about 10% more drag 
than with normal spacing. 

The most important form of interference is that which occurs 
when anything is attached to an airplane wing, particularly on the 
upper surface. This may be illustrated by the results of some 
tests on a twin engine airplane model at the Washington Navy 
Yard. The drag in the first test appeared high and the nacelles 
were removed and tested alone. At 40 mi./hr. the nacelle drag 
was found to be .078 Ib. in free air, and 0.120 lb. on the model. 
A. plasticine fairing with generous fillets between the nacelle 
and the wing reduced the interference drag from .042 lb. to .o02 lb. 
This model was a biplane with the nacelles mounted directly 
on the upper surface of the lower wing. A still greater effect 
probably would be found if the nacelles were located on the 
upper surface of the upper wing of a biplane or on the upper 
surface of a monoplane wing. When it is necessary to mount 
nacelles above a wing, wind tunnel tests with and without 
nacelles are advisable in order to determine whether or not the 
fairing is adequate, or to find the fairing which gives a minimum 
interference drag. 

Any projection from the upper surface of an airplane wing 
such as strut fittings, bolts, brackets, etc., will show interference 


20 A. F. Zahm, ‘‘Flow and Drag Formulas for Simple Quadtics” (1926). 


21 “The Determination of the Forces on Two Struts in Close Proximit , 
Br. A.C. A. R.&M. No. 204 (1915). ; Sete ricer cit 


Ch. 4] PARASITE DRAG DATA gi 


_ drags ranging from 50% to 200% of the free air drag of the part 
_ in question. Any form of fairing between wings and fuselage 
that interferes with the natural flow of air over the wing may be 
expected to show a pronounced interference effect. Horizontal 
tail surfaces fitted just above the fuselage may also be expected 
to show interference drag effects if the gap is not filled. 


CHAPTERS YV 
CONTROL SURFACE DESIGN 


Axes, forces, and moments. The three major axes of an air- 
plane are the longitudinal or X, lateral or Y, and normal or Z. 
Motion about the X, Y, and Z axes is designated as roll, pitch, 
and yaw, respectively. Control about these axes is secured by 
the ailerons, elevators, and rudder, respectively. By universal 
agreement in this country the positive direction along the axes is 
to the rear for X, to the left (as seen from the rear) for Y, and up 
for Z. The positive direction for angles or moments is from + Y 
toward + Z (or right wing down) in roll, from +Z to +X (or 
nose up) in pitch, and from +X to + Y (or nose swinging to 
right) in yaw. 

Stability. An airplane is statically stable if any displacement 
from a given attitude sets up forces and moments tending to 
restore the original attitude. It is dynamically stable if the result- 
ing motion is stable, that is, if any oscillations due to static 
stability are quickly damped. Static stability may be easily 
measured by wind tunnel tests on an airplane model; it is directly 
proportional to the slope of the moment curve. Dynamic 
stability, on the other hand, must be laboriously calculated from 
rather extensive wind tunnel tests, using assumptions which are 
questionable. It is not surprising, therefore, to find that practic- 
ally all aeronautical engineers ignore dynamic stability as a 
design factor. 

In pitch, a fair degree of static stability with normal design 
proportions is usually accompanied by dynamic stability. 
Stability in pitch is called “longitudinal stability.” Stability 
in roll and stability in yaw are not easily separated. They are 
always treated in combination as “lateral stability.” Static 

92 
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stability in both roll and yaw does not insure dynamic lateral 
stability and in fact lateral instability may easily result from too 
much “directional stability.”” In this connection, the use of 
moderate dihedral is probably desirable in all seaplanes and in 
most land planes. The prejudice of many pilots against the 
use of dihedral in pursuit type airplanes is not substantiated by 
full-scale rolling and turning data. 

The study of dynamic stability is of greater value to the 
aeronautical engineer than is generally realized. It is an ex- 
cellent method of analyzing the effects of changes in design 
proportions and will probably be used more extensively for 
this purpose in the future when the demands for more refinement 
in design methods make it necessary. 

The literature on stability calculations is extensive. The 
most important papers are: 

L. Bairstow, “Investigation into the Stability of an 
Airplane, Hte:” Br. A-C.A. R:&M. No. 77 (1913). 

L. Bairstow, “‘The Experimental Determination of Rotary 
Coefficients.’ Br. A.C.A. R.&M. No. 78. 

L. Bairstow, ‘‘An Examination into the Longitudinal 
Stability of a Monoplane.” Br. A.C.A. R.&M. 
No. 79. 

J. C. Hunsaker, ‘Experimental Analysis of Inherent 
Longitudinal Stability for a Typical Biplane.” 
N. A. C. A. Technical Report No. 1. 

Excellent treatment of stability is also given in the following 
books: 

L. Bairstow, ‘‘ Applied Aerodynamics.’’ Longmans, Green 
fCo. (1020). 

E. B. Wilson, ‘‘ Aeronautics.’”” John Wiley & Son. (1920). 

Control. One of the most important problems in airplane 
design is to determine the sizes and proportions of the control 
surfaces so as to obtain a satisfactory degree of control and 
stability under all conditions of flight. The remainder of this 
chapter is devoted to the general question of control surface 
design. 
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Horizontal tail area. First Method: The horizontal tail area 
required to give static longitudinal stability may be obtained by 
writing the equation for pitching moments about the center of 
gravity. If the resultant air force be assumed normal to the 
wing chord and equal to the lift (which is a close approximation 
except at very low lift coefficients), the slope of the wing moment 
curve against angle of attack is 


dC. Cp 
Be = gSuC [2- (Co ray a) (43) 


where S, is the total wing area, C the mean wing chord, a the 
fore-and-aft center of gravity location on the mean chord, C, 
the center of pressure, a, the absolute angle of attack measured 
from zero lift, and g the dynamic pressure. The slope of the 
tail moment curve is 


di da da (44) 


where @, is the angle of attack, Cz; the lift coefficient, and S, 
the total area, of the horizontal tail surfaces. J is the distance 
from the c.c. of the airplane to the center of pressure of the 
tail surfaces, and may be taken for convenience as the distance - 
from the c.c. to the hinge axis. 
The slope of the total resultant moment curve for the air- 
plane may be written 
aM _ Kqwc (45) 
where W is the gross weight, C the mean wing chord, and K a 
constant depending on ee Meee 
dM dM, 


Assuming that —- = ne + oe 
da 


(43), (44), and (45) gives 


~ and combining equations 


KqWC = qS.C% | % — (c, + a2) | —, oe 
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ae d . 
Dividing by (g5eC Wu, =) and arranging terms 
da; da 
ek I af i) a ( ay dCr 
da, da ; 
pC doy _ ( dC,\ - dC, 
_ letting iS Big et 24 Cota at) =F; 5 and ae 


= F,, then equation (46) becomes 


onmn|-k()+G-5)%] on 


F, and F, are the slopes of the lift curves for tail and wings 
respectively and depend on the effective aspect ratio and section. 
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Figure 58. Variation of Slope of Lift Curve with Effective Aspect Ratio. 
Stability Factors Fr: and Fy 


_ Values of F, and F, may be read from Figure 58. For example, 


if o = .072 at aspect ratio 6, then ot = .065 at aspect 
4 


ratio 4, or .o76 at aspect ratio 8. The variation with aspect 
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Figure 59. Variation of the (Downwash) Stability Factor F2 with Effective 
Aspect Ratio and Tail Length 
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ratio is along a curve similar to those shown and if rs is known 
a 


for any aspect ratio, the value is found at any other aspect ratio 
by following along or in between the curves given. The average 


value of a is about .072 at aspect ratio 6. Note that this 
% 


is the effective aspect ratio n = re “(See Chapter II for 
values of k.) | U7 

F, is a wing section stability factor depending on the camber 
and trailing edge curvature. Values of ot and F, for many of 
the wing sections in common use are given in Table 8. 


F,, is a downwash factor depending on the effective aspect 
ratio of the wings and the tail length /, as shown on Figure so. 
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The static stability constant K has been obtained from wind 
tunnel tests on a great number of satisfactory airplanes and 
varies between —G005 and —.oo10. The former value is 
about the minimum for static stability at all speeds, while 
values greater than — .ooto give excessive stability or ‘“‘stiffness.” 

Equation (47) has given uniformly excellent results for design 
purposes. It is of interest to note that the left-hand side is the 
familiar “‘th.”” This is not a constant as many designers have 
assumed; it varies with design proportions. The author has 
found values varying from 0.25 to 1.0 and the extreme limits 
were certainly not attained. 

The effect of each variable on the horizontal area required 
for constant static stability has been calculated. Figure 60 
shows the effect of wing aspect ratio, Figure 61 the effect of 
tail aspect ratio, Figure 62 the effect of tail length, and Figure 
63 the effect of fore-and-aft c.c. location. 
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Figure 60. Effect of Wing Aspect Ratio on Horizontal Tail Area Required 
’ for Constant Static Stability 
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Figure 61. Effect of Tail Aspect Ratio on Tail Area Required for Constant 
Static Stability 


For constant static stability the area of the horizontal tail 
surfaces continues to decrease as the c.c. is moved forward, but 
in order to provide control, the decrease cannot be utilized. 
Experience indicates that control and stability combine to give 
a minimum area for a c.c. location at about 30% of the mean 
chord. 


Second Method: The horizontal tail area may also be calcu- 
lated by writing moments about the c.c. at zero lift, when the 
downwash is zero. The wing moment is 
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where Cy, is the moment about the leading edge (or any point) 
at zero lift, S,, the total wing area, and C the mean chord, as 
before. The tail moment is 


OG: 
M; = Cuq Sil = ao Go g Sil (49) 


where «a, is the “‘effective longitudinal dihedral,” measured be- 
tween the no-lift lines of the wings and tail, and the other terms 
as in equation (44). The resultant moment on the airplane at 
zero lift should be given by 


MM, = Eg W-C (50) 
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Figure 62. Effect of Tail Length on Tail Area Required for Constant Static 
Stability. 2 = Effective Aspect Ratio 
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Figure 63. Effect of c.c. Location on Tail Area Required for Constant Static 
Stability 


from equation (45). Equating these moments gives 


dCi 
dai 


ee eras 
#-2- Z| x. (2 Cue | (51) 


Analysis of wind tunnel test data shows that 


K, = Ka,’ 


Ay 


OST CingSeC = Rear 


from which 


where a,’ is the absolute angle of attack for balance of moments 
and K a constant varying from about .coo5 to .ooro according 
to the stability desired. K has the same value as in equation 
(45), but it is now positive, because a,’ is considered positive. 
The analysis also shows a definite relation between the “stabilizer — 
setting” «, and the absolute angle of attack for balance. The 
average value is 


% = (3.0° + 0.25 a’) 


Equation (51) may therefore be written 
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S; hie I : W 
Sw» CF, (3.0 +0.25aé) | Koa =) ti Cus | (52) 


The moment coefficient at zero lift, Cy, is normally negative, 
and may be found in Table 8 for many of the wing sections in 
common use. 

As an example of the calculation for horizontal tail area, 
consider an airplane having the following characteristics: 


Sw 


Effective wing aspect ratio = 4.55 
Effective tail aspect ratio = 3.35 


oe = 11.20 lbs./sq. ft. Clark Y section 


Tail length ? = 2.85 


Balance at « = 0° ee SAS 
c. G. location at 32% of mean chord. 


TABLE 8. WING SECTION STABILITY DATA 


Slope of 


2 Bel Moment 
eee ( C EN Coefficient 
- Wing Section Ratio 6 eS aa Pe at 
GG. F; Zero Lift 
da Cu 

PREM HTS... ce. 077 25 —.052 
MOAI Oso ste ate Succes 080 24 — .060 
CE ee es O71 23) — .080 
Gottingen 387...... .072 LM —.095 
Gottingen 398...... FO72 23 —.079 
Gottingen 420...... .075 .21 — .090 
pe eineen ABOn eee .072 23 —.078 
US eee o71 23 — .086 
OAK 35A. nes 073 20 —.120 
Mlaivey IN-Ov.. nie ees: 072 24 — .060 
Bee IN=TO siete a 26s « 080 23 —.075 
Navy N-22 074 22 — .083 
Pew Wie6) 05) 072 27 +.010 
: — .005 
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From these data 


F, = 0.0613 F, = 0.520 

F, = 0.23 F, = 0.0662 
from which by equation (47) J 
a 5 “ = Seer [.0006 X 11.20 + (0.32 — 0.23) X 0.662] 

= 0.397 
and »y equation (52) 
2 : ‘ ere rere art [.0006 X 5.4 X 11.20 + .080] 
= 0.434 


a coefficient of about 0.42 is probably a good compromise. 


Elevator and stabilizer areas. The best distribution of area 
between the stabilizer and elevators obviously depends on 
several design factors, but in general the elevator area should lie 
between 35% and 50% of the total. The recommended range 
is between 40% and 45%. In a seaplane, values below 40% 
usually lead to deficient control in take-off. 


Tail-surface sections. If tail surfaces are to give maximum 
stabilizing and control effects, the section must be as thin as 
practicable. The chief disadvantages of thick surfaces are: (1) 
low value of ee. (2) low full-scale maximum C,, (or Cy); (3) un-— 
steady air flow leading to irregularities in moment curves and 
control effect; and (4) sensitiveness to slight irregularities in 
section due to method of construction or to handling in service. 

It is particularly desirable to avoid rapid changes in section. 
The sections should preferably be geometrically similar, but 
slight deviations for structural reasons are allowable. Over- 
hanging balances must have reasonably thin sections. 

Wind tunnel tests indicate that aerodynamically there is no 
one section that is notably superior to the others. In the 
matter of minimum drag, the Gottingen sections Nos. 443, 444, 
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and 445 show up very well, but these are not entirely satisfactory 
in regard to beam depths. A section similar to that shown in 
Figure 64 is apparently the best compromise between aero- 
dynamic and structural requirements. This section is shown 


THICKNESS /CHORD 


“i 
° 


70 ~=«.80 
DISTANCE FROM LEADING EDGE 


Figure 64. Airfoil Section for Tail Surfaces 


with sharp leading and trailing edges, but any reasonable 
radius may be used on either, or both, leading and trailing edges. 
The sharp edge is not recommended. The maximum thickness 
should not be greater than about 8% of the chord. 


Plan-form of horizontal tail surfaces. Figure 65 gives the 
most frequently used plan-forms for horizontal tail surfaces. 
Form A is approximately rectangular with a modified elliptical 
tip. Form B has a circular or elliptical tip with raked, but 
straight, leading, and trailing edges. Form C is elliptical, and 
Form D is approximately parabolical or triangular in shape. 


Vertical tail surface required. The only certain method of 
designing vertical tail surfaces is by wind tunnel tests, but a 
close approximation may be secured by a simple formula. Let 
the full-scale yawing moment at, say, 10° yaw, be N,, lbs-ft., 
the total vertical tail area Sy sq. ft., and the tail length (from 
c.c. to rudder hinge axis) / ft. Then it may be expected that 


Na _ (2M) 
re aa ee RngSyl (53) 
If &,, as obtained from wind tunnel tests be plotted against J, 
the points are found to fall into three fairly definite groups 


representing landplanes, seaplanes, and flying boats. 


penny. in Bas. ON ot cal 
fn % y¢ mm im Sf \ ww? 


[ 
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Tigure 65. ‘Tail Surface Plan-Forms 
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Figure 66. Effect of Tail Length on Yawing Moment Slope Coefficient 
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« 
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Figure 66 is a plot of k, against 7 for about fifty airplanes 
for which wind tunnel data are available. k, appears to vary 
directly as J and has average values of .oor 5 for landplanes, 
.oort for seaplanes, and .00075 for flying boats. From which 
k, = K,l and for 


Landplanes (=) = .0015 gSyi? (54) 
dN 

Seaplanes (F )= .OOII gSyl? (55) 

Flying Boats (a) = .00075 gSyl? (56) 


It will be noted that the points for part of the flying boats 
are on line B in Figure 66. These have a comparatively short 
length of hull forward of the wings, as compared with the older 


types. i is the effective slope of the yawing moment curve 


for an appreciable angle of yaw, and is used in preference to the 
actual slope at o® yaw on account of the difficulty in measuring 
the latter. 

It is logical to assume that the desirable slope of the moment 
curve should be proportional to the type of airplane, its gross 
weight, and maximum ee That is _ 


ai = Kyqwb (57) 
Analysis of wind tunnel test data on about fifty airplanes, in 
the light of flight test reports on the same airplanes, shows that. 


_ the best values of Ky are between 4 X 10 —5 and 8X 107%. 


The desirable limits for various types of airplanes appear to be 


Pursuit, Ky between 4 X 10-5 and 6 xX 10 5 

Observation, Ky between 5 X 10-5 and 7 X IO—5 

Bombing, training, commercial, etc., Ky between 
60:30 **-and'to X 10% 


(Fy). from equation (53) and ( a from equation (57) may 


now be equated, giving. 
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K,qSvV?’ = KyqWb 


from which 


_ Ky Wo 
Sy == K, P : ‘ (58) 


The values of Ky have just been given, Avhile K, is .oo15 
for landplanes, .oo11 for seaplanes, and,doo75 for flying boats. 
For equal directional stability, a seaplane requires about 35% 
more, and a flying boat 100% more, vertical tail surface than a 
landplane of the same gross weight, span, and tail length. 

As an example of calculation for vertical tail area, consider a 
pursuit type landplane for which W = 2,800 lbs., 6 = 32 ft., 
andl = 15.0ft. K,= .ocors and Ky should be about 5 X 10 —5. 
Bes 107° \ 2,800 SO:32 


se .OO15 225 


= 13.3 sq. ft. 
As a seaplane, this design would probably weigh about 10% 
more, or W = 3,080 lbs. For equal directional stability, since 
K, = .oor1 for a seaplane: 
=—5 
ay y¢ 3080 X 32 
.OOII 225 


= 19.9 sq. ft. 


If the same vertical surface is to be used on both landplane 
and seaplane, the area should be about 15 or 16 sq. ft. 


Rudder and fin area. The correct division of the vertical 
area between the rudder and fim depends upon the type of 
airplane; landplane, seaplane, or flying boat. It also depends 
upon the length and shape of the fuselage, hull, and floats, and 
upon the wing arrangement. The only satisfactory method 
is to make wind tunnel tests, but analysis of available wind 
tunnel test data shows that a majority of airplane designs will 
be satisfactory with the following proportions: 


Rudder Area 
Type Total Vertical Area 
Bandplanes® cc cede suai taee -55 to .60 
SEADIANES cc acc Mae Cee 50, tO) .35 


Flying Boats 0.0. chur ees -45 to .50 
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Aileron area required. The effectiveness of an aileron, as 
measured by the rolling moment due to a given setting, varies 
with the aileron span, and chord. It may be shown? from test 
data that the variation of efficiency with chord is given by 


Capers ‘) 
Sates 2.00 e (59) 


where ¢ is the aileron chord and ¢ the total wing chord, including 


«TE se i Ci iD 
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Figure 67. Aileron Proportions. Relative Span and Chord 


the aileron chord. Ina similar manner the variation of efficiency 
with span is 


= 1.20 —0.60(F) (60) 


where / is the aileron span and 0 the wing span. . 
Assuming that satisfactory control is given by an aileron 


of the proportions Lae 0.25 and iB = 0.40, it is possible to calculate 
C 


I W.S. Diehl, ‘‘Notes on the Design of Ailerons,” N. A. C. A.iTechnical Note No. 144 
(1923). 
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Figure 68. Aileron Proportions. Relative Area and Chord 


the proportions of all other ailerons having the same effectiveness 


in producing rolling moments. The assumed aileron has Sie 


oe 
0.10, where S, is the aileron area and S, the wing area. The 


2 = 3.00 and “2 =0.96, and the “effective 


No No 


area”’ is the product of these efficiencies by the actual area, or 


relative efficiencies are 


Se = 0.10 X 1.00 X 0.96 = .096 


The proportions of all ailerons having this effective area are 
given by 


|r. 50 — 2.00 ( alli it 20 — 0.60(5) | = .096.c 
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If the satisfactory aileron proportions be assumed as ee 0.20 and 
C 


l s : : A 
= = 0.40, the proportions of equivalent ailerons are given by 


b 
if 50 — 2.00 ( (2) | [. ches) og 


Figure 67 is a plot of . against “ and Figure 68 is a plot of 2 


w 


against : as calculated from the two equations. These are not 


the extreme limits for satisfactory ailerons, but the best results 

have been obtained from approximately the proportions given. 
The aileron when displaced from neutral affects the lift dis- 

tribution over the whole wing, but chiefly over the part of which 

it forms the trailing edge. If the area thus affected is S., with 

the center of this area at distance d from the fore-and-aft axis, 

then a convenient measure of control effectiveness is given by 

a S. 


K,= ba Sy ; (61) 
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Figure 69. Aileron Coefficient Ka 
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Figure 69 contains a plot of K, against f for the two effective 


areas previously considered. The points on this figure represent 
recent designs. Airplanes which were reported to have ex- 
cellent lateral control are indicated by circles with check marks 
(,/). A plain circle indicates no particular comment, while a 
square indicates unfavorable comment. An aileron design 
based on K, between the limits shown on Figure 69 will cer- 
tainly be satisfactory, in so far as area is concerned. 

The aileron chord should lie between the extreme limits of | 


a | 


= 0.16 and ‘ = 0.28. The most desirable chord, aerodynam- 
F ; t t 
ically, is between 7 = 0.20 and a 0.25. 


Aileron types. Figure 70 illustrates eight types of wing tips 
and ailerons. Types A and C are simple ailerons on elliptical 
tips. These are very satisfactory, and should be used when- 
ever practicable on all high-speed airplanes. Type B is a 
simple aileron on a modified elliptical tip. This aileron is 
slightly improved if the outer tip is cut away as shown by the 
dotted line. Type D is an inset aileron on a modified elliptical 
tip. This type does not extend to the wing tip. It is desirable 
on wings of great span and unusual flexibility. 

Types E, F, G, and H are undesirable for various reasons. 
All of these tend to be “heavy” on control owing to the peak 
in loading at the tip. Types G and H should never be used 
on very high-speed designs owing to probability of flutter. 

Care must be taken in the design of hinges for long, narrow 
ailerons in order to avoid binding when the wings deflect under 
load. 


_  Aileron moments. The best data available at this time on 
rolling moments and aileron hinge moments are given in R. & M. 
Nos. 550 and 615. These data were obtained at a test speed of 


2**An Investigation of the Aerodynamic Properties of Wing Ailerons,” Part I b 
H. B. Irving, E. Ower, and G. A. Hankins (1918), and Part II by Irving and Ower. (1910) 
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DESIRABLE FORMS 


Figure 70. Aileron Types 


40 ft./sec. on a biplane, gap = chord, no stagger and both 
wings 6” x 36”, R. A. F.—15. The wing tips were square 
(with a small radius at each corner, however) and ailerons were 
fitted on both wings. From these tests it appears that the 
hinge moment does not change very much with angle of attack. 
Figure 71 is a plot of the hinge moment coefficient Ky against 


the ratio - for various aileron angles at an angle of attack «,,=8°. 


The hinge moment for any condition is 
H = Ky-£--CS, Vv? (62) 


where C is the wing chord, S, the total aileron area, and V the 
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Figure 71. Aileron Hinge Moment Coefficients 


air speed. Ky must be taken as the average for the “up” and 
“down” aileron angles. This method of plotting allows dif- 
ferential movements to be studied. 

In R. & M. No. 515 the coefficient for rolling moment per 
unit aileron area is defined by the equation | 


L= Ky, £bS.V" (63) 
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where 6 is the span and S, the total aileron area. This form 
is about as satisfactory as any that can be written and it has 
the advantage of showing that the decrease in efficiency is 
almost exactly equal to the decrease in moment arm from the 
center line of the airplane to the center of the area affected by 
the aileron. That is, if instead of the span 6, the moment arm d 
were used, K; would show no appreciable variation with aileron 
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Figure 72. ' Aileron Rolling Moment Coefficients a = 8° 


span. This may be done by using instead of 0, the factor d = 


2 — a} where 0, is the total span of the ailerons. 
2 


Values of K, for various aileron angles are plotted against : 


in Figures 72 and 73 for angles of attack of 8° and 16° respec- 
tively. These values of Kz are based on ailerons of 6” span ona 
wing of 36” span. Consequently, the actual effective arm was 
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Figure 73. Aileron Rolling Moment Coefficients a = 16° 


$b instead of 6 and for exact results with any aileron span, the 
rolling moments should be obtained from 


L = Kr £dS.V? (64) 


2b — b, 
2 
tween d and the variation of the efficiency with span, equation 


(60). 


Balanced controls. Elevators, rudders, and ailerons are 
often “balanced” in order to reduce control forces. For small 
or moderate-sized airplanes the chief reason for balancing is to 
improve maneuverability, but for large airplanes it is necessary 
in order that the control forces do not exceed the pilot’s strength. 
The most desirable type of balance varies with the use and size 
of the control surface. For example, what is probably the most 


where d = ( ) : = (1.26 — 0.6b,). Note the relation be- 
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satisfactory type of balance for a rudder is unsuited for an 
aileron on a high-speed airplane. 


Figure 74. Balanced Controls 


In general, no type of aileron balance should be considered 
satisfactory unless it allows the use of an efficient wing tip, and 
tends to bring the aileron c.c. near to the hinge axis. Likewise, 
no type of rudder balance can be considered satisfactory if it 
reduces the maximum rudder control, or in any way reduces the 
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rudder control. Elevator balance requirements are similar to 
those for rudders. 


Types of balanced controls. Figure 74 shows six types of 
balances. Types A and B were formerly used on all controls, 
but they are now used only on rudders and-elevators. There is 
little choice between them, but type B is the more common at 
present. Neither A nor B should be used on ailerons owing to 
the poor wing tip and to the high peak loadings on the balanced 
portion. Type C is liable to overbalance and flutter. Type D 
has no advantages and is rarely used. Types E and F have 
been used to a limited extent, but have no outstanding 
advantage. 

Figure 75 shows three widely used forms of balances. Type I 
is usually called the ‘Handley Page balance.” It is satisfactory 
for control surfaces that need not be moved through large angles. 
With a normal airfoil section, the leading edge of the control 
emerges from the wake of the fixed surface at a control angle of 
about 15°. Above this critical angle no appreciable increase in 
control is obtained. For best results, the hinge axis should be 
located at approximately 25% of control chord from its leading 
edge. Type II is the familiar “‘paddle” balance. Full balance 
cannot easily be secured with this type but it has been fairly 
satisfactory on all surfaces. The auxiliary surface should be of 
large aspect ratio, and as far forward as practicable. The area 
should be generous for good balance. A symmetrical section is 
generally used, and on ailerons it is set at an angle of +3° to + 5° 
to the wing chord. On rudders and elevators, the auxiliary 
surface would be set between 0° and 1° to the main surface axis. 
The amount of balance is readily calculated. 

Type II is the Frise balance. This is the most satisfactory 
form of balance now available for ailerons. When the aileron 
is up, the leading edge projects below the wing and thus adds a 
drag which tends to counteract the yawing moment due to the 
drag of the down aileron. The hinge axis is preferably located 
below the aileron center-line and at about 22% to 25% of the 


Ch. 5] CONTROL SURFACE DESIGN 117 


aileron chord. The hinge axis should not be further aft than 

25% of the aileron chord. ; 
The effective area of a balanced control surface of the 

Handley Page or the Frise type is measured between the hinge 
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Figure 75. Balanced Controls 


axis and the trailing edge. That part of the surface forward of 
the hinge axis does not contribute to the control effect. 


Calculation for simple balances. For a long time naval arch- 
itects have balanced rudders by assuming: (1) uniform loading, 
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(2) that the center of pressure on a flat plate is at 0.20C, (3) that 
the center of pressure on a movable surface of chord C, which 
trails a fixed surface is at 0.331. With these assumptions the 
average center of pressure may be calculated and the axis placed 


Figure 76. Ilustrating Calculation for Rudder Balance 


forward of this point at a distance sufficient to avoid overbalance 
as determined by accummulated experience. This method is 
applicable to air as well as water, and it gives excellent results. 
As applied to a rudder balance, for example, consider Figure 
76. The rudder is divided into a number of strips A, B ,C, D, etc., 
spaced as necessary in order to secure accuracy. The center of 
pressure of each strip is located as indicated by the circles, 
assuming Cy at 0.20 C for all of the overhanging strips and at 
0.33 C; for all of the trailing strips. The area of each strip is 
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then multiplied by the distance of its C, from the axis, consider- 
ing distances forward as negative. The algebraic sum of these 
moments divided by the total area is the average distance of the 
C, from the axis. This distance should be between 4 in. and 
6 in. for satisfactory balance, the distance decreasing slightly 
with the size of the airplane. Four and one-half inches to five 
inches is probably the best location for any airplane having a 
gross weight between 2,500 and 10,000 lbs. If the first calcula- 
tion shows the center of pressure outside of the desired range, 
balance is added or subtracted as required and additional 
calculations made until the desired location is obtained. 


CHAPTER VI 


ENGINE AND PROPELLER CONSIDERATIONS 


B.Hp. at given r.p.m. : 
—__ et — be 
B.Hp. at 1,800 r.p.m. 
plotted against r.p.m. for various modern airplane engines, it is 
found that a single curve may be drawn, from which no point 
differs appreciably more than the variation between individual 
engines on test. Such a curve, based on the average power 
curves for nine engine types, is given in Figure 77. This curve 
may be expanded as in Figure 78 to take care of the maximum 
r.p.m. actually used. Figure 78 is of great value in performance 
estimation when specific engine data are not obtainable. 

_ It should be noted that most of the engines designed prior to 
about 1920 delivered their maximum power at what is now a 
normal operating r._p.m. For example, the Liberty ‘“‘peaked”’ at 
less than 2,o00or.p.m. Curve three on Figure 78 should be used 
in this case, and the value of NV taken as about 1,750r.p.m. instead 
of the 2,400 r.p.m. there given. 


Variation of B.Hp. with altitude. The variation of B.Hp. 
with altitude may be obtained either from altitude chamber tests 
or free flight tests. Using altitude chamber test data it can be 
shown’ that at constant r.p.m. 

nee = le) aH Temperature constant 
BHp.. \P, 
B.Hp. (F 


—0:5 
and Bip. = 7) Pressure constant 


© rap - (ey, Me) (65) 


tW.S. Diehl, ‘‘Engine Performance and The Det ination f Absol iling,’’ 
N. A. C. A. Technical Report No. 171. paren a emer 


General power curves. If the ratio 


120 
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In the standard atmosphere, pressure and temperature are con- 
nected by the relation 
‘ay a a 
(7.) = (=) ve 


and pressure and density by the relation 


(5) =(E) (67) 


Therefore, in the standard atmosphere 


B.Hp. ou < (ace (68) 
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Figure 77. General Brake Horsepower Curve 
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Figure 78. General Brake Horsepower Curves 


As a result of the N. A.C. A. flight tests on a DH-g airplane 
fitted with a Bendemann dynamometer hub on the Liberty 
engine?, Gove and Green recommend the relation 


Bip. 2s ef 
Bp. = 1.088 Fr .088 (69) 


which is given in the “Handbook of Instructions for Airplane 
Designers.’’3 


2W. D. Gove and M. W. Green, ‘‘The Direct Measurement of Engine Power on an 


fon in Flight with a Hub Dynamometer,’’ N. A.C. A. Technical Report No. 252 
1926). 


3 Engineering Division, Air Service, McCook Field (4th edition, 1925). 
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A replot of the N. A. C. A. data on logarithmic paper shows 
that it would probably be better to write 


tit @” 


(70) 


0 


and in fact (1.088 : — .088) is almost exactly equivalent to 


(®) ; % In view of the fact that these data do not extend to 
altitudes above 13,000 ft., while pressures and temperatures were 
not recorded in the report, there is some question as to what was 
actually measured, or rather to what altitudes or densities the 
measured data correspond. 

The actual difference between equations (68) and (69) is not 
great, as the following comparison will show: 


Altitude 2 (=) ae a 
ft Ro Po Po Ah 

1) I .000 I .000 I.000 ° 
5,000 .862 .836 .824 — 400 
10,000 .738 .695 .673 — 800 
15,000 .629 BS 7S -547 —I,100 
20,000 2533 .471 -442 —1,600 
25,000 .448 2373 1255 —1,300 


Ah is the difference in altitude represented by the exponent. 
A difference of about 8% in altitude or 6% in density is indicated. 

Equations (65) and (68) are believed to give the variation of 
B.Hp. with altitude (at constant r.p.m.) as accurately as required 
for any engineering purpose. The variation of T.Hp. with 
altitude is given later. 


Specific fuel consumption. The specific fuel consumption 
(lbs./B.Hp./hr.) varies with the condition of an engine and the 
carburetor adjustment, but the effects of throttling, compression 
ratio, and mixture control are fairly definite. 

The effects of throttling with and without mixture control 
are given on Figure 79, which is based on the average of a great 
number of ground tests. The variation of specific fuel con- 


124 ENGINEERING AERODYNAMICS [Ch. 6 


sumption with compression ratio is given on Figure 80 for the 
two conditions of normal mixture and lean mixture. The latter 
is about the minimum that can be obtained in flight operation, 
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Figure 79. Specific Fuel Consumption at Part Throttle 


although it is possible to secure slightly lower values on the test 
stand. 

The variation in specific fuel consumption with altitude is 
given on Figure 81. This curve is based on the Bureau of 
Standards altitude chamber tests. 4 


Correcting test-stand data to standard conditions. Since the 
brake horsepower of an engine varies with pressure and tem- 
perature, it is necessary, in general, to correct the observed data ~ 
for deviations from the standard P= 29.92” Hg., t=15°C. While 


4N.A.C. A. Technical Reports Nos. 102 and 103. 
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altitude chamber test data indicate that the B.H Pie tS T= Ss: 
the deviations from standard pressure during a test are normally 
quite small and by general agreement the correction is made 


50 


.46 
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SPECIFIC FUEL CONSUMPTION 


38 
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Figure 80. Variation of Specific Fuel Consumption with Compression Ratio 
and Mixture Control 


on the basis of B.Hp.<«PT—°*5. If the tests are made with a 
dynamometer, this correction may be applied directly, but if 
the tests are made with a propeller or a test “club,” allowance 
must be made for the effect of air density on the power absorbed 
by the club. 

At constant r.p.m., the power absorbed by the club varies 
directly as the air density, or 
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BHp.= Bilpm 5 “ (71) 


where B.H p.m is the B.Hp. absorbed by the club at P = Po, and 
T=T,.. That is, B.Hp.m is the value read from the power 


"@ 8000 10000 15000 20000 25000 
Figure 81. Variation of Specific Fuel Consumption with Altitude 


curve for the club. The power delivered by the engine is the 
same as that absorbed by the club but the engine power is 


pe 2 Te 
B.Hp.a = BHp.s- V oa (72) 


equating B.Hp., and B.Hp., and solving for B.H.,. 


Piln P fe 


B.Hp., = “Ea = BHAp.m-: jj (73) 


This is on the assumption that the air temperature is the 
same at the carburetor and propeller. If these temperatures 
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differ and have the values T, and T, respectively, equation (73) 
becomes 


Te, 1, 
B.Hp.s = B.Hp.m ye (74) 


Propeller diameter. The best diameter for a propeller is not 
sharply defined since it varies with design constants, such as 
blade area, blade sections, plan-form, etc., which are also variable. 
However, the propellers commonly in use do not differ greatly 
in the minor proportions so that it is possible to calculate the 
diameter by formula and secure results satisfactory for general 
design purposes and performance calculations. 


If the power coefficient C, be plotted against ne for a family 


of geometrically similar propellers, it is found that the points for 
maximum efficiency lie on straight line passing through the 
origin. That is 


Sa ee 
a oe nD 


from which 


, Ban VERE Oe 
D= KM ap Kye Comyv 7) 


For V in mi./hr. the values of K in common use for ogee 
purposes are approximately as follows: 


WoopDEN PROPELLERS 


BIG SOE SAGES ele. 5 civic |e scaler a e+ 2 3 4 
K Speed (optimum slip)............. 285 255 240 
f Glimbi(Zero sup) hee sss ee ee 320 290 270 
a COMDEOMISE Se uta ede cle t cat a9 wel 305 275 260 


ING wOfeIAdeSte amare hicle seitiere et ee ot 2 3 


Speed (optimum slip)............. 280 260 
f CUMAD AGEING SUD) cc. s aiers oe 00s oo 8 315 300 
UC ORG DEGINISE) <cs title sles c'e-v ewes ee 300 290 
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The solution when, ® = 1.00, is conveniently made if 


0 


equation (75) is written in the form 


This is easily handled on the A, B, C, and D scales of any slide 
rule; divide K? (average value is 90,000) on D scale by r.p.m. on 


A : KT he 
C scale, set runner to index which reads op on A scale, 


move slide until V (mi./hr.) is under hair line, (dividing by V), 
move runner to B.Hp. on B scale (multiplying by B.Hp.). The 
hair line now indicates D4 on the A scale and D? on the D scale, 
from which the diameter is determined. This operation may be 
performed very quickly and is often necessary in making perform- 
ance estimates. An average value of K? = go,ooo is sufficiently 
accurate for all general purposes, such as performance calcula- 
tions, preliminary design studies, etc. 

A very convenient method® of igre diameter fae 


maximum efficiency is to calculate F = — y/o eV? and use 
Figure 82 to find 4 for maximum Pr The diameter is 


Us 
obviously D = a If V is in mi./hr., 2 in r.p.m., and P in 
nD 


oe 


B.Hp. 


‘ 


B.Hp., then F = 0.325 7 


3 


General efficiency curve. The curve of efficiency against ae 
n 
is similar for all conventional propellers,® and if the ratios of 


5 For additional data see W. S. Diehl, “The Application of Propeller Test Data to 
Design and Performance Calculations,” N. A. C. A. Technical Report No. 186 (1924). 

6 W.S. Diehl, “The General Efficiency Curve for Air Propellers,” N. A. C. A. Techni- 
cal Report No. 168 (1923). 
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Figure 82. Chart for Determination of Propeller Pitch, Diameter, and 
Maximum Efficiency 


= be plotted against ee a (5). where (5). corresponds to 
the maximum efficiency fm, it is found that the points lie on 
a single curve. This curve, which may be designated as a 
“general efficiency curve,” is given on Figure 83. From it the 


V 
efficiency curve for any propeller may be drawn when (|) 
and }m are known. i 
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—— 


Figure 83. General Efficiency Curve for Propellers 


Maximum efficiency. The maximum efficiency for any 


propeller of a geometrically similar series depends on slip 


V 
nD’ 
stream obstruction, and tip speed. Figure 84 gives the normal 
maximum efficiencies for conventional two-bladed wooden and 
metal propellers with unrestricted flow. 

When the slip stream is obstructed, as by a fuselage or nacelle, 
there is a reduction in efficiency depending on the size and shape 
of the interfering body. The approximate value of this reduction 
is given by Figure 85.7 The average reduction is about 2% or 
3%. The body diameter should be taken at a point about one 
propeller diameter aft of the propeller. 

The effect of high tip speed is very noticeable on a propeller _ 
having thick blade sections. No appreciable effect is found below 
about 0.8a, where a is the velocity of sound. Between o.8a and 
1.0@ the efficiency decreases rapidly with a total change which 


__ 7 Derived from data in British A.R.C. R.& M. No. 830, ‘Experiments on Airscrews 
with Tractor and Pusher Bodies—Part II,’’ by Fage, Lock, Bateman, and Williams. 
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Figure 84. Variation of Maximum Propeller Efficiency with aa 
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Figure 85. Effect of Body Interference on Propeller Efficiency 


70 


EFFICIENCY REDUCTION FACTOR 
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may be as much as 30% according to British tests. There is 
considerable evidence now available to show that the efficiency 
of a metal propeller with thin blade sections is not seriously 
affected for tip velocities below o.ga. 


Relation between s and the value of ~. for maximum effi- 
ciency. Definite relations exist between the nominal 4 ratio 


and the value of - for maximum efficiency. Designating this 


value of a by the subscript o, Durand’s data on wooden pro- 


pellers, show that for general service 


V 2 
(=), ~ 0-945 (5) ~ 46 (5) ™ 
(+) is sometimes called the “design (ae This equation 
nD/ o nD 


may be written in the form 


gyre Mla 


we nails 8 
(78) 
0.945 — -053 (55), 


Best climb is obtained when 


(=), = @) 


and maximum speed when 


(5), 0. (f) - o6() 


Adjustable pitch propellers. A large number of metal pro- 
pellers now in use have detachable blades which can be adjusted ~ 
on the ground to any desired blade angle, or pitch. The usual 
method of specifying the pitch of these propellers is to give the 
blade angle at a definite station. Figure 86 gives the pitch in 
feet corresponding to blade angles at the 42” radius. 

The maximum efficiencies obtained from metal propellers 
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v2] 8 9. wo 


4 5 & 
Figure 86. Relation Between Blade Setting and Pitch of Adjustable Blade 
| Propellers 

with adjustable blades are, for all practical purposes, identical 
with those given on Figure 83. Some slight differences are 
- certain to occur when a blade is used at a pitch angle widely 
different from its designed value, but data now available indi- 
cate that the characteristics are still determined by the effective 
| i, and the value of ae for maximum efficiency. 

Static thrust. The static thrust of a propeller is given 


closely by the empirical equation 
= oe 4 B.Ap.. 
T. = 6,000 Ez 9.5 ie ep D Ibs. (79) 


Where # is the nominal pitch and D the diameter. 
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This equation is based on data® for two-bladed, wooden 
propellers. No similar test data for metal propellers %available 
at this time, but the indications are that little difference exists. 

Negative thrust. The ratio of Us for zero thrust to ali for 

nD nD 
maximum efficiency, varies with blade area and P/D as shown 
on Figure 87. At values of a higher than that giving zero 


thrust the propeller acts as a brake and the thrust is negative. 
2.20 


1.80 


140 


: V 
Figure 87. me for Zero Thrust and Maximum Efficiency 


The value of this negative thrust may be expressed as a drag 
coefficient Cp defined by the relation 
b 


Negative thrust = Cp = 1B (2) J? (80) 


where V is the forward velocity, 6 the maximum blade width, 
and D the diameter. Values of Cp for 6 of Durand and Lesley’s 


propellers? are plotted against saa in Figure 88. It is of great 


8 From W. F. Durand and E. P. Lesley, “E i i ihe 
WAC A ye bucl eee ee so Crees). xperimental Research on Air Propellers, II,” 
9 Propellers Nos. 3, 4, 7, 8, 11, and 12. Data from W. F. Durand and E. P. Lesley, 


Gane Research on Air Propellers, II,’ N.A.C.A. Technical Report No. 30 
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interest to note that the drag of a propeller is, in general, con- 
siderably greater when idling than when “‘dead.’”’ The increase 
in drag due to idling is greater as the pitch decreases. 


25 


es) 
Figure 88. 


of Pitch and Blade Widths 


The drag coefficient of a dead propeller varies with the 
pitch/diameter ratio and, from the data on Figure 88 has the 
value 

Cp = 1.16 — 0.44 5 (81) 

General power coefficient curves. A general curve similar to 
the general efficiency curve may be drawn for the power coeffi- 
ae . Figure 8g is a plot of the ratio et against 


64) hi (=) for a number ‘of Durand’s propellers. (C.)o is 
N. nN 0 


ment.C, = 
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: V Via be : 
the value of the coefficient at aD (|). and y = m- Figure 


89 is a folded logarithmic plot, the three lines being parts of a 
single curve, with the limits indicated. This ee is 


necessary if the scale is to be open. Values of ——— ( < C,). read from a 


large scale plot are given in Table 9. 
A similar general curve may be drawn for C’, = nC, as in 
Figure 90. This curve necessarily follows from Figure 89 and 


the general efficiency curve. Values of ass obtained from a 


large-scale plot are also given in Table o. 


n \ n Te 3 4 S 6 aid 8 2 10 
1.0 LS 20 


Figure 89. General Power Coefficient Curve 
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Figure 90. General Thrust Power Coefficient Curve 


Figures 89 and go lead to the general solution of a number of 
common problems involving the relations between P, V, and n; 
some of these will be given. 


Variation of thrust power with altitude. The general power 
coefficient curves provide a simple method for determining the 
variation of either B.Hp. or T.Hp. with altitude. It is necessary 
that a climb be made to ceiling at constant true air speed and the 
variation of r.p.m. with altitude determined. Then, at any 


altitude, the ratio (~.) vee (=). is known, consequently 
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TABLE 9. NON-DIMENSIONAL PROPELLER DATA 


V V 
nD Cs nC2 0 : nD Cz Ca an 
=) (Cao | (aCs)o | 1m (% ) (Ca)o | (C20 | 1m 
nD/ o nD] o 
.25 75.00 33.70 .450 . 80 2.300 2.200 -956 
.30 45.50 | 23.60 . 522 85 1.850 1. 800 975 
-35 30.00 17.60 . 586 .90 1.490 1.470 .988 
-40 20.50 13.20 .645 .95 1.210 | 1.210 .998 
-45 14.50 10.20 . 700 1.00 1.000 1.000 | 1.000 
-50 10.60 7.90 ES 1.05 . 830 . 830 -995 
-55 7.90 6.30 -798 Testo) .680 .670 .982 
.60 6.05 5.08 .840 1.15 - 523 - 500 -956 
-65 4-73 4.15 WO7E 1220 - 404. .370 -917 
.70 3.70 330 .908 1.25 «324 .280 . 864 
-75 2.90 2a75 - 934 | 1.30 .250 . 198 -792 


Ce oe ee 
(Cd; Gc. i 
at sea level is equal to the product of the density ratio by & aa 


o 


isknown. The ratio of B.Hp. at altitude to B. a 


at altitude and R, is the value at 


Cc. 
(Cc ae 
sea level. The T.Hp. ratio is determined in the same manner, 
by using the general curve for nC, instead of C,. 

Table 10 contains the ratios of zone 

TAP.» 

calculated in this manner from flight test data. These data 

extended up to an average altitude of about 22,000 ft.; above this 

altitude the values in the table have been obtained by careful 

extrapolation. The relation between T.Hp. and density in the 
standard atmosphere is approximately 


Fp (am 
Pu Pen eee 


where R, is the value of 


at various altitudes as 


(82) 


These data have been checked by comparison of calculated 
performance with observed flight test data for a number of air- 
planes and excellent agreement found in every case. 
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TABLE 10. VARIATION OF THRUST HoRSEPOWER WITH ALTITUDE AT 
Constant AIR SPEED (STANDARD ATMOSPHERE) 


Altitude T Hp. Po Altitude TH. y@ 
Ft. T.Hp.o Pp Ft. T.Hp.o ? 
ro) 1.000 I.000 16,000 508 1.282 
2,000 925 1.030 18,000 465 1.325 
4,000 854 1.061 20,000 425 1.370 
5,000 820 1.077 22,000 387 1.418 
6,000 787 1.094 24,000 352 1.468 
8,000 723 1.128 25,000 336 1.521 
10,000 663 1.164 26,000 320 1.577 
12,000 608 1.201 28,000 291 1.635 
14,000 .556 1.240 30,000 261 1.697 
15,000 «532 1.261 32,000 234 1.763 


B.Hp. and r.p.m. at various speeds. A general relation pe- 
tween full throttle, B.Hp., r.p.m., and V is readily determined 
from the general power coefficient curve. Holding V constant 


and varying 7, the ratio oe ~ may be calculated from the 


(oo) 
(ro), 


If this procedure is carried out, using ratios instead of definite 
values, a family of curves may be plotted as in Figure 91. Each 


Hp. 


of the steeply inclined curves represents the variation in Blip. 


where B.Hp.. is the B.Hp. at 


known ratios 


C, 


and (Cayo Ej; 


‘ ore 
with ~ at a constant ratio —, 
No Vo 
r.p.m. = nm and at the designed speed for maximum efficiency V. 
If an engine power curve is plotted on this figure, in the form of 


el. against ", the intersections with the velocity curves 


B. p-o No 

determine the B.Hp. and r.p.m. at any speed. The power curves 
from Figure 78 may be plotted on Figure 91. These curves 
represent normal maximum r.p.m. of 1,800, 2,100, and 2,400 on 
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cee uniene 


a 
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| AT 
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UA 
CAPA EEE 
CMTE TEE || ies 
MA MT /aeteeep. | | 


A 
2) 1.00 110 1.20 
Rihiboo gi. General Curves for B.Hp., r.p.m. and V 


engines which deliver maximum B.Hp. between 2,500 and 2,700 
r.p.m. The curves might also be considered to represent approx- 
imately 75%, 85%, and 95% of the r.p.m. for maximum B.Hp. 
In this manner, allowance may be made for the fact that certain 
engines are operated nearer to their maximum power than others, 
or the same engine may be used with widely different r.p.m. in 
different airplanes. 


B.H V 
Figure 92 is a plot of BH = against V, and Figure 93 isa 
r.p.m 
r.p.M.o 


~ against ue as determined by Figure gr. 


plot of 
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Figure 92. General Curves for B.Hp. and V 
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Figure 93. General Curves for r.p.m. and V. No = 2,400 r.p.m. for Curve 
A, 2,100 r.p.m. for Curve B, and 1,800 r.p.m. for Curve C 


eee 
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Figure 96. General Curves for T.Hp. and V 


T.Hp., r.p.m., and V. General relations between T.Hp., 
r.p.m., and V in throttled flight may be determined from the 


general curve for nC, by holding ” constant and varying Ui 


No Ve 
ve 
and Gc, 
n ° 


are therefore known and the thrust power ratio is 


cee si (5) 


(83) 
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Figures 94 and 95 give the relation between 7.H>., r.p.m., 
and V, enabling the solution of any of the common problems 
involving these factors. 

The relation between T.H>p. and V with full throttle is found 
from the general curves of B.Hp., r.p.m., and 4 against V. 
This relation given on Figure 96 is of great value in performance 

V 


: f Wed seth anh : aes DA poe 
timation since it directly th t f th: =>: 
estimation since it gives directly the variation o rapa V, 


T. Hp. 
h ] 
When the value of TA 


is known at the designed speed V, for 


maximum propeller efficiency, the curve of maximum thrust 
power available is obviously determined by one of the curves on 
Figure 96. 


CHAPTER VII 
PERFORMANCE CALCULATION 


The complete performance of an airplane comprises its 
maximum and minimum speeds, rate of climb at various altitudes, 
time to climb to various altitudes, and range and endurance at 
given speeds with specified fuel loads. While it is possible to 
calculate performance by purely analytical methods, a combina- 
tion of the graphical and analytical will be found more simple and 
direct. The usual method is to calculate, for horizontal flight 
at various speeds and at a given altitude, the thrust power 
required for horizontal flight and the maximum thrust power 
available. If the curves of power required and power available 
are plotted against speed, their intersection obviously will 
determine the maximum speed at the altitude under considera- 
tion. The difference between power available and power 
required at any given speed is the excess power available for 
climb, and the maximum rate of climb occurs at that speed at 
which the excess power is a maximum. 

Before going into a detailed explanation of the general 
methods for calculating performance, a brief statement of the 
basic assumptions will be given. In steady horizontal flight the 
forces acting on an airplane are: lift, drag, gravity, and propeller 
thrust. Since the motion is steady, the vector sum of the forces 
or their components in any direction, and the resultant moment | 
about any axis must be zero. The equations of motion are 
easily obtained. Vertical and horizontal components give 
respectively 


L+Tsn0-w-=% =o (84) 
D — Tcos6 = ad =0 (85) 


146 
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where 6 is the angle of pitch of the thrust axis, L the lift, D the 
drag, and T the thrust. The thrust correction T sin 6 in equation 
(84) is negligible except at fairly large angles, and it is usually 
neglected, although the effect of a moderate throttle opening 
may be a reduction of the order of 5 mi./hr. in landing speed for 
airplanes of the pursuit type. 6 is ordinarily less than 20°, so 
that cos @ = 1.0 for all practical purposes. Therefore, unless 
extreme accuracy is required, it is always assumed that L = W 
and D = T. 

With the foregoing assumptions, the equation connecting 
weight and speed is 


L=W=C1£ sv (86) 
from which 


V= W = Boy 
i Cup SOS e (87) 
2 


since Ky = Cz, = It must be noted that when the engineering 


coefficient K, is used, the density ratio must be introduced, other- 
wise the equation holds only for standard density. Since 
absolute coefficients will be used more and more in the future, 
the engineer should familiarize himself with them. It is often 
convenient to use equation (87) in the form 


ax re 
V = 29.00 (=) 4/2 (V in ft./sec., W in Ibs., 
av an ° Sin sq. ft.) 


CL 


or V= uy el a (V in mi./ hr.) (88) 


CL 
The minimum speed, or “‘stalling speed” is obviously determined — 
by using the maximum value of Cz, in equation (88). 


Power required for horizontal flight. The power required 
for horizontal flight is 
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DV : 
1 11D = V in it'/ sec. 
pr= =e (Vin tt:/see.) 
or, T. Hp., = ad (V in mi./hr.) (89) 


D being the total drag at the air speed V. The drag D can be | 
divided into two parts, each with two subdivisions. The major 
components are wing drag and residual, or parasite drag. Wing 
drag is further divided into induced and profile (or section) drags. 
The coefficient of the former depends only on the lift coefficient 
and effective aspect ratio, while the coefficient of the latter varies 
with the wing section. Parasite drag must be divided into two 
parts, the first variable with angle of attack and designated here as 
P ,, the second independent of angle of attack and designated as P 3. 
P,, includes such items as nacelles, floats, hulls, tail surfaces, 
wing section drag, and fuselages having square or rectangular 
cross-sections. P, includes struts, wires, fittings, wheels, tail 
skids, ordnance equipment, air-cooled engines, radiators, and 
fuselages having circular or elliptical cross-sections. If there is 
any doubt as to the classification of the fuselage drag, it should 
be put in P,. 

P, and P. vary directly as V’, while induced drag varies 
inversely as V?. 


Wing drag, induced. The induced drag is that part of the 
total wing drag which is due to a virtual inclination of the lift 
vector as a result of the downwash, the effect being the same as 
if the wing were operating along an upward slope. The induced 
drag is independent of wing section and varies only with lift 
coefficient C; and effective aspect ratio n. 


CriS'. Cx? 


d i p= Se = ( 
Induced drag coefficient Cp; x(kb)? sae (90) 


where S is the total wing area, b the maximum wing span, and 
k Munk’s factor for equivalent monoplane span. For a mono- 
plane k = 1.0. Fora biplane & varies with the ratio of gap/max- 
imum-span, and with the ratio of spans. 
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Figures 26 and 27 (Chapter II) give the values of & for any 
conventional arrangement. 
Equation (90) is most conveniently used in performance 
calculations if written in the form 
2W? 
OF sO Oe 


W being the gross weight. 
For standard sea-level density this reduces to 


D; = eOTee ay (V in ft./sec.) 


VON 
een. Gs ee re 
or D?;= v2 \ap (V in mi./hr.) (92) 


These equations reduce to the form D; = = for any particular 


case. 
It is often desirable to calculate directly the induced power 
required. The equations at sea-level density are: 


TTT ee ae pe OT i) (V in ft./sec.) 


550 V 
eS eae Sees f 
or T.Ap.; = Seer 2 as NEE (V in mi./hr.) (93) 
In any particular case, these reduce to the form T.Hp.; = 


Parasite drag P, —variation with angle of attack. A study of 
the drag of nacelles, floats, hulls, and fuselages, having square or 
rectangular cross-section, shows a very definite dependence on 
angle of attack. (See Figure 3, N. A.C. A. Technical Report 
No. 236.) A study of the profile drag of ten of the most widely 


: V : 
used wing sections shows that at any given value of ve the ratio 
s 


(Cp o) 
(Cp a) t) 


previously found for variation with angle is converted to average 
speed range, the agreement will be found very close, so that a 


is practically independent of section. If the relation 
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18950 1.20 1-40 1.60 1.80 


Figure 97. Factor Fy for Finding Value of Variable Parasite Drag Px at 
Any Speed 


single curve of 2 against Le is sufficient, but in most cases 
0 S 
Di VN ae 
it will be found more convenient to use F, = dD. (;-) which 
0 S 
is plotted in Figure 97. The values most frequently used are 


as follows: (") (2) me (2).(F) 


VS 

1.00 4.100 4.10 
1.05 22210 2.44 
TeLO L756 a0) 
TE 1.500 1.98 
I.20 I.350 TAOS 
130 T2L90 2.01 
I.40 I.IIO 2.18 
I.50 1.060 2.38 
1.60 12,025 2.63 
1.70 I.007 2.91 
1.80 I.000 Boy! 
2.00 I.000 4.00 
2.20 I.000 4.84 
2.40 I.000 5.76 
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By the use of F,, the accurate calculation of performance is 
considerably simplified since P,; = (P;)o- Fy, where (P;:)o is the 
value of P, (without correction for angle) at stalling speed, Vs. 
P, is usually calculated for the condition of thrust line horizontal. 

An example of the calculation of power required for horizontal 
flight will follow. 


Example of calculation of power required. Given an airplane 
of the following characteristics: 


Gross weight: W = 2,100 lbs. 

Wing area: S = 290 sq. ft. 

Wing section: R.A. F.-15 

Maximum span: Upper = 33.2 ft. Lower = 33.2 ft. 
Span ratio= 1.00 Gap = 4.65 ft. 

Ratio (gap/span) = 0.139 Chords equal. 

Span factor k (from Figure 26) = 1.125 

Effective aspect ratio = ik = 4.8 

Parasite drag coefficients: P, = .020 P; = .026 


The maximum lift coefficient of the R.A.F.-15 from wind 
tunnel tests is Cz = 1.02. The stalling speed is 


= ok 
V = 19.77 / s A = 52.6 mi./hr. 
CL max 


The induced drag from equation (92) is 


D. = 1245 (oe _ 397,000 
nim ee enti d Af) 


kb 
The basic value of P, at stalling speed is 


(Ps)o = .020gS = .020 X 002378 X (52.6 X 1.467)? 
< 200 = 41.2 Ibs. 


The value of P, at any speed is 
P, = (Pr)o: Fs = 41.2 Fy 
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It is somewhat easier to obtain the basic value by the use | 
of the dynamic pressure at too mi./hr., that is, groo = 25.58 


: Vine 
Ibs./ft.2, the formula being (P:)o = Cop, X Yx00 X S (—) é 
Using this method, the basic value of P, is 
(P2)o = .026 X 25.58 X 290 X (0.526)? = 53.4 lbs. 


The value of P, at any speed is 


P, = (Po (7) = 53-4 (G) 


The calculations for thrust power required at sea level are 
given in Table 11. 


2 


TABLE 11. EXAMPLE OF CALCULATION FOR THRUST POWER REQUIRED 
AT SEA LEVEL 


WEIGHT = 2,100 lbs. 


DRAG 
Ad 
Speed Air oes 
Ratio | Speed | Factor Parasite Induced| Total |Required| JL 
V V Fy at WD) 

Vs |mi./hr. Sea 

P, 2 Di D Level 

lbs lbs lbs lbs 
1.00 52.6 | 4.10 168 53 143 364 51.0 5.77 
1.05 55.2) |) 2.44 100 59 130 289 AZas 7.28 
1.10 BU/fosei ll Dette 87 65 119 271 41.8 776 
1.20 OZR Le e1eO5 80 77 100 257 43.3 8.18 
1.40 Giese? || raaite’ 90 105 73, 268 52.7 Teoh 
1.70 89.4 | 2.91 120 154 50 324 Gh has 6.64 
2200)))|105..2) |= 4500 164 213 36 413 116.0 5.08 
2.30 T2120) | 5229) 217 283 27 527 170.0 3.99 


Note that < has been tabulated along with T.H.,. in Table 


tr. This is not necessary, except as a general check on the 
calculations, unless the power required for two or more weights 
is to be calculated. If the change in weight does not affect the 
parasite drag coefficients, as for example with a varying fuel load, 
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L i ; ; 
the value of D at any given speed ratio af remains constant, 
Ss 


and it is unnecessary to carry out the full calculation to find the 
new Value of T.Hp.,... At any given value of - and at constant 


S 
density, the following relations hold 


“vw W. (94) 


7 
W, W, 

D; = one = Bal iie (95) 
D 


(LHP), _ ae (96) 


(T.Hp.y); 
Equations (94), (95), and (96) are often used in performance 
calculations. 


Maximum thrust power available at sea level. The next step 
is to calculate the maximum thrust power available at various 
air speeds. To do this, the engine and propeller characteristics 
must be known. Let it be assumed that the engine is the Wright- 
Hispano, Type E, developing 195 B.Hp. at 1,800 r.p.m., and 
that the propeller is a two-bladed wooden type designed for 
general service which is a compromise between best speed and 
best climb. ‘The constant K in the diameter formula of equation 
(75) will be 300. The maximum speed will be somewhat greater 
than 115 mi./hr. as shown by Figure 138. The diameter will 
therefore be 


as YY Bp. _ Oa OS La ee 
Dr Se yf 22, ae (1,800)? X 115 
For 10% slip the pitch must be 


II5 X 1.467 X 1.10 
30 


From Figure 82 the maximum efficiency of this propeller will 


IPs = O. Bitty 


V : shite , V 
pba d —. 
occur at aD 0.72, approximately, and this is the designe aD 
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The corresponding maximum efficiency from Figure 83 is about 
78%: 
Using Figure 96, the thrust power is readily calculated for 


various speeds, by finding the values of z e ~ at the correspond- 


ing values of M where Vz = 115 mi./hr., the designed speed for 
fa 


r 


maximum propeller efficiency. At this speed the value of 
T.Hp.. is T.Hp.. = B.Ap. X qm = 195 X .78 = 152.0. Since 
1,800 r.p.m. is well below the peak of the power curve for the 
Wright E engine, the 1,800 r.p.m. curve on Figure 96 should be 
used. The calculations for maximum 7.H)., at sea level are 
given in Table 12. 


TABLE 12. CALCULATION OF MAxiwum THrust Horse- 
POWER AVAILABLE AT SEA LEVEL 


Air Speed 
V 


Re T.Hp. T.Hp. 
mi./hr. Ve T.Hp.o 

50 435 633 96.3 
60 .522 -713 108.2 
7O . 608 -770 118.5 
So . 696 . 846 128.5 
90 .782 . 897 136.2 
100 . 870 -945 143-5 
110 -957 -984 149.3 
15 1.000 1.000 152.0 
120 1.043 1.013 154.0 


Graphical calculation of performance. If the values of 
T.Hp., and T.Hp., from Tables 11 and 12, are plotted against 
air speed, as on Figure 98, the sea-level performance is easily 
obtained by graphical solution. The intersection of the curves 
shows that the maximum speed is 116 mi./hr. The maximum 
difference between T.Hp., and T.H$., occurs at 7o mi./hr. where 
(T.Hp., — T.Hp., = 70 E.Hp. The rate of climb corresponding 
to any excess thrust power is 


= E.Hp. X eo (97) 
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where W is the gross weight of the airplane. In this case 
E.Hp. = 7oand W = 2,100, so that 


dh _ 9 x 33:000 


aon so 1,100 ft./min. 


This is the rate of climb at sea level, usually designated as 
‘initial rate of climb.” 


AIRSPEED ‘V-MI./HR. 
100 1Q 120 430 


Figure 98. Curves of Power Available and Power Required, Illustrating 
Graphical Solution for Absolute Ceiling 


Absolute Ceiling or the altitude at which the rate of climb is 
zero may now be determined accurately by a simple graphical 
construction. At the absolute ceiling, the airplane will be flying 
at an angle of attack near that giving minimum power. If flight 
were started at sea level at this angle of attack and continued all 
the way up to the absolute ceiling without change, the maximum 
rate of climb would not be obtained at all altitudes, but the 
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absolute ceiling would be the same as that obtained if the correct 
air speeds for best climb were used. In this climb at constant 
angle of attack, both the air speed and T.Hp. required will vary 
inversely as the square root of the density. That is 


V= IV, Po 

i) 

and PAP pet pe Go 
? 


or, when V = KV.,T.Hp.,= KT.Hp.,. Referring to Figure 
98, it will be seen that the minimum value of T.H.,,. is 41.5 at 


V.= 57-5 mi./hr. Hence at the altitude where i 5 = 2.0, 


V=115 mi./hr. and 7.Hp.,= 83. A straight line passing 
through this point B and the initial point (V = 57.5 mi./hr., 
T.Hp., = 41.5, or the origin, V = 0, T.Hp., = 0) is the locus of 
minimum power at all densities. Consider now the power 
available. 

At V. = 57.5 mi./hr., T.Hp.. = 105.3. This is the mitial 


value C. At the altitude where V = 60 mi./hr., yf | Sone 
e 


57-5 
1.043, and = = 0.918. The altitude and the corresponding ratio 
a T.Hp. 
T.Hp.o 


It is more convenient, however, to plot the ratio = a BS ~ against 


Ve = as in Figure 99, and eliminate the intermediate step. 
? 


° 


In either case the T.Hp., at 60 mi./hr. and at the altitude where 


60 1. 
fer. —— ie det dE Pee 
f= one is determine igure 99 gives T Hp. 0.895 at 


- 1.043. Figure 98 shows that 7.Hp... = 108. Hence, 
Ss 
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T.Hp.. = 108 X 0.895 = 96.7, which is to be plotted at point 
D on Figure 98. Repeating this process for various speeds gives 


V T.Hp. 

V v. T Hb. Tepes T.Hp. 
rials I .000 I .000 105.3 105.3 
60.0 1.043 -895 108.0 96.7 
65.0 T2130 .718 113.6 81.6 
70.0 T2207 588 IIQ.0 70.0 
7ERO I. 303 .490 123.6 60.6 
80.0 I.391 . 409 128.0 Koes 


Plotting these values of T.Hp. against V on Figure 98 and 
passing the curve CDEO through the points, it is found to 
intersect the line AB at V = 78.0mi./hr. This signifies that at 
the density where the air speed for minimum power is 78 mi./hr., 
the curves of T.Hp.,and T.H>., are tangent, which is the con- 


Vo 57-5 
sponding altitude obtained from standard atmosphere tables and 
charts, or the altitude may be read directly from Figure 99 where 
it has been included for convenience. This altitude is 19,500 ft. 
The points C, D, E,—0O, just calculated are each a single point 
on a curve of 7.H>., for the altitude represented by the relation 


oes e In the example given, 6 points were used to deter- 
e : to) 
mine the curve CDEO, but 3 points usually are sufficient. 
Service Ceiling is the altitude at which the rate of climb is 
100 ft./min. If the rate of climb be assumed linear with altitude, 


that is 


iC ah (98) 


then the service ceiling is given by 


Ce ne 
hs = H | el (99) 
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where Z is the absolute ceiling and C, the initial rate of climb. ' 


For the case just calculated, H = 19,500 ft.and C. = 1,100ft./min. 


1,100 — 100 


Fie | = 17,730 ft. 


hs = 19,500 | 


Time of Climb to Any Altitude. On the assumption that 


- = (C, — ah), the time of climb to any altitude / is 


ET ED * H H 
T= C, 108e E = ‘| = 2.303 C, 108r0 laa (100) 


LJ 
WS 
S/lJ 
E\< 
—!| lil 
qin 
i Lae 
Lig 
a. | 0: 
a i os 
Fie 


lO 86120 k30 
Figure 99. Relation Between T.Hp., Density and Standard Altitude 
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Altitude Climbed in a Given Time, T. The preceding equation 
may be converted to the form 


Co 
hw e- #) (101) 
for calculating the altitude climbed in time T. 


. Calculation of power required at altitude. In steady hori- 

zontal flight Cz and Cp are constant at a given angle of attack, 
if the variation in slip stream effects on wing and tail is neglected, 
and this variation is normally negligible. Then at any given 


angle of attack, C; and af being constant, the following ratios 
Ss 


hold in passing from density 9, to pz. 


V2 O; 
Wath 4/2 (102) 


Dp: ae D (103) 
(PAPp aa) Ox 
(hip. VG. (104) 


_ These equations supplement equations (94), (95), and (96); 
they are often used in calculating performance at altitude. 

Table 13 contains the calculations for V and T.Hp., at 5,000, 

10,000, 15,000, and 20,000 ft., using the sea-level data from 


Table 11. 


Calculation of power available at altitude. The decrease in 
power available with altitude is given by Figure 99. ‘This 
decrease includes the effect of changes in B.Hp.,r.p.m., and 
propeller efficiency at constant air speed. Hence, if T.Hp., is 
known for a given speed at sea level, the corresponding value 
- for the same speed at any other altitude is 7.Hp. = T.Hp.a. X 


LUD, : : ‘ 
wae where the value of TH ° for the altitude in question 
- is read from Figure 99. 

Table 14 contains the calculation for T.Hp., at various 


altitudes, based on the data in Table 12. 
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TABLE 14. Maximum Turust HorSEPOWER AVAILABLE AT ALTITUDES 


T.He. at ALTITUDE 


Air 2 hege T.Hb.0s 


mi. /hr Sea Level 5,000 ft. 10,000 ft. | 15,000 ft. | 20,000 ft. 
: : F = 0.820. F = 0.663 F = 0.532 F = 0.425 


——q—x | ——_—_q—l qKj [ue ii “oj .__ 


50 96.3 79.0 63.9 51.2 40.9 
60 108.2 88.8 7158 57.6 46.0 
70 118.5 97.2 78.6 _ 63.1 50.3 
80 129.5 105.4 85.2 68.4 54.6 
90 136.2 Gi 7, 90.3 72.5 57.9 
100 143.5 L747. 95.2 76.3 61.0 
110 149.3 122.4 99.0 79.4 63.5 
120 154.0 126.3 102.2 sehen 
te T.Hp.a 
T.Hp.a0 


Calculation of performance at altitude. If the values of 
T.Hp.,and T.Hp., for various altitudes now be plotted against 
speed as in Figure 100, speeds and climbs are readily obtained by 
the graphical solution. Following the method previously used, 
one obtains: am 


Climbing 
Altitude Maximum Maximum Speed Rate of 
Speed E.Hp. Ve Climb 
° 116.0 70.0 70 1,100 
5,000 112.8 49.0 72 779° 
10,000 108.3 30.4 74 478 
15,000 100.6 14.0 78 220 


20,000 (85) 


The rates of climb in the fifth column are plotted against 
altitude in Figure 101, and a smooth.curve drawn through the 
points. This curve intersects the altitude axis (rate of climb = o) 
at about 19,55oft. The absolute ceiling obtained by plotting 
rates of climb, therefore, checks almost exactly with the 19,500 
ft. obtained by the short method illustrated on Figure 98. 
Figure 101 shows, however, that in general it is not permissible 
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to assume a linear relation between rate of climb and altitude. 
This subject will be discussed at length in the following chapter. 

Figure 102 is a plot of maximum air speed and air speed for 
best climb, against altitude. The variation of maximum speed ~ 


“|2000 


Be 
= : Tiga 
Gees Ce, 
B: a TMZ 
Sere 
fe 0094 4K 
et ft 
LZ , = 
100 
Figure 100. ieefaact of oe Sanne and Power Rese fe Varig 
Altitudes 


it 
ANS 
a AN 


2o5 


with altitude is much the same for all airplanes, but there may 
be considerable difference in the form of curve for best air speed 
in climb. Whether the best air speed in climb remains constant, 
or increases with altitude depends chiefly on the parasite resist- 
ance and only slightly on aspect ratio, as shown in Chapter IX. 


Effect of slip stream and miscellaneous corrections. It is 
entirely unnecessary to apply any corrections other than those 
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given. In particular, the effects of slip stream on drag are either 
negligible or may better be handled by an appropriate reduction 
in propeller efficiency. It should be noted that the factor & (the 
ratio of the span of the equivalent monoplane to the span of a 


uJ 
(a) 
=) 
= 
- 
= 
< 


RATE OF CLIMB-FT/MN. 


oO 400 800 1200 
Figure tor. Variation of Rate of Climb with Altitude 


biplane) includes the effect of the factors formerly listed as 
biplane interference, gap/chord ratio, and aspect ratio. The 
effect of stagger is considered negligible from a theoretical view- 
point. 

The effectiveness of any method is in the consistency of the 
results obtained. The method outlined has been thoroughly 
. tested and found far superior to any other so far developed. 
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Figure 102. Variation of Maximum Speed and Best Climbing Speed with 
Altitude 
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CHAPTER VIII 
VARIATION OF RATE OF CLIMB WITH ALTITUDE 


In the past it has been customary to assume that the rate of 
climb decreased linearly with altitude, that is, 


dh 


oF = C, — ah (98) 


and this assumption was apparently justified by practically all 
of the early test data. In 1924 the author made a theoretical 
analysis which led to the conclusion that the relation should be 
of the form 

= = —a-+ 2bh (105) 
which means physically that it is the slope of the rate of climb 
curve, instead of the curve itself, which is linear with altitude. 
Flight test data were analyzed and ample evidence of the existence 
of this relation found. Subsequently all of the accurately con- 
ducted Navy performance tests have been so analyzed and a large 
number of performances calculated for the same purpose. No 
exceptions being found, the conclusion is unavoidable that all 
future work must be on ‘the basis of equation (105). 

This does not lead to any great amount of extra work. 

Integrating equation (105) gives 

dh 


ioe C, — ah + bh (106) 


A second integration gives 


of ue ae ee eel 
Va? — 45C, loge E a — Va? — 4bC,) + 2C (107) 
165 
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I 
Letting K = Ve ae: 


C, = (-—a + Va? — 4bC,) 


and C, = (—a — Va? — 46C,) 
the equation reduces to 
hy Crh + 2G, 
T = KE log berrae a sal (108) 


At the absolute ceiling a =o. Hence 


bH? —aH+C,=0 


+a —Va? — 40C, 


and H = ok (109) 


The time to climb to the absolute ceiling is infinite and this 
occurs when the denominator of the log term in equation (108) 
is zero. Hence 


C,H + 2C,=0 
ot H = “3S = PE (110) 
The service ceiling is given by bhs? — ahs + C. = ah = 100, from 
which e 
fig Ge 


The time to climb to any altitude is readily obtained from 
equation (108) when written in the form 


ok Crh +2C, 
Czh + 2C, 


and solved for 4, giving 


T 
_ 2C, (eK — 1) 


T (112 
C, — Cz, eK 
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Values of the constants a and 6. a and b may be calculated 
from equation (105) if the slope to the rate of climb curve be 
measured at two or more altitudes, but this method is not 
practical because a slight error in reading the slope will be 
enormously magnified in the values of aand b. It is much better 


-.07 


Figure 103. Graphical Solution for a and b in True Rate of Climb Equation 


to plot the slopes against altitude and pass a straight line through 
the points. The intersection of this line with the slope axis 
(h = o) gives directly the value of a, while 0 is readily determined 
by the slope of the line. tras "9 

To illustrate by an example, the values of qh (5) obtained 
from Figure 1o1 are as follows: 


Altitude 2,000 4,000 6,000 8,000 10,000 12,000 14,000 
(—a+2bh) —.0660 —.0635 —.o610 —.0580 —.0550 —.0520 —.0495 


When (—a-+ 2dh) is plotted against altitude, Figure 103, a 
straight line passed through the points is found to intersect the 
axis at @= —.0693. To find 0, take the reading at some 
altitude say 10,000 ft. where (—a + 20k) = .o551. Hence 


b= 2093 — 055" = +.000000710. To obtain maximum ac- 
2 ) 
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curacy this operation should be performed at the highest altitude 

on the plot; in this case 16,000 ft., where (—a@ + 2bh) = —.0466. 

iapresiaes .0693 — .0466 
32,000 


the first value. 


= +.000000710, checking exactly with 


Absolute ceiling. Substituting these values, a = —.o693 and 
= .000000710, into equation (109) gives 
_ ta + Va? — 4bC, fr CONSE aa V .00481 — .00312 Be 


2b 2 X .000000710 
19,500 ft. 


156 


which checks the value indicated on Figure rot. 


Service ceiling. To illustrate the calculation for service ceil- 
ing, the same data will be used. The initial rate of climb was 
1,100 ft./min. From equation (111) the service ceiling is 


ee Goe Va? — 40(Ge =" 100) 1 £6603 7— V00481 — .00284 
a 2b i .000000710 X 2 
= 17,500 ft. 


hs 


which also checks with the value indicated on Figure ror. 


Time of climb. The values of the constants in equation (108) 
are: 


I I I 
ba REY = ee ore ew = 24. 
Va? — 40C. V.00481 — .00312  .O4tI i 
C, = (-a + Va? — 46C,.) = —.0693 + .0411 = —.0282 
C,.=(—a- Vee 4bC,.) = —.06903 — .O411 = —.1104 


The time of climb to any altitude is, therefore, given by 


2,200 — .0282 | 


T = 24, | queee sed ee 
punenes Ee — .1104h 


Figure 105 contains curves for the true time of climb curve 
as calculated by this equation. On the same figure are given the 
climb curves for the two straight-line assumptions shown on 
Figure 104. Assumption No. 1 is that the rate of climb is linear 
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with altitude from the initial climb to the absolute ceiling. This 
assumption results in an entirely too optimistic time of climb 
curve and must be rejected. Assumption No. 2 is that the rate 
of climb is linear with altitude and that the slope is determined 


K\S~ASSUMPTION NO.1 

AN ASSUMPTION ee 
\\ 

a te 
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Figure 104. Comparison of True Rate of Climb with Assumed Rates 


by the rates of climb at moderate altitudes. This assumption 
is very common in reducing performance to standard conditions, 
for obvious reasons. It leads to a reasonable approximation to 
the true times of climb, but the initial climb, absolute ceiling, 
and service ceiling are reduced 5% or more. If accuracy is 
desired, the true rate of climb curve must be drawn and the values 
of a and 6 determined for each case. 
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At this time it is desired to point out that the curvature in the 
rate of climb curve which has so far been considered, is due 
entirely to the relations between power available and power 


o' : - 
fo) 10 20 30 40 50 60 70 
Figure 105. Comparison of True Time of Climb with Times Obtained from 
Assumptions Illustrated in Figure 104 


required, and not to decrease in weight from use of fuel in the 
climb. The decrease in fuel load is normally sufficient to increase 
the service ceiling but slightly. The two effects are combined in 
the flight test data. 


Average values of aand 6. The value of a depends largely 
on the relation between absolute ceiling H and initial rate of 
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climb C., and may vary from .o4o to .or1, according to data now 
available. The value of 6 varies from about 4.0X 10 7 to 
9.0 X 10—7 with an average value around 7.0 X 1077. 


When H, and C, are given, a is readily calculated from the 
equation 


a= 4bH 


In performance estimates the value 


b = 7.0 X 10-7 may be used. 


CHAPTER TX 
ASPECT RATIO AND PARASITE DRAG 
The performance of an airplane is determined almost entirely 


by aspect ratio and parasite drag, yet few designers have more 
than a very general understanding of the relative importance of 


Beear 
4eAE 
a 


FULL. SCALE POLAR 
Co 
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Figure 106. Full-Scale Polar for VE-7 Airplane 


these two variables. Since much could be written on the subject, 

this short chapter supplies only an outline. In order to obtain 

the most striking results, the procedure adopted is as follows: 

An airplane for which the full-scale polar is available, is selected, 
172 
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and the performance calculated for a series of aspect ratios with 
the original parasite coefficient Cpp = .046. This process is 
repeated with the total parasite coefficient reduced to Cpp = .031, 
corresponding to a reduction of about one-half in the original 
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Figure 107. Curves of Power Required at Sea Level for Various Aspect 
Ratios with Normal Parasite 


structural parasite, and also with the total parasite increased to 
Cpp = .076, corresponding to twice the original structural para- 
site. ‘‘Aspect ratio,” here as elsewhere in this volume, is to be 


__ (kb)? a 


the “‘effective aspect ratio,” , as explained in Chapters 


II and VII. The parasite bee vary with Cz, and speed 

range but in referring to any condition, the constant coefficient 

at low lifts will be used. It includes the wing section drag. 
The VE-7 airplane has been selected as the basis for the 
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calculations because the full-scale polar is available,* and also 
because it lies near an average value for both aspect ratio and — 
parasite. Figure 106 is the full-scale polar as given on Figure 11 
of N.A.C.A. Technical Report No. 220. On this figure are 
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Figure 108. Curves of Power Required at Sea Level for Various Aspect 
Ratios with Low Parasite 


also drawn the curve of induced drag coefficient for the VE-7 
effective aspect ratio, 7 = 4.8, and the curve of parasite drag 
coefficient Copp = Cp —Cp;. The full-scale maximum lift co- 
efficient is approximately C; = 1.20. The wing area is 290 sq. ft. 
With a gross load of 2,230 lbs., the stalling speed is 50 mi./hr. 

The method of calculation was as follows: At a number 
of selected air speeds the corresponding values of Cz were cal- 
culated and Cpp read from Figure 106. The parasite drag and 
parasite T.Hp., were then calculated from 


IW. F. Durand and E. P. Lesley, “Comparison of Tests on Air Propellers in Flight 


with Wind Tunnel Model’ Tests on Similar Forms,” N.A.C.A. Technical Report No. 
220 (1925). 
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Dp = Cop£ SV? 
IDES NEA 
375 


The same values of 7.Hp.,p were used with all aspect ratios. 
Since induced drag was not required, the induced power required, 
was calculated directly from the equation 


222 eS a eee) Gee ,690 
T Hip = 382 (WY. -382W" _ 5.69 


and 7.Hp.,p = 


nSV nV 
n= being the effective aspect ratio, and V the air speed 
in mi./hr. At any air speed, total T.Hp., = T.Hp..p + T.HAp.,;. 
160 ; 
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Figure 109. Curves of Power Required at Sea Level for Various Aspect 
Ratios with High Parasite 
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Values of T.Hp., were calculated at various speeds for aspect 
ratios 3, 4, 5, 6, 8, and 10, with each of the parasite coefficients. 
Power curves were drawn for a series of altitudes so that speeds, 
rates of climb, and ceilings could be accurately determined. 
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Figure 110. Rates of Climb for Various Aspect Ratios with Normal Parasite 
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In these calculations the Wright E engine developing 195 B.Hp. 
at 1,800 r.p.m. was used. (See Figure 10 of N. A. C. A. Technical 
Report No. 220.) Figures 107, 108, and 109 contain the sea-level 
power curves while Figures 110, 111, and 112 contain the rates 
of climb against altitude. Let us now consider each item of 
performance. . 


High speed at sea level. It is immediately apparent from an 
inspection of Figures 107, 108, and 109, that the effect of parasite 
is tremendously greater than that of aspect ratio. When Cpp 
= .046, the high speed increases from 114.0 to 117.3 mi./hr. 
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_ Figure 111. Rates of Climb for Various Aspect Ratios with Low Parasite 
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Figure 112. Rates of Climb for Various Aspect Ratios with High Parasite 
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as the effective aspect ratio is increased from 3 to 10 and a 
similar increase is obtained at the other parasite values.. In 
contrast to the small effect of aspect ratio, the effect of changing 
Cpp from .046 to .031 is to increase the high speed (at aspect 
ratio 5) from 115.7 to 134 mi./hr. The relative effects are 
more clearly brought out by Figure 113. 

It must be understood that these curves merely show the 
effect of aspect ratio and parasite on the high speed of an airplane 
having 290 sq. ft. of wing area, weighing 2,230 lbs. and fitted 
with an engine developing 195 B.Hp. at 1,800 r.p.m. They 
do not apply to any other wing area, weight or power. 
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Figure 113. Comparative Effect of Aspect Ratio and Parasite on High Speed’ 
at Sea Level : 


a | a 


Ch. 9] ASPECT RATIO AND PARASITE DRAG 179 


[EAL 


/ 


2 
Z 
17 


° ABSOLUTE PARASITE DRAG COEFFICIENT—Copo. 
‘ 04 ‘ 


02 eX) 0S 06 07 08 


Figure 114. Comparative Effect of Aspect Ratio and Parasite on Initial 
Rate of Climb 


Initial rate of climb. Figure 114 is a plot of initial rate of 
climb against aspect ratio and parasite coefficient. This figure 
shows that. the higher the parasite, the greater the effect of 
aspect ratio. In general, aspect ratio has somewhat more 
effect than parasite on the initial rate of climb. 

The percentage changes shown on this figure may be decep- 
tive. The actual increments shown will hold approximately for 
any basic rate of climb. For example, if the power be increased 
so that the initial climb is 2,000 ft./min. instead of 1,000 ft./min., 
then all of the other climbs will be increased by approximately 
1,000 ft./min. Consequently aspect ratio becomes increasingly 
important as the initial rate of climb is decreased. 


Absolute ceiling. Figure 115 is a plot of absolute ceiling 
against aspect ratio and parasite coefficient. The effect of a 
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change in aspect ratio is to add or subtract a constant amount 
from the ceiling regardless of the parasite coefficient. Aspect — 
ratio is of more importance than parasite coefficient, and par- 
ticularly so when the normal ceiling is low. 
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Figure 115. Comparative Effect of Aspect Ratio and Parasite on Absolute 
Ceiling 
These curves, like the high-speed curves, apply only to the 
one particular wing area, weight, and power, used in the cal- 
culations. They cannot be used directly to estimate the abso- 
lute ceiling for any other weight, wing area, and power. 


Maximum 5 Values of . may be calculated from the 


relation 
L eat WV 
D 375: T.Ap., (113) 


where 7.Hp., is the total thrust power required for horizontal 
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flight at the air speed V mi./hr. Values of - were obtained in this 


manner and plotted against V in order to determine the maximum 
and the corresponding air speed. These maxima are plotted 
against aspect ratio and parasite coefficient in Figure 116, which 
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Figure 116. Variation of Maximum . with Aspect Ratio and Parasite 


is of a general nature and holds true within narrow limits for all 
airplanes having a stalling speed of approximately 50 mi./hr. 


Figure 117 is a plot of the ratio (air speed at maximum -) /(stall- 


ing speed), against aspect ratio and parasite coefficient. 


Speed for minimum power. Figure 118 is a plot of the ratio 
(speed for minimum power)/(stalling speed), against aspect 
ratio and parasite coefficient. This information is of considerable 
value in estimating maximum endurance and absolute ceilings. 
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Figure 117. Effect of Aspect Ratio and Parasite on Air Speed for Maximum f 


Speed range constant K. The writer showed in N. A.C. A. 
Technical Report No. 173,” that the speed range of the average 
airplane was given with considerable accuracy by 


Vs ( ay: (114) 


where K had an average value of 10.2 with negligible variation 
due to changes in wing loading and power loading, at constant 
aspect ratio, and parasite coefficient. As a by-product of the 
present investigation, the variation of K with aspect ratio and 
parasite coefficient has been determined and is plotted on 
Figure 119. The original value (K = 10.2) need not be changed 
for the average airplane, but the extreme cases likely to be met 


2°'W.S. Diehl, ‘Reliable Formula for Estimating Airplane Performance and the Eff 
of Changes in Weight, Wing Area or Power,” (1923). i 
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range from about 9.5 to 10.8, with a normal variation from 


about 10.0 to 10.4. . may be obtained from Figure 116. 
Initial climbing speed factor. It was shown in N.A.C.A. 
Technical Report No. 173, that the initial rate of climb was given 
by 
tel Kahin = Vs peest (Vin 7 Vs) 
Co = 33,000 Ww ee eee (115) 
Hp. 375 5 
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Figure 118. Effect of Aspect Ratio and Parasite on Speed for Minimum Power 
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where K, had an average value of about o.30 for the cases then 
considered. As a matter of convenience it was assumed that 
K,=3 so that the climbing speed Ve = Vs + 3 (Vin — Vs). 

Values of K, have been determined as a part of this study 
and are plotted against aspect ratio and parasite in Figure 120. 
It appears that K, may be expected to vary from 0.25 to about 
0.33 under normal conditions and between 0.20 and 0.35 under 
extreme conditions. 0.28 is a better average value than 0.33. 

Figure 120 may be used to estimate the sea-level air speed 
for best climb. 
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Figure 119. Effect of Aspect Ratio and Parasite on Speed Range Constant K 
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Figure 120. Effect of Aspect Ratio and Parasite on Air Speed for Best Climb 
at Sea Level 


Climbing speeds at altitude. This study will be completed 
by Figure 121 showing the effect of aspect ratio and parasite 
on the air speed for best climb at various altitudes. This and 
the preceding Figure should be of some value to test pilots in 
laying out air speeds for best climb. 

The maximum speeds at altitude will be analyzed in 
Chapter X. 
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Figure 121. Effect of Aspect Ratio and Parasite on Air Speed for Best 
Climb at Altitudes 


CHAPTER X 


REDUCTION OF OBSERVED PERFORMANCE TO 
STANDARD CONDITIONS 


Since the power required for horizontal flight at a con- 
stant angle of attack varies as i regardless of the actual 


pressure and temperature while the thrust power available 
varies approximately as P™ and T~°5, it follows that the 
performance obtained at any given pressure varies with the 
temperature at that pressure. “Reduction to standard con- 
ditions” is descriptive of the corrections that must be applied to 
observed performance in order to bring into agreement test 
data obtained under different temperature-pressure relations. 
The standard condition generally agreed on is an arbitrary 
specified variation of pressure and free-air temperature with 
altitude, designated as the “standard atmosphere.” If the 
method of reduction is correct, the reduced performance in 
standard atmosphere will be identical, regardless of the tem- 
perature-pressure conditions during the flights, providing neces- 
sarily that the temperatures are within reasonable limits im- 
posed by the engine-cooling system. 

Two general methods have been widely used in performance 
reduction. In the older of the two, the ‘‘density” method, an 
observed rate of climb or an observed air speed is plotted at 
the altitude in standard atmosphere at which the density is 
that determined by the observed pressure and temperature. 
In the other, the “‘pressure”’ method, the altitude is determined 
by the pressure only. Neither method is satisfactory, but in 
general, the pressure method is the better of the two. If the 
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temperatures do not depart more than a few degrees from 
standard, both methods give very nearly identical results, but 
as the temperatures diverge from standard, the reduced data 
become unreliable. Consider two climbs, one made with tem- 
perature 10°C. above normal, the other with temperature 
10°C. below normal at all altitudes. On the density basis, 
the reduction gives two approximately parallel climb curves, 
one starting at + 1,170 ft. and the other at — 1,230 ft. On 
the pressure basis, the reduction gives two diverging curves, 
both starting at the same point. It can, therefore, be definitely 
stated that the density method will not give consistent results, 
and that the pressure method is the better, but not entirely 
satisfactory. 

Several modifications to each method have been used at 
one time or another. A method recently proposed and investi- 
gated in England ist based on the use of an engine power-factor 
varying as P3ot and according to the British reports this method 
is superior to the pressure method. It is of interest to note that 
the relation Pip? is equivalent to PT—3 which may be compared 
with the relation P™’5 T~°5 used by the author. 

Instead of trying to fit various theoretical or empirical 
methods to observed performance, it appears more logical to 
determine the variation of performance with temperature and 
pressure and thus determine an accurate method of perfor- 
mance reduction. An extensive series of systematic perfor- 
mance calculations have been made for this purpose and very 
definite results obtained. Owing to the great number of calcu- 
lations and curves involved, it will be possible to give only a 
brief outline of the methods followed and the final results. 

In the first series of calculations, use was made of the per- 
formance data for aspect ratio 4.8 and normal parasite, as 
calculated for standard conditions in the study of aspect ratio 
and parasite in Chapter IX. Powers required and available 


‘s R. & M. No, 1080, ‘“‘The Reduction of Performance Tests on the Standard Atmos- 
phere,” by R. S. Capon (1927); and R. & M. No. 1099, ‘A Discussion of the Law of Varia- 
tion of Engine Power with Height,” by H. Glauert (1927). 
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were calculated for absolute temperatures 83.3%, 90.9%, 110%, 
and 120% normal at constant pressures corresponding to alti- 
tudes of 4,000, 8,000, 10,000, 12,000, and 16,000 ft. The parti- 
cular low temperatures used were selected on account of the 
simplification in slide-rule calculations. Powers available were 
assumed to vary as P™'ST~°5, The data were plotted on a 
fairly large scale so that maximum speeds, climbing speeds, and 
actual rates of climb could be determined accurately. 

In the second series of calculations the data on the six aspect 
ratios and three parasite coefficients considered in Chapter IX 
were used. This investigation was limited to temperatures 
83.3% and 120% normal at 10,000 ft. only, the first series having 
shown the effects to be linear with temperature and independent 
of altitude. The results of this investigation are outlined in 
the following paragraphs. 


The reduction of climb to standard. Table 15 contains a 
summary of the results of the first series of calculations at 
12,000 ft. pressure altitude. It was immediately apparent from 
an inspection of the preliminary data that a definite relation 
existed between the pressure altitude 4», the density altitude hz, 
and what may be called the ‘‘equivalent altitude” in standard 
atmosphere #. (At an “equivalent altitude,” under standard 
conditions, the rate of climb is the same as that obtained under a 
given non-standard condition.) As shown by Table 15 the 
final relation was found to be 


vas hp — K (hp — ha) (116) 


The values of K found at altitudes other than 12,000 ft. are as 
follows: 


ak 833 909 I.10 I.20 
Ts. 
Values of K at 4,000’....... 354 BE 350 362 
S5000 gatas ales 337 328 331 382 
TO,OOO( etre 336 2.31 309 377 
) 


[Ch. ro 


ENGINEERING AERODYNAMICS 


190 


o'I+ 6'0+ 
¢Clo-o 9Sf‘o 
oz9‘S — o0S6'c— 
oo1‘z— o0So‘I — 
OoI‘rI oSo‘f1 
C9z (oh 
£°61 Loce 
o'9L o'¢l 
Z°So1 9 ‘Sor 
o0z9‘Lr oS6‘v1 
LLS‘o0 of9°0 
000‘ZI 000‘zI 
00z ‘I OOI ‘I 


000°! 


C‘o— 
Zgl'o 
o10'e 


oSi‘t 
oS$g‘or 


Lby 
Zoe 
0°89 
Z ‘vor 
066‘g 
z9L‘o 
000‘ZI 


6060 


po y— = ee ee eee SUeLe) #86) 6: @e, @ 'em8: Buei6) é. OL . ne eee 
sgt oO yu ty ~ 2 
oS9'S soe eee rene i  ) (Py—Fy) 
oSi‘z Ce . . . ed (Yy—?Y) 
0S$3‘6 oe eee enews ee  Y ° “WY ‘quilg 
jo oyey enjoy JOJ “WY “puis Url opnyny 
oos 6. ig, O68 6) 0 0) Bl oh ne: uru/ aa ‘quil[9 jo a4e yy jenpy 
v . cc G11G) 8) 6 le. si wile) Swe) islie, db) syle. w) ate fale asS (xeut) ‘GH I SSOOXAL 
¢ “99 Ce ee? Jy /*rur DA ‘peeds SsuIquiyyy 
o°CSor Cr sy / "tu uf ‘pseds “XPT 
oS1‘9 CC ee Sd Py ‘apninyty Ayisued 
zg Oy ne hes eusnetle Money er 54: sito al araeueus o- od /d ‘oryey Ayisuaq 
ooo'zl ec cy cin Fy ‘apninyty ainssolg 
CeO. Olt pe wanes ser eta dwey, prepueis/eimyeiodue y, 


Se RIES WAN CM CDS CEE IER LA iin Scat a i A Sa IE Dt Bt SI, DA Ae Ed Ot ee Me ie eae 
‘Id 000'ZI Ly ANITD {0 ALVY IVALOY NO GaALVAaANAL 10 Lloaaay ‘SI ATAV], 


REDUCTION OF OBSERVED PERFORMANCE I9I 


Ch. 10] 


gI¢ 6z¢° cce Loc’ 
cee vce: 1ge cove 
gre: Col 6g¢° Ivc: 
oOg9f° TLCS 6g¢° LS¢° 
glo: oot: Zev: ac 
10v° giv: oS: 16¢° 
oz'I ££8°0 
eseloAy ST ST 
aseloAY ee = — 
[etouay L L 
oyseleg Ys 


oo¢: gic: Sgz° 
LIS: bre: b6z° 
Ife: ZC’ boc: 


Qbe- zoe" ro) toi 
gce: Bgc° [hese 
Zee’ gov" 19¢° 
Oz ‘I ££38°0 

esereay| 1p!) Sr 
pal ae ob 


ayIseleg MOT 


SINGIOMAAIOD ALISVAVG GNV SOILVY Loddsy saoravA wor 


oz'I £lg°0 
aseioAy | _ es om cE 
L cE 


= 


C9 FINO CO © 


o1ey 
qoodsy 
aATOOT A 


192 ENGINEERING AERODYNAMICS [Ch. 10 


Values of K found in the second series are given in Table 16. 
The variation in K is so slight that the value 0.36 may be used 
without appreciable error for any normal airplane. 

Two methods of plotting are available. A curve of equiva- 
lent altitude h = hy — 0.36 (4, — ha) against time may be 
drawn and the rates of climb at various altitudes determined 
from the slope of the curve, or the geometrical rate of climb 
may be calculated and plotted against h. 


The reduction of maximum speed to standard conditions. 
Referring to Table 15, it is seen that the maximum speed at 
12,000 ft. pressure altitude is decreased 1.4 mi./hr. by tem- 
peratures 83.3% normal and increased 1.0 mi./hr. by tempera- 
tures 120% normal. These two conditions correspond to density 
altitudes of 6,150 ft. and 17,620 ft., respectfully. It is there- 
fore obvious that speeds should not be plotted on the density 
basis. j 

Table 17 contains the maximum speed data for the six 
aspect ratios and three parasite coefficients at temperatures 
83.3% and 120% normal at 12,000 ft. pressure altitude. These 
two temperatures correspond to density altitudes of 4,040 ft. 
and 15,670 ft. respectively. A study of the speed data shows 
that the effect of a temperature change at a given pressure 
increases very slightly with increase in aspect ratio and with 
decrease in parasite coefficient. Low temperatures decrease the 
maximum speed, high temperatures increase it, at a given 
pressure altitude. Up to 70% of the absolute ceiling, the 
change in high speed for normal temperatures is less than the 
experimental error, so that actual measured true air speed may 
be plotted, without correction, against the pressure altitude. A 
considerable error might easily be possible at altitudes above 
70% of the absolute ceiling, if the measured speeds are followed 
too closely. The relations are such in this region that in all 
probability no satisfactory method of exact reduction can be 
found. It is recommended that this part of the data be ob- 
tained by extrapolation, to the absolute ceiling, of the curves of 
maximum speed and climbing speed plotted against altitude. 
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Example of climb reduction to standard. Tables 18 and 19 
contain the observed data and calculated altitudes for two 
climbs on a typical airplane. In climb No. 1 the ground tem- 
perature was 26.5°C. and at all altitudes the temperature was 
higher than normal. In climb No. 2 the ground temperature 
was 12°C. and the temperature at altitudes was sometimes 
below, sometimes above normal. The altitudes h, calculated from 
equation (116), are plotted on Figure 122 and the rates of climb 
determined by the slopes of the curves. In spite of the con- 
siderable difference in initial conditions, the rates of climb are 
in excellent agreement when reduced by the new method. 
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Figure 122. Comparison of Reduced Performance from Two Climbs by 
Different Pilots 
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TABLE 18. REDUCTION OF CLIMB TO STANDARD, Ciimp No. 1 
OBSERVED DATA ALTITUDES (ft.) 

. Pressure} Air . 
Time Pressure |Density 
Min. ae toe a he hp —ha| Ah h 
re) 762 26.5 —70 1,170 |—1,240 | —446 +376 
2 697 18.5 Bafa aigaOet — 900. B45 et ia. 7 35 
4 642 12.0 4,590 5,290 —700 —252 4,842 
6 594 8.0 6,660 7,430 | —770 —277 6,937 
8 563 5.0 8,070 8,780 |, —710 —255 8,325 
10 533 3.0 9,490 | 10,270 —780 —280 9,770 
12 509 0.5 10,680 | 11,460 —780 —280 | 10,960 
14 493 —3.0 11,490 | 12,050 | —560 —202 | 11,692 
16 477 —6.0 12,330 | 12,730 | —400 —144 | 12,474 
18 458 =7;,0 13,350 | 13,890 | —540 | —194 | 13,544 
20 446 Oi 14,020 | 14,470 —450 —162 14,182 
25 421 |—I10.0 15,450 | 16,140 —690 —248 | 15,698 
30 ZO3' 5 — 10 16,230 | 16,950 —720 —259 | 16,489 
A h = 0.36 (hp — ha) h=hp — (Ah 
TABLE 19. REDUCTION OF CLIMB TO STANDARD CLimB No. 2 
OBSERVED DATA ALTITUDES (ft.) 
ame ike heats Bee hehe | Ak h 
in. Hg. on p 
re) 764 12.0 —140 —520 | +380 +137 —277 
3 670 6.0 || 3,440 | 3,180 | +260 $030) 3847 
4 645 5-9 || 4,470 | 4,340 | +130 +47 | 4,423 
6 605 5.0 6,170 6,470 —300 —108 6,278 
8 565 1.0 7,970 8,160 —190 —68 8,038 
10 531 mono 9,580 9,700 | —I120 —43 9,623 
12 507 —3.0 || 10,780 | 11,150 | —370 —133 | 10,913 
16 | 464 —6.5 || 13,030 | 13,330 —300 —108 13,138 
18 447 —9.5 || 13,960 | 14,330 | —370 —133 | 14,093 
20 435 —10.5 || 14,640 | 14,950 | —310 It | 14,751 
25 414 —15.5 || 15,870 | 15,980 | —I10 —40 | 15,910 
30 398 —17.5 16,830 | 16,950 —120 Ag 16,873 
35 388 —19.5 |} 17,460 | 17,460 fe) Oo | 17,460 
40 381 —20.5 || 17,900 | 17,880 20 7 nL 7,007. 
Ah = 0.36 (hp — ha) h=hp — Lh 


CHAPTER “XI 
NOTES ON FLIGHT TESTING 


The following notes on flight testing have been prepared for 
the information of test pilots. Some of the material is not new 
but appears desirable in explanation of the principles involved. 
Other parts, such as the conditions for best climb and methods of 
specifying these conditions, are not available elsewhere. 


Calibration of air-speed indicators. In general, air-speed in- 
dicators must be calibrated by runs up and down wind over a 
measured course. Special methods, such as calibration against 
the reading of an instrument suspended well below the airplane* 
or against previously calibrated readings in a second airplane, 
are not available except at well-equipped flight test organizations. 
Calibration runs should not be attempted when the cross course 
component of the wind exceeds 14% (4h) of the airplane speed, 
since this value reduces the measured speed 1%. 

Assume that a series of runs in pairs has been made over a 
measured course, each pair consisting of a run up and down 
wind, during which the following data are recorded: time over 
course, air-speed indicator reading, r.p.m., air temperature, and 
pressure. The time over the course gives the ground speed. 
Averaging ground speeds (not the times over course) in each 
pair of runs gives the true air-speed for the given r.p.m. and 
indicator reading. The indicator reading is a function of 4V? 
and V is the true air speed. The curve of readings plotted 
against eV? is the calibration curve. ‘True speed may be 
obtained from the instrument reading and this curve, when ¢ is 
known. While the true speed is determined by the density, 
the speed itself must be plotted against pressure altitude. . 


7 See W. G. Brown, ‘‘ Measuring an Airplane’s True Speed in Flight Testing,” N. A.- 
C. A. Technical Note No. 135 (1923). . : i: seine ee 
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Maximum speed. In order to determine the actual maxi- 
mum speed in horizontal flight, care must be taken to avoid 
errors due to following sources: (1) failure to maintain hori- 
zontal flight, (2) starting over a speed course before the air- 
plane has settled down to steady flight, (3) variable winds, or 
cross wind component too large. There are many other pre- 
cautions to be observed, but those enumerated constitute the 
chief sources of error. 

A change in altitude of 40 ft./mi. will give 1% change in the 
measured high speed of an average airplane. For this reason 
high-speed runs are usually made at a very low altitude and, if 
available, a statoscope is carried, in order that the flight path 
may be maintained horizontal. If the ratio of maximum to 
stalling speed is greater than 1.5, the error due to ground effect 
is negligible. 

All high-speed runs should be started at a distance from 
the first’ marker on the course which will allow from 30 seconds 
to 60 seconds for settling down to steady conditions before 
crossing the line. This precaution is very important. The last 
two or three miles of high speed are attained by a comparatively 
slow acceleration, and any control movements or changes of 
altitude and direction will affect the high speed. In the first 
stages of the approach to the starting line of the course, the 
pilot should fix his altitude and direction. Under no conditions 
would a diving start be permissible. 

High-speed runs should be made in a calm if pers onic 
although good results can be obtained in any steady wind which 
does not have a cross course component greater than 7 of 
the speed to be measured. Tests should never be made with 
gusty, variable, or unsteady winds, no matter now light. 


Importance of correct air speed in climb. Figure 123 gives 
the calculated variation of rate of climb with true air speed and 
altitude for a typical airplane. Each heavy curve gives the 
variation of rate of climb with air speed at constant standard 
altitude. The central cross curve (long dashes) passes through 
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Figure 123. Variation of Climb with Air Speed and Altitude, Showing 
Necessity for Making Climb at Correct Air Speeds 


the maximum rate of climb at each altitude and therefore gives 
the variation of best climbing speed with altitude. The two 
outside cross curves (short dashes) are marked 95%, and pass 
through rates of climb 95% of maximum at all altitudes. These 
curves show very strikingly the importance of correct air speed 
in climb. ; 
Variation of best climbing air speed with pressure and 
temperature. It is very important that the test pilot know the 
proper variation of climbing air speed with altitude. How to 
express the variation of climbing air speed with altitude so as to 
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allow for variation from standard atmosphere conditions, has 
long been an insistent question. The answer to this question 
now appears simple. 

In the study of performance reduction, the best climbing 
speeds for a typical airplane were determined at various pressure 
altitudes for temperatures 83.3%, 90.9%, 110%, and 120% 
normal. A large effect was found as follows: 


Temperature; % Normal... 22. ...024.5. 83.3 7.90.9 100 110 120 
4,000". . .66.0 67.5>— 70109) 7075" 27 h-0 

eats . 8,000’...66.0 68:0 70.0 73.0 74.0 
east ond A tie Sas 10,000’...65.0 68.0 70.0 72.0 74.0 
P 12,000’...66.5 68.0 70.0 73.0 76.0 
10;000/4-150 720" 70,0), 972.0; 7760) 97020 


Taking the ratio of these climbing speeds to the correspond- 
ing stalling speeds gives: 


memperature, 9, Normal... .......2003 3263 ORO 100 IIo 120 
LAGOON Serhan Nag Ane) Pee ine) 

Ratic Best climbing speed SOOO se Tee On ein 20 = 2 24= el oe Oo 
atio Stalling speed 10,0007. Gib 220) 2522 F-02018 (48 Pasi 
at given pressure altitude 12,000". eet OU ER tOus ot. D7 ET Ooee a Keo 
C000 Deion sles 5.) tar 4An 1.203 1-13 


It is obvious that for all practical purposes the ratio of best 
climbing speed to stalling speed may be considered constant 
at any given pressure altitude. This means that the relation 
between angle of attack and pressure altitude is not affected by 
temperature variations. Consequently, the conditions for best 
climb at altitude must be specified as the relation between the 
readings of a sensitive angle of attack indicator and an aneroid. 
The common method of specifying the indicated air speeds at 
various indicated altitudes may give erroneous results. It is 
satisfactory, however, to use indicated air speeds and aneroid 
pressures. 

Determination of best climbing air speeds. The procedure 
formerly followed was to make a series of ‘‘saw-tooth”’ climbs 
at various air speeds at a given altitude and thus determine the 
air speed giving the greatest rate of climb. ‘This was repeated at 
various altitudes and a curve of climbing air speed plotted 
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against altitude. Such a curve holds true only for the particular 
pressure and temperature conditions existing during the saw- 
tooth climbs. In order to obtain a general relation, it is neces- 
sary to specify either angles of attack or indicated air speeds 
in terms of aneroid pressure. 

The recommended procedure is to determine the angle of 
attack, or indicated air speed, which gives the greatest rate 
of climb at each aneroid reading, and plot this against aneroid 
reading. For example, the angle of attack that gives the least 
time required to climb between aneroid readings of p = 510 mm. 
and P = 490, would be taken as the best angle of attack for 
P= soo. The actual increment in pressure used to determine 
the best reading should be so varied according to the rate of 
climb as to give a minimum time interval of at least 30 seconds, 
and preferably more. Otherwise experimental errors may affect 
the results. The desirable increment in pressure will vary with 
altitude (or pressure) but it is constant in all of the runs at a 
given pressure. 


Climbing tests. After the relation between angle of attack 
and pressure altitude has been determined, the actual climbs may 
be made. Under average conditions, the best method of start- 
ing a climb is to fly horizontally in a convenient direction and 
at a low altitude, long enough for the airplane to settle down to 
steady conditions at the proper angle of attack for the existing 
ground pressure. The time required is determined by the lag 
in the aneroid and air thermometer, and the climb should not 
be started until these instruments record ground conditions. At 
the proper time, the climb is started by opening the throttle 
while holding the angle of attack constant. 

The climb should be made with a minimum number of turns, 
which should be rather wide for best results. It is advisable that 
a climb be repeated by one or more different pilots, if possible. 
The composite results of three climbs should be quite definite. 


CHAPTER XII 


RANGE AND ENDURANCE 


- 


Maximum range. The maximum range for a given fuel load 
is obtained by continuously flying, until the fuel is exhausted, at 
the air speed giving the most miles per pound (or the minimum 
pounds per mile). For a given airplane this air speed, usually 
designated as the ‘‘most economical air speed,” or “‘economical 
speed,” varies almost linearly with gross weight. The fuel 
consumption in Ibs./mile also varies almost linearly with gross 
weight. It therefore follows that the most economical speed and 
corresponding fuel consumption need be calculated for only three 
gross weights in order to determine the maximum range. These 
weights should obviously be at, or near, full load, half-fuel load, 
and no-fuel load. 

Fuel consumption is obtained at a given speed and weight by 
calculating thev.p.m. and B.H>p., using the data given in Chapter 
VI. The various steps in the calculation may be tabulated as 
follows: 


1. Air speed V. 


2. Ratio Va = propeller design speed. In absence 
d 


of specific data, assume Vg = Vy. 
3. T.Hp.required at speed V. Read from curve of T.Hp,,. 


ivi tred pe © ES 
4. Ratio TH. T.Hp.. = T.Hp.y at speed Va. 
N r.p.m. 


5. Ne ae Read from Figure 94 or 95. 


6. 7.p.m. = = X 7.p.m.o. 1.p.m.. designed, or rated 7.p.m-, 


at speed Vg. 
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1) anaes 
ND ND oe 4 N° 


8. ‘e from Figure 83. 


Fs 


9 N=mxX - \m is the maximum propeller efficiency. 
m 


10. B.Hp. = ee 

11. Specific fuel consumption C. Obtain from engine data 
or from Figure 79. 

12. Fuel consumption, lbs./hr. = B.Hp. XC. 

BAP. XE 


13. Fuel consumption, lbs./mile = V7 


The fuel consumption in lbs./mile is plotted against speed 
in order to obtain the minimum consumption and the corre- | 
sponding speed (most economical speed). This operation is 
repeated for the other gross loads. ‘The minimum fuel consump- 
tions in Ibs./mile and the corresponding economical speeds are 
then plotted against gross weight. If there is a pronounced 
curvature in the line passed through the minimum fuel consump- 
tion points on this plot, the range must be determined either 
by a step-by-step integration, or by a planimeter. In the step- 
by-step method the average fuel consumptions are determined 
for convenient increments in weight, say 500 lbs. or 1,000 lbs. 
The distance flown for the AW under consideration is obviously 

AW Ibs. 
AR= aoe aaa 
crement AW. The endurance may be calculated along with 


where C, is the average lbs./mile for the in- 


the range from AE hrs. = = 


, where V, is the average speed 


for the increment AR. 

If it is desired to use a planimeter, the reciprocal of the mini- 
mum fuel consumption, that is, the maximum miles per pound 
of fuel must be plotted against gross weight. The area under 
mi 


-X lbs. = miles, and the 


this curve has the dimensions of ib 
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distance flown while using a given weight of fuel is proportional 
to the area under the curve and between the two corresponding 
gross weights. If the curve is plotted to the scales of 1” = 1,000 
Ibs. weight, and 1” = 0.10 (mi./lb.), each square inch under the 
curve will represent 100 miles. The endurance may be found 
in the same manner by plotting (hrs. endurance /Ib. fuel) 
against gross weight and measuring the area under the curve. 

The range may be calculated very quickly by formula when 
the fuel consumption in lbs./mile plots as a straight line against 
gross weight. If the slope of the fuel consumption curve, 
Ibs./mile, against weight is 0, and the equation of fuel consump- 
tion C=a+6W, then 


ut Wao ere 
ss (a +oW) 
W. 


° 


pale logro be + al miles (117) 


b a+bW, 
Since (a + 6W.) = initial fuel consumption, C. (lbs./mile), and 
(a + bW;) = final fuel consumption, C, this becomes 


ian BE 208 Cy lbs. /mile 
hn ee ie ibs. /mnile os 


The endurance is found directly from the fuel consumption in 
Ibs./hr., if C is of the form (m + nW): 


be Poe hay. 
= (m + nW) 
W, 


2.303 hte & +n2W, 


nN m+nW, 
305 C3 lbs. / hr: 
ge OE E (payor ee) 


Theoretically, the specific fuel consumption should include 
both gas and oil consumption, and the final weight should be 
based on the amount of oil remaining when all of the gasoline is 


204 ENGINEERING AERODYNAMICS [Ch. 12 


gone. Practically, the oil consumption may be neglected when 
modern engines are used and the final weight Wx taken as Wo 
less the weight of gasoline W,. That is, W. = W. — Wg. 


Calculation of maximum endurance. Maximum range is ob- 
tained by flying at a speed about 40% greater than the stalling 
speed, but maximum endurance is obtained at a considerably 
lower speed, 15% to 20% greater than the stalling speed. The 
calculation for maximum endurance is exactly like that for 
maximum range except that the factor now desired is the maxi- 
mum (hrs./Ib. fuel), at a given gross weight. Plotting (max. 
hrs./Ib. fuel) against gross weight, the area under the curve will 
be proportional to the maximum endurance in hours. 

The maximum (hrs./lb. fuel) should be obtained by calcu- 
lating (hrs./lb.) at a number of speeds from 1.05 Vs to 1.50 Vs 
and plotting against speed. 


Range and endurance formulas. Breguet’s formulas for 
range and endurance are: 


Range = 863.5 z. a logro = miles (120) 


* Ven Ley I I 
Endurance = 750 Ve D ea — S| hrs, (421) 


where 5 is the ratio of lift to drag for the entire airplane, 4 the 


average propeller efficiency, C the average specific fuel consump- 
tion, W. the initial gross weight, W, the final gross weight, and 
Vc the cruising speed at the gross weight W,. 

The author has shown’ that the range is given quite accurately 


by 
LG W, 0-45 
Riso pata | (122) 
when the mixture control is not used, and by 
L W, 0.42 
R = 892.9 aa (7) | (123) 


«W. S. Diehl, ‘‘ Three’ Methods of Calculating R: d End f Ai wy 
N. A.C. A. Technical Report No. 234 (1925). pie debarmiaiy aie Cia 
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when mixture control is used. The corresponding equations 
for endurance in hours, are 


cane ee | W,\ °° 
T = 7,500 Roy ge. re = | (124) 
when mixture control is not used, and 
rs eS W,\:°8 
T = 4,687.5 cD val (Ww) =e (125) 


when mixture control is used. In equations (122) to (125), C 
is the initial instead of the average, specific fuel consumption, and 
Vc is the initial cruising speed in mi./hr. 


1.4 I8 20 22 24 26 28 30 


Figure 124. Variation of Specific Fuel Consumption with Speed Ratio 


Ranges and endurance calculated by these formulas are in 
excellent agreement with step-by-step integration, provided 
that the proper values of C and C, are selected. Figure 124 en- 
ables the values of either C or C, to be calculated accurately. 
To calculate the average specific fuel consumption for a given 
fuel load, find the average ratio of maximum to minimum speed _ 
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i and read the value of < from Figure 124. C, is obtained 
Ss 0 


from the ratio ra corresponding to the initial value of iH : 
0) S 


240 
200) 
160 
120 

80 


40 


RATIO -LVEL LOAD — We 
TOTAL LOAD Wo 


0 10 20 .30 “40 i) £0 30 
Figure 125. Factors for Simplified Range Formulas 


Equations (120) to (125) may be written 


pane eed RAE 
R= Brey5 co (126) 
re Ege t 
P= Been ies (127) 
L ‘ 
R= Kr: Dd’ a (No mixture control) (128) 
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Set ee , 
T=Ke eye ige (No mixture control) (129) 
R= K’r 5a (With mixture control) (130) 
Poe ee” ity ciztire Gontnol) 2) (Ea) 
DC. Va ‘ _ 


The factors Br, Bz, Kr, Kz, K'p, and K's will be found on 
Figures 125 and 126. 


440 


s/n 
At 


Figure 126. Factors for retina Endurance baie 
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: : L ; 
In using these equations, the maximum D and maximum 


propeller efficiency may be used if maximum range is desired. 
The specific fuel consumption at full throttle C. is obtained either 
from engine test data or estimated from general engine data. 
The oil consumption may be included in C, if it is so desired, 


24000 


20000, 


tg000 


Figure 127. Graphical Solution for Bombing Range 


but the final weight W, must then allow for the oil remaining 
when all of the gasoline is exhausted. The difference in range 
obtained by allowing for oil consumption is hardly worth the 
trouble if the engine is average and modern. Some of the early 
engines used enough oil to make an appreciable difference. 


Bombing range. The distance that an airplane can fly with 
a given fuel load, drop a bomb, and return to the starting point, 
is easily obtained by a simple graphical solution. Let Wo be 
the initial gross load including full bomb and fuel loads, and W, 
the final load without fuel or bombs. Calculate the range for 
two or three fuel loads assuming constant initial load = Wo. 
The largest fuel load should be somewhat more than 3 (W. — W,). 
Repeat the process assuming constant final load = W,. Plot 
as two curves of range against weight, one starting from W. 
and showing increasing range with decreasing weight, the other 
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starting from W, and showing increasing range with increasing 
weight. The intersection of these curves gives the maximum 
radius with all fuel and no bombs. Figure 127 illustrates how 
the range is found for any desired bomb load for an airplane 
weighing 24,000 lbs. and carrying 10,000 lbs. of fuel and bombs. 


700 


Figure 128. Effect of Wind on Most Economical Air Speed 


For 10,000 lbs. of fuel and no bombs the curves intersect at 
1,110 miles; hence the range is 2 X 1,110 = 2,220 miles. Fora 
4,000-lb. bomb load and 6,000 lbs. of fuel, the airplane can fly 
out to the point A at 660 miles, drop the bomb load AB = 4,000 
Ibs. and return to the point of departure. 


Effect of wind on range. The most economical speed is af- 
fected considerably by wind direction and force, being increased 
by head winds and decreased by following winds. The most 
economical speed is readily found for any wind velocity, by a 
simple graphical solution. If the fuel consumption in Ibs./hr. 
is plotted against air speed in mi./hr., with both scales starting 
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Figure 129. Effect of Wind on Fuel Consumption for a Typical Flying Boat 
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Figure 130. Effect of Wind on Most Economical Speed for a 
Typical Flying Boat 
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from zero, as in Figure 128, then a line drawn from any point 
on the speed axis is tangent to the curve at the most economical 
speed for the wind condition determined by the starting point. 
Since the ground speed is the difference between the air speed and 
the wind speed, the starting point or origin for ground speed is 


ADIUS WITH WIND Vw 
RADIUS IN CALM 


RATIO 


WIND VELOCITY 


1 _Vw 
RATIO 37ERAGE ECONOMICAL SPEED = Ve 


% 40 20 30 40 50 40 70 80 ete) too 
Figure 131. Reduction of Range in a Steady Wind 


moved to the right or to the left for head winds and following 
winds respectively. That is, with a 4o mi./hr. head wind, the 
ground speeds will be 40 mi./hr. less than the air speeds, and 
the origin will be at 40 mi./hr. A line drawn from this point 
in Figure 128 is tangent to the curve at 99 mi./hr., which is the 
economical speed for a 40 mi./hr. head wind, at the gross weight 
represented by the fuel consumption curve. Similar tangents 
drawn for zero wind and 4o mi./hr. following wind indicate 
economical speeds of 86 mi./hr. and 81 mi./hr. respectively. 
It is obvious that this operation determines the minimum pounds 
of fuel per ground mile. 

It is somewhat longer, but perhaps more accurate, to calcu- 
late the lbs. fuel per ground-mile for various air speeds and 
wind speeds. If the Ibs. fuel per ground-mile at a given wind 
speed be plotted against air speed, the most economical air speed 
is accurately determined. Figure 129 shows the variation in 
fuel consumption and Figure 130 shows the variation in the 
most economical speed for a typical large flying boat. 


212 ENGINEERING AERODYNAMICS [Ch. 12 


A head wind of a given velocity reduces range more than a 
following wind of the same velocity increases it. The maximum 
distanee that can be flown out and return to starting point is 
less in a steady wind than in a calm, and decreases rapidly as 
the steady wind velocity increases. The radius is greater if 
the head wind is encountered going out with full load than if 
the head wind is met on the return flight. Figure 131 gives the 
average effect of a steady wind on radius of action. 


CHAPTER XIII 
SPECIAL FLIGHT PROBLEMS 


There are a number of problems of considerable interest and 
some practical value connected with special flight conditions. 
For most solutions extreme accuracy is not required, and reason- 
able simplifying assumptions are allowable. A few of the more 
important problems in this class will be outlined briefly. 


Gliding flight without power. For equilibrium in a glide 
without power, at an angle 6 to the horizontal 


W sin = D (132) 

Lcoos§+ Dsin§ = W (133) 

and © L.sin 0 = Dcos®@ (134) 
From equation (134) 

= = cot 0 (135) 


or the angle of glide is equal to the angle whose cotangent is a 
Since L = K ay V? and D = Kp ake, it follows from equation 


(133) that the speed along the flight path is 


va y/2 4/ = (136) 
° Kp (& cos 6 + sin 0) me 


For the limiting speed in a vertical dive Kz cos 6 is negligible 
and sin @ is substantially unity. Hence 


DMS Ae 


Vi / : i Kp (137) 
213 
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The altitude lost per minute in a glide without power at a 
fixed angle of attack must be such that the rate of change in 
potential energy is equal to the power required for horizontal 
flight at the same angle of attack. That is 


dh _ 33,000 1 .Ap., 


Veee ey W 


ft./min. (138) 
From these equations it is possible to solve practically all 
of the problems of gliding flight. 


Gliding flight with power. Assuming that the thrust T acts 
along the flight path and letting the angle @ between the flight 
path and the horizontal be positive upward, the equations of 
equilibrium are 


T—Wsind =D (139) 
Lcos§ — (D —T)sin6 = W (140) 
and Lsin§@ = (D — T)cos8 (141) 


The inclination of the thrust to the flight path may be 
allowed for, but little is to be gained in accuracy. 

Assuming that the thrust varies with velocity and that the 
drag is D = KV’, it can be shown * that the velocity along the 
flight path is 


V= Var |b +4/-2.25 = (5) sin | (142) 


if the thrust at V = 2Vy is zero. This condition applies to 
propellers of low pitch/diameter ratio. For propellers of high 
pitch/diameter ratio the thrust holds up better. Assuming 
that T= T./2 when V= 2Vy, (T= thrust at V= Vy,) the 
velocity along the flight path is 


v= Vul — 4+ 4/1.5605 - (&) sin 5625 — ( a 5), sin 0 | (143) 


I W.S. Diehl, ‘‘ The uae of Flight Path Inclinati 
maa Bast ee a g ination on Airplane Velocity,’ N. A. C. A. 
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In these equations (=) is the value of “ for the airplane 


Ds D 
at the maximum horizontal speed Vy and has the value 
5) = Vu W 
(5 0 375 BHp. Be) 


(=) has no direct relation to the maximum value of s. Figure 


132 isa plot of i against 6 for various values of (5) according 
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Figure 132. Air Speed in a Dive with Power On. Thrust Zero at V = 2 Vu 
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to equation (142), and Figure 133 is a similar plot according to 
equation (143). The effect of thrust on diving speed is com- 
paratively small. 


BSE 
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Figure 133. Air Speed in a Dive with Power On, T = = at V =2Vm 


Circling flight. For equilibrium in a horizontal turn, the 
weight must be balanced by the vertical component of the lift 
and the centrifugal force must be balanced by the horizontal 
component of the lift. If 9 is the angle of bank, the equations 
for equilibrium are 


W = Locos 9 (145) 


and centrifugal force = 7 = Lsing (146) 
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2 


from which tan 9 = a (147) 


V and r must be consistent units, e.g., ft./sec. and ft. Since 
W is constant, L must vary as sec 9 and the acceleration in a 
turn is 


L 

w= se? (148) 
Since L varies as V’, the relation between the speed in straight 
horizontal flight and the speed in a turn at constant lift coef- 
ficient or angle of attack is 


Vi 


vy = sec? (149) 


At constant lift coefficient, or angle of attack, the drag will vary 
as V? and the thrust horsepower required as V%. If curves of 
T.Hp. available and T.Hp. required are plotted against velocity 
on logarithmic scales as in Figure 134, the solution for minimum 
radius of turn at any given lift coefficient is made very simple. 
The power required varies as V3, along straight lines AA’ and 
BB’. The intersection of such a line with the curve T.Hp., 
gives the maximum air speed at which turning flight can be 
made with the initial lift coefficient. Consider the line BB’. 
B is at 100 mi./hr. and B’ at 192 mi./hr. Hence the angle of 


bank is 
- Gy = ey - 
seco = ie = CTA Peay 3.68 
or 9 = 74° — 13’ 
From equation (147) the radius of the turn is 


a BOT Gy 382 as 605 ft: 


g 32.2 

This process may be repeated for as many points, lift coef- 

ficients, speeds, or angles of attack, as desired. The results 
would be as tabulated on page 219. 
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Figure 134. Graphical Solution for Speed in a Turn 
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ie Air Speed in mi./hr. | Angle of 
Ang e _ Seah ere | Bank Radius of 
Seo Horizontal Ina ? Tura 
egrees Flight Turn deg.-min. fe 
I 198.2 250.0 51 — 06 3,370 
2 166.8 244.7 625-22 2,100 
4 130.3 227.5 70 — 52 1,200 
6 III.o 207.5 es <P Phi 861 
8 98.5 190.0 74 — 25 673 
10 89.2 171.5 (4A 17 Soo 
12 $227, 152.0 42-48 478 
14 77-4 132.3 69 — 58 427 
18 73.8 103.8 59 — 42 422 
20 | 73-5 93.6 559 457 


If the time required to turn through a given angle is calcu- 
lated, it is found that the minimum time corresponds with the 
maximum bank. The foregoing data are based on the wind 
tunnel tests of a 1924 racing airplane. 

‘“‘Squashing”’ in a turn is due to an attempt to fly under a 
condition requiring more power than is available. For example, 
any shorter radius of turn than 422 ft. at « = 18°, in the case 
calculated above, will require more power than is available and 
the airplane will lose flying speed and ‘‘squash”’ from centri- 
fugal force since the centripetal force cannot be obtained from 
the reduced lift. The same condition exists in a pull-up from 
a dive. The lighter the load and wing loading and the greater 
the power, the more difficult to “squash” on a turn. High 
aspect ratio and low parasite resistance give low minimum power 
required and therefore reduce the tendency to “squash.” 


Spiral gliding flight. The spiral glide is a combination of 
gliding and circling, the path of the center of gravity being a 
true helix. If the angle of bank is 9 and the helix angle is 6, 
then assuming that the radius 7 is large in comparison with 
the span, so that velocity over wings is substantially uniform, 
it follows that 


Lcos¢gcos8 + Dsin0 = W (150) 
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W sin 9 = D (151) 

and W cos 6 = Lcos@ (152) 
, L cot 6 

from which Ta Cee (153) 


which reduces to 5 cot 6, equation (135), when the angle 


of bank 9 = o, that is, in a straight glide. 

The proper angle of bank is determined as in simple circling 
flight by the balancing of the centrifugal force with the radial 
component of lift. That is, 

WV? cos? 0 


centrifugal force = Bee ae Lcos 6 sin 9 (154) 


I 
°. 


which reduces to equation (148) when 6 


Effect of a diving start on speed over a measured course. 
If an airplane is dived at high speed and pulled out into level 
flight, it will have a speed in excess of that which can be main- 
tained in horizontal flight. If a high-speed flight is made over 
a measured course after a diving start, the speed will asymptoti- 
cally approach the normal maximum, and the average over 
the course will be greater than the normal maximum. This 
method was formerly used extensively in races, and sometimes 
by test pilots in order to get the maximum speed possible. 
Assuming that the propeller thrust is constant and the 
flight over the course is at constant altitude, the effect of a 
diving start on the average speed may be calculated. The 
speed at any time ¢ is 
C,e*# + C, 
aes pee Ts (155) 


and the distance flown is 


Ww W ad 
Sa 7K loge (C; e% - C,) =a Vile gK loge (2VTK) (156) 


2 W.S. Diehl, ‘‘The Effect“of a Diving Start on Airplane Speed,” N. A. C. A. Tech- 
nical Report No. 228 (1925). 
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where V, is the speed (ft./sec.) at the time of crossing the 
starting line (¢ = ¢,), V the speed at time #, S the distance flown 
in ¢ seconds, W the gross weight, T the thrust at normal maxi- 


mum speed Vy, and K the drag coefficient (K = ae The 
other constants have the values 
C.=T+V.VTK 
C= V.VTK —T 
C;=VITK+4+KY, 
C,= VTK — KY, 
2g VTK 
W 
The persistance of excess velocity depends largely on the 
value of a which does not vary over as wide a range as might 
be expected. Representative values for various types of air- 
planes are: 


I 


and a 


Type W Ibs. V uw mi. /hr. a 
RAGIN D a recen et een vspiesciorers: Siakets 2,100 250 .050 
PUIESULE Rime ies teens ple lens So eh GH 25 O00 160 .046 
Obsenvatioutac: sade csns 2,300 125 .076 
SOMMER yp elute sce s Warnvesare 6S eso 7,000 120 . IOI 


A racing or pursuit airplane requires about two minutes to 
settle down to normal speed after a diving start and a 10% in- 
crease in velocity at the start of a one-mile course will give 
about 6% increase over the normal maximum speed. 


Limiting speed in a dive. The maximum speed in a dive 
with full throttle can be obtained from Figure 132 or Figure 133. 
If the propeller is idling, the thrust will be negative, and the 
value will depend on the pitch/diameter ratio. Data on pro- 
peller drag are given in Chapter VI. When propeller thrust is 
neglected, the solution gives the limiting speed at sea level 


vi =4/% (157) 


where K is the drag coefficient defined by D = KV’, and W is 
the gross weight. 
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While the time required to attain the theoretical limiting 
velocity is infinite, the initial approach is very rapid. If the 
fall be from rest, as from a stall or a loop, the relation between 
velocity and time is 


Vin Valea | (158) 


28 
where a = —.- 
T : 
The time required to attain any speed V is obtained by solv- 


ing equation (158) for é, and is 


t= x log. as Vv (159) 

The relation between velocity and altitude is 
V = Viva — e&*) (160) 
where b = re and h is the total altitude lost. If / is in feet 


and ¢ is in seconds, Vz must be in ft./sec. to determine a and 6. 
It is of interest to find how quickly an airplane picks up speed 
in a dive. Figure 135 gives the speed against altitude lost for 
limiting speeds of 200, 300, 400, and 500 ft./sec. The following 
table may be used in plotting similar curves for particular 
airplanes: 
a == 


— 0.2 0.4 0.5 0.6 0.7 0.8 0.9 
Va 


a > = 00064 .0027 .00448 . 00693 .O105 .O159 .0258 
ye 


It is possible to calculate the effect of density on diving speeds 
but the result is of little practical value. 

At some moderate altitude, an airplane can attain a speed in 
excess of the limiting velocity at that altitude if the dive is 
started at a very high altitude. If the dive is continued to 
lower altitudes, the speed will decrease slightly owing to the 
increasing air density. 
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Figure 135. Variation of Air ita with pee Velocity and Altitude 
Lost in a Vertical Dive 
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Take-off run. The distance run by an airplane in taking off 
is readily calculated if the following assumptions are made: 


Air forces vary as V?. Thatis, L = K,V?and D = Kp V’. 

Average coefficients may be used when the angle of attack 
is variable. 

Thrust increases uniformity with velocity. T = T,’—aV 
(a = thrust gradient, 7,’ = static thrust). 

Friction drag varies as load on wheels and skid with 
constant coefficient of traction wu. 


The derivation of the equations from these assumptions is 
straightforward but too lengthy to be given here. The velocity V 
(ft./sec.) at the time ¢ (sec.) from the start is 
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2T, (e% — 1) 
v= a eee (161) 


and the total distance run is 


OI, 


‘ Ww 
S= C, t+ a loge | (162) 


C,e— C, | 
2 Va? + 4KT, 


The values of the constants are 


T,=T,’ —pW.  T.’ = static thrust, Ibs. 


Lig So 

a = thrust gradient. = i 
Cy = 0+ Ve Pais, 
C, =a— Va? + 4kT, 
W = gross weight, lbs. 
K = Kp — pK 

ot! Wea 
b W a?+4kKT, 


The time to reach any velocity V may be obtained by solving 
equation (161) for ¢. 


I 27T,— GV 


The foregoing equations are for calm. When the take-off is 
into a wind velocity Vy, the time required is ¢ = ¢, — ¢,, where - 
t, is the time required to attain the speed V = Vy, and #, is 
the time of the getaway speed, both ina calm. The air distance 
traveled is S = S, — S,, where S, and S, are the distances cov- 
ered in times ¢, and #¢, respectively. The ground distance, or 
ground-run is S,= S — Vp-t. 

If V and S are calculated at various intervals and plotted 
against ¢, the times and distances may be read directly for any 
wind speed and getaway speed. A convenient method of pre- 
senting these data for any given airplane is in the form of a 
curve of ground-run against wind speed. The calculations for 
this curve may be tabulated as follows: 
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1. Wind velocity, Vy ft./sec. 
2. Time saved by wind, Af. Read from curve of V 
against ¢. 
3. Time to take off, ¢ = 4, — At. ¢, is time to take off in 
calm. 
4. Air distance at start of take-off, AS. Read from 
curve of S against ¢. 
5. Air distance during take-off, S’=S,— AS. S, is 
total distance in a calm. 
6. Wind effect, Sw = Vy: t. 
7. Ground-run, S = S’ — Sy. 
The static thrust is calculated by the empirical equation 
pb] BAHb.. 
T, = 6,000 18.7 ~ 9.5 5 (r.p.m.)-D (79) 


where # is the pitch and D is the diameter. The thrust at any 
speed V may be found from 
T = 315 T.Hp. 
V mi./hr. 
in which T.Hp. is the thrust power at the speed V. This can ° 
be obtained from Figure 96. 

Approximate take-off run. A very close approximation to 
the take-off run may be obtained from a simple formula. The 
distance S required to attain a speed Vs with an average accelera- 
tion a is 

S= = (164) 


The average acceleration force on an airplane in take-off is 
proportional to the thrust power, or F = Kyn-B.Hp. Hence 


sks (165) 
(a5) 


The average efficiency in take-off is not greatly affected by 
propeller characteristics so that 


S= KV 5? Es 


Kn -B.Hp. = 


oss 


(166) 


226 ENGINEERING AERODYNAMICS [Ch. 13 


Values of K have been determined from a number of ob- 
served and calculated take-off runs. If allowance be made for 
the known B.Hp. instead of the nominal B.Hp., and if values of 
S are taken from Figure to instead of Table I in N. A.C. A. 
Technical Report No. 249,3 the observed value of K varies from 
.o11g to .0132 with an average value of about .o125. Values of K 
from eight calculated take-off runs range between .o1ro to .o126 
with an average value of .o118. The following values may be | 
used: - 


PURSUE Misch creer aes a ee KG wons 
Observation 204 moa cea Ko = 30122 
Bight-Bombers ss. 5 oe eee K = .o120 
eave BOMDELSianwNa-ss see ee K = .o118 


If the engine is geared or if the propeller does not have the 
normal pitch and diameter for maximum efficiency at or just 
below high speed, it’ is better to write equation (166) in the 


form 


2.25 
S = .0215 Vs i (166a) 


Where Vs is the take-off speed in mi./hr., W the gross weight 
and T, the static thrust (See equation (79) ). Equation (166a) 
is more accurate than equation (166). 


General curve for take-off run. If the ground-run required 
to attain the getaway speed Vs is So in a calm, and Sj in a wind 
velocity Vw, then the ratio Sw is substantially the same at any 

; Vw, ea ; : 
given value of ae for all airplanes. Figure 136 is a plot of li. 
Ss ) 


against a as determined from six calculated take-off runs for ~ 
S 


widely varying types of airplanes. This curve may be used to 
construct a curve of ground-run against wind velocity for any 
airplane for which the ground-run in a calm, So, and the stalling 
speed, Vs, are known. If great accuracy is not required, S$, 
may be obtained from equation (166). 


3 Thomas Carroll, “A Comparison of the Take-off and Landing Characteristics of a 
Number of Service Airplanes’’ (1927). 
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Figure 136. General Curve of Effect of Wind on Take-off Run 


Landing run. It may be shown?’ that the landing run is 


given by 
S= ea) ae( (7) ) (167) 


L 


where Vz is the landing speed in feet per second, ? the reciprocal 


of the - at the landing angle, and v. the coefficient of friction. 


4H. Glauert, ‘The Landing Speed of Airplanes,’’ Br. A.C. A., R.& M. No. 666, 
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The value of p. is higher than that in the take-off owing to the 
high value of the tail skid drag. According to Glauert, the 
value of » for the skid alone is 0.50, and the average effective 
value of » for both skid and wheels is about 0.12. Obviously, 
this value is much affected by the shape of the skid and the load 
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Figure 137. Relation aoe ies Run, Landing Speed, Coefficient of 


Friction and 5 


that it carries. Since there is considerable variation in landing 
speed according to the method of landing, a similar variation in 
run is to be expected. 


equation (167). A study of this figure in connection with 
known landing runs leads to the conclusion that the average 
value of » is probably between .o6 and .10, and that the value 
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S 
of — Ve is normally between .oo8 and .or4o. The value of 2 is 
L 


considerably reduced by increasing the angle of attack of the 
wings when the skid is on the ground. 


Descent of a parachute. The terminal velocity of a para- 
chute, in ft./sec., is 


V= DOS = Me (168) 
Cp eS K&s 
Po 


where W is the gross weight in lbs., 9 the air density (lbs.-ft.-sec. 
units), and S the total area in sq. ft. Cp varies according to the 
shape of the parachute, between 0.70 and 0.80. The corre- 
sponding values of K are about .o017 and .oo1g. (See Chapter 
IV.) 

The usual method of determining K by timing the drop from 
a known altitude may lead to values as much as 10% too high 
unless correction is made for air density and its variation with 
altitude. A simple method of making the correction is available. 
For any increment less than 6,000 ft. in altitude, the density 


ratio - is linear with altitude 
eae (1 — ah) (169) 


The parachute may be assumed to descend with its instantaneous 
limiting velocity at all altitudes. The equation of motion is 


dV Per et Oy 2 
pa a wo ah) V (170) 


a is very small in comparison with other terms since it has a 


Be age ONE dV 
value of the order of .o1. Substituting Zi for V, neglecting ee 
and integrating, gives 


Lo V/ = . 2a — ah,)'5 — (1 — ah.) | (171) 


h, being the higher altitude. 
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The value of a must be determined by the densities at h, 
and h,, that is, from pressure and temperature readings. Since 


ae 0.3789 where P is the pressure in mm. of Hg and T 


is the absolute temperature in °C., the value of @ may be ob- 


tained from 
0.378 122. /2e 
ge sane E Z TF | ere? 


In the standard atmosphere below 6,000 ft., a@ = .000027. 


For a parachute having a limiting velocity / us = 20 ft./sec. 


in air of standard density, the time of fall from 4, = 2,000 ft., to 
h, = 0 ft., (sea level), is 98.6 sec. giving an average velocity of 
V = 20.28 ft./sec., which reduces K about 3%. The error is 
approximately proportional to the altitude. 

It is interesting to note that limiting velocity of a parachute 
may vary more than 5% with the normal sea-level temperature 
changes. 


The “flat spin.”” To the casual observer, the only difference 
between a conventional spin and the so-called ‘‘flat spin” is in 
the attitude of the airplane. In a normal spin, the longitudinal 
axis is inclined well downward and is nearer vertical than hori- 
zontal. In the flat spin the longitudinal axis may be almost 
horizontal and the motion almost pure yawing. Very few air- 
planes can be kept in a flat spin and once in a flat spin recovery 
is always difficult. The reasons for this are now well known as 
a result of independent and correlated investigations carried out 
by the Army, Navy, and National Advisory Committee for 
Aeronautics.’ These investigations included a large number of 
actual spins carried out by the Army and Navy. 


s Montgomery Knight, ‘‘ Wind Tunnel Tests on Autorotati in,’” 
N. A.C. A. Technical Report No. 273 (1927). ‘ me AsttorotetlOn 90. Sa ee ae 
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The flat spin is due chiefly to three major factors and perhaps 
to two or more minor factors. The major factors are wing 
arrangement, mass distribution, and center of gravity location. 
The minor factors are control size, control efficiency, etc. 

In order that a flat spin may be maintained, it is necessary 
that rotation at high angles of attack be a stable condition. 
This condition apparently exists in all biplanes having wings of 
equal chord without stagger. On the other hand, it is probable 
that no monoplane will autorotate at very large angles of attack. 

The second major condition is mass distribution. A solid 
body always tends to rotate about its principal axis having 
the greatest moment of inertia. An airplane which has a 
much larger moment of inertia about its normal axis than about 
its longitudinal and lateral axes is subject to large precessional 
moments in the combined roll and yaw of a spin. The major 
precessional moment is a stalling moment which holds up the 
nose of the airplane. If the c.c. is well aft, say beyond 35% of 
the mean chord, the precessional moment may be greater than 
the elevator moment, and recovery from a spin is very difficult. 

‘The obvious remedies in biplane design are to use stagger 
and overhang if possible, preferably with some taper in the 
wings, and to keep the wing weights low in order to reduce the 
moment of inertia in yaw. Large tail surfaces are of consider- 
able assistance in preventing a dangerous spin or in assuring 
quick recovery from any spin. A majority of the airplanes 
which have been reported to recover slowly from spins are de- 
ficient in horizontal tail surface area. 

The results of the British investigation of spinning are 
given in R. & M. No. 1oo1, “The Spinning of Aeroplanes,” by 
W. B. Gates and L. W. Bryant (1927). This is a very extensive 
and valuable report, covering all phases of the spin. 


CHAPTER XIV 
PERFORMANCE ESTIMATION 


The aeronautical engineer frequently finds it desirable to 
estimate the change in performance due to slight changes in de- 
sign characteristics, not great enough to justify a complete 
performance calculation. Methods for the rapid estimating of 
performance are also of value for other purposes where an 
approximate solution will serve. An experienced estimator can 
usually obtain results which agree very well with accurate 
calculations or flight test data. Numerous methods for estimat- 
ing performance have been published from time to time, but 
none can be considered as a substitute for the complete cal- 
culation, if accuracy is essential. The methods described in 
this chapter have been thoroughly tried by the author and 
are believed to be sufficiently reliable for most purposes. 


Stalling speed. The stalling speed, in miles per hour at 
sea level is calculated from the equation 


/ 


where * is the wing loading in lbs./sq. ft. If the wing section 


is normal in thickness and camber, the maximum value of C, 
or K, may be read directly from wind tunnel test data, but the 
model chord should be greater than 4 in. and the wind speed 
greater than 50 ft./sec. in the test. If the section is thin, or 
“double cambered,” the lift coefficient may increase consid- 
erably with Reynolds Number. If the section is thick or 
deeply cambered, the lift may decrease with Reynolds Number. 
(See Chapter I for additional information.) 
232 
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I Por derivation of this and other performance formulas see W. S. Diehl, ‘Reliable 
i i i .C. A. Technical Report No. 


ormula 
173 (1923). 
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in which Vy is the maximum speed, Vs the stalling speed, Z 
the maximum overall lift/drag ratio, y, the maximum pro- 
peller efficiency, ee the power loading in Ibs. per B.Hp., and K a 


constant which depends on the effective aspect ratio and parasite 
resistance as shown on Figure 119, Chapter IX. The average 


x 
value of K is about 10.2, and the average value of K @ . nm)? is 


between 19.0 and 20.0. Figure 138 is a plot of equation (114) 
enabling maximum speeds to be read directly. 


Initial rate of climb. The initial rate of climb in ft./min. 
may be calculated from the formula 


dh _ = Ks Ym Ve 
(4 5 = Ae a > | (173) 


in which & is the power loading in lbs./B.Hp., 7m the maximum 
propeller efficiency, Vc the climbing speed in mi./hr., 5, the 
maximum overall ratio of lift to drag, and K, the ratio of power — 
available at climbing speed to power available at maximum 


speed. Analysis of power curves for a number of airplanes 
shows that an average value of K, is 


Vu -0.27 


K. = (7) Ba aise) 


K, is plotted against a in Figure 139. This average value of K, 
s 
is based on a climbing speed Vc = Vs + 3 (Vu — Vs) = 2VstVu 
3 


Actually, the climbing speed varies considerably with aspect 
ratio and parasite as shown by Figure 120, but the effect on the 
rate of climb is not great. However, if desired, a more exact 
value of Vc may be used and the value of K, determined from 
Ke 


Figure 96, using the ratio 
Vu 
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Figure 139. Thrust Power Coefficient K2, in Equation (173). Values of Kz 


are Above and Values of o are Below the Base Lines 


Absolute ceiling. When plotted against speed, the curves of 
T.Hp., and T.Hp., are tangent at the absolute ceiling. The 
point of tangency is very near the minimum power required, and 


normally at a speed approximately ae Vup, where Vyp 


is the speed for minimum power at sea level. These relations 
enable the absolute ceiling to be given in terms of the ratio 
Hn 
AP. 
and Hp. is the power available at sea level and at the air speed 
of the tangency. It may be shown that 


where Hp.,, is the minimum power required at sea level, 


L 
AD. ae Ye Tm ° Hp. a K, ae (Z 


Hh.r0 ke aa K,vs(#) (175) 


373 (5) 
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Analysis of power curves shows that K, 


Vu 
V. 


of 


stituting the value 


from equation (114) gives: 
s 


ny Lee 


Hp 
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substituting this into the equation (175) gives 


ice 16; 0.80 
H: LLP. a apy 
ee RC Reb 
( 5) Ym : aup, 


Neither K nor K; vary enough to affect the value of the expres- 
sion in the brackets, so that the average values of K? = (10.2)? = 
5.7 and (K,)* = (1.10)? = 1.13 may be used. The product 
K? (K,,)* then has the average value 6.40. 

Equation (176) is most conveniently used when plotted as 
in Figure 140. To use this plot, first calculate the value of the 


“Se eee yee 
goth 


l2 | 16 20 24 28 32 36 40 
Figure 141. Absolute Ceiling Constant Kw in Equation (177) 
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denominator | S40 Vs: re | then pass vertically along this 


value to the proper - and then horizontally to the absolute 
ceiling. 

A slightly more reliable method of estimating absolute ceiling 
is based on the equation 


adh ae el wiry | (177) 
ke, (- ap 


where uf is the wing loading (lbs./sq. ft.), 7 the power loading 


(lbs./B.Hp.), 1 the maximum propeller efficiency, and Ky a 
factor which varies with aspect ratio and parasite as shown on 
Figure 141. This formula has given very satisfactory results. 


Service ceiling. If the rate of climb decreased linearly with 
altitude, the service ceiling (at which the rate of climb is too 
ft./min.) would be given by 


Cy — 100 
hs = H (ee) (99) 


0 


Owing to the curvature in the rate of climb curve, the actual 
service ceiling will be slightly below the value given by equation 
(99). Figure 142 is a plot of Avs against absolute ceiling Z, 
Ahs being the difference between the straight line value and 


the exact value of fs. Ahs is closely approximated by the 
relation 


HON 
A hs = 40 (—*) (178) 


consequently the exact service ceiling is given by 


wy TI) ( H. 3 
hs =H ( Gs mig ad) (179) 


fo) 8000 16000 24000 32000 


Figure 142. Difference Between Service Ceilings for Linear and for Non- 
Linear Rate-of-Climb Curves 


Time of climb. With a linear relation between rate of climb 
and altitude, the altitude / climbed in T minutes is 


per sere oy (101) 


where ZH is the absolute ceiling (ft.), and C. the initial rate of 
climb (ft./min.). A simple relation may be found between 
the altitude given by equation (ror), and the actual altitude 
for a non-linear rate of climb. 


When the ratio of < is plotted against - as in Figure 143, 


the various points (representing about 20 airplanes) lie very 
near a single curve, which indicates that all climb curves may 
be similar. Figure 143 may be used for drawing up estimated 
rate-of-climb curves when the initial rate of climb C. and the 
absolute ceiling H are known. It may also be used to determine 
some general relations. 
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Plotting from Figure 143 against i and integrating the 
area under the curve gives factors proportional to the times 


required to climb to any altitude #. A similar plot of a from 
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Figure 143. General Rate-of-Climb Curve 


equation (101) gives similar factors for linear rate of climb. 


The ratio of these factors gives the ratio between the times Te 
S 
to climb to a given percentage of the absolute ceiling, T¢ being 


the time for non-linear and 7’; the time for linear rate of climb 


Ch, 14] PERFORMANCE ESTIMATION 241 


1.00r 


-80 


© 
60] |< 
l}_s 
Oi 

=) 
Lobe 
elu 
aie 
4o| “3 
JO 
(Vp) 
(ea) 
<x 


PAH 
4 ; 
See eee 


1.00 1.04 1.08 112 bI6 1.20 


Figure 144. Relation Between the Times to Climb toa Given Altitude with 
Linear and Non-Linear Rate-of-Climb Curves 


against altitude. a: is plotted against Ag Figure 144. This 
s 
figure may be used to determine the value of Tc from a known 


or easily determined value of Ts. 
Ts Tc 


) () 


data in Table 20. This figure may be used to calculate either 
the time required to climb to a given altitude or the altitude 
climbed in a given time, for either a linear or non-linear rate of 
climb. For example, assume H = 20,000 ft., Co = 1,600 ft./min., 


against us using the 


and 


Figure 145 is a plot of 
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then a= = 12.5 min. The times to climb to 10,000 ft. ys 0.5, 


are Ts= 06) == 8.62 min., and Te= 015 = 9.37 min. 


The altitudes climbed in 10 minutes are found from Ts = 


(<) 


ae tO =0. 80, giving # = 0.55 from which hs = 11,000 ft., 


12.5 


and e = 0.52 from which hc = 10,400. 


TABLE 20. TIME OF CLIMB FACTORS 


em cela ae) * joe Fe 

CG; c. 3 
Ce) fo) fe) I.000 
.10 . 1053 1068 1.014 
.20 2232 2300 1.029 
.25 .2875 2982 1.037 
-30 . 3563 3723 1.045 
-35 -4305 4540 1.054 
-40 .511 544 1.063 
-45 .598 642 1.072 
.50 . 693 .749 1.081 
.55 .797 .870 1.091 
.60 . 916 1.009 1.102 
.65 1.051 ign aie TLL 
.70 1.203 1.352 1.124 
.75 1. 386 1.574 1.136 
. 80 1.609 1.849 1.148 
.90 2.303 2.703 Lanz 3 


The climb in 10 minutes is used extensively as a measure 
of performance. ‘Three charts based on the linear relation 
have been prepared for the solution of various problems involvy- 


ing this factor. Figure 146 is a plot of Are against Ge from which — 


h,, may be accurately determined. This value may then be 
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Figure 145. Time to Climb to a Given Altitude. Ts for Linear, Tc for 
Non-Linear Rate of Climb Curves 
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corrected for the non-linear relation, if desired, using one of the 
preceding figures. Figure 147 is a plot of /,, against H and C,, 
_ from which /,, may be read directly. Figure 148 is a similar 
plot with the service ceiling 4s instead of the absolute ceiling H. 
The latter figure may be used to determine initial rate of climb 
when only h,. and hs are given. 

The climb in 10 minutes is easily obtained from Figure 145, | 
as illustrated by the example on page 242. In general, Figure 
145 will be found most convenient for this purpose, but there 
are many problems for which Figures 146-148 are especially 
adapted. 
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Figure 146. Climb in 10 Minutes, from Absolute Ceiling and Initial Rate of 


Climb (linear). Values of tee are Above and Values of 2 are Below Base Lines 
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Figure 147. Climb in 10 Minutes from Absolute Ceiling and Initial Rate of 
Climb (Linear) 


AWN] | 


Va 


Ch. 14] PERFORMANCE ESTIMATION 245 


Range and endurance. The endurance at maximum speed 
is obtained by dividing the full load by the full throttle fuel 
consumption in lbs./hr., using the known specific fuel consump- 
tion and actual B.Hp. if known, otherwise the rated, or nominal 
B.Hp. If the specific fuel consumption is not known, it may 


10000 
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Figure 148. Climb in 10 Minutes from Service Ceiling and Initial Rate of 
Climb (Linear) 


be estimated: for modern engines it varies from 0.46 to 0.55 
Ibs./B.Hp./hr., depending largely on the treatment the engine 
receives. An average value of 0.50 will not be far off. The 
range at maximum speed is obviously the product of the corre- 
sponding endurance and speed. 

After the endurance and range at maximum speed have 
been determined, the endurance and range at any other speed 
are readily found by the factor method’. This method is based 
on the factors 


2 W.S. Diehl, ‘‘Three Methods of Calculating Range and Endurance of Airplanes,” 
N. A.C. A. Technical Report No. 234 (1925). 
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Figure 149. Range Factors in Terms of Speed and Speed Range 
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Figure 150. 


‘ MAXIMUM SPEED 
! MENACE ada STALLING SPEED 
46 18 20 


22 24 26 28 3.0 


Figure 151. Maximum Range Factors 


Ch. 14] PERFORMANCE ESTIMATION 249 


Endurance at speed V 
Endurance at max. speed 


Endurance Factor Fr = 


and 
Range at speed V 
Range at max. speed 


Range Factor Fr = 


ENOURANCE AT MAXIMUM SPEED 
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Figure 152. Endurance Factors at Speed for Maximum Range 


Fr and Fz vary with the ratios me and re, as shown on Figures 
Ss M 


149 and 150. These figures are based on the assumed use of 
mixture control, the effect of which is shown elsewhere. 
In Figure 151 the value of Fg for maximum range is plotted 
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Vu 
Vs 
sponding values of Fz are given. These figures enable the maxi- 


against the speed range ratio —, and in Figure 1 52 the corre- 


MAXIMUM 


ENDURANCE 


RATIO 


Figure 153. Maximum Endurance Factors 


mum range and corresponding endurance to be calculated 
rapidly. 
The maximum endurance is given by the plot of maximum 


Fx against i in Figure 153. 
Ss 
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Example of performance estimation. Consider the airplane 
for which performance was calculated in Chapter VII. The 
characteristics were as follows: 


W = 2,100 lbs. S = 290 sq. ft. B.Hp. = 195 


Whe 2 YE ean 6 = 0.78 
time ak. B.Hp. — 7 nm = 0.7 


(=) = 8.2 Vs = 52.6 mi./hr. 
max 


Effective aspect ratio = 4.8 


Maximum Speed: Use equation (114). From Figure 1109, 
K = Io.Io. 
Vu _ 10.10 (8.2 X 0.78)3 
Vou 8 215256:X 107 76)8 
Vu = 2.24 X 52.6 = 118 mi./hr. (116) 


= 2.24 


Initial rate of Climb: Use equation (173). TM 2.24. 
5 
From Figure 139, K, = 0.805. 


x 0.805 X 0.78 _— (2 X 52.6 + 118) 
Fo 33,000 | 10.76 (FX375 268.2) 


= 1,130 ft./min. (1,100) 


Absolute Ceiling: Use equation (177). From Figure 141, 


Ky = 23,600 for a = 8.2, aspect ratio = 4.8. 


10,000 
1.78)" 
faa (2 78 


H = 23,600 log, I | = 20,200 ft. (19,600) 


Service Ceiling: Use equation (179). 


1,130 — | (Oe 
1,130 10,000 


3 
hs = 20,200 ( ) = 18,070 ft. (17,500) 


Climb in 10 Minutes: 


C 1,130 10 

Bey — y =, — . 6 

i, 20,300 G50 i tae 
Co 
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: h 
From Figure 145, ae 0.41. 


Ivo = O.41 X 20,200 ='8,280 It, (8.000) 


The values obtained in the exact calculation are given in 
parentheses alongside the estimated values for comparison. 


& 


CHAPTER XV 
SEAPLANE FLOATS 


The naval architect describes seaplane floats with a number 
of technical words and phrases peculiar to his profession. For 
the benefit of the student and engineer who is unfamiliar with 
these terms, a short list of definitions has been prepared, limited 
to the most frequently used words and phrases. 


Afterbody. That part of a float between the main step and the stern. 
Bottom. The area included between chines and keel from bow to stern. 
Bow. The extreme forward point, or portion of a float. 

Buoyancy. The displacement (in lbs. of sea water) to a given water line. 

Buoyancy, center of. The center of gravity of the displaced volume of 
water. 

Buoyancy, excess. The difference between the total or submerged and 
normal or load water line displacements. Usually expressed as a 
percentage of the normal displacement. 

Chine. The line of intersection of the bottom with the sides or deck. 

Deadrise. The angle which each side of the bottom makes with the 
hopizontal, as measured in a transverse plane. 

Deck. The upper surface between the sides. If the sections are rounded 
without flat or vertical portion, then all of the upper surface between 
the chines is called the deck. 

Deck line. The upper boundary of the float in a side elevation. 

Displacement. The weight of the sea water displaced to a given water 
line, or simply the load carried by a float under given conditions. 

Displacement submerged. The weight of the sea water displaced when the 
float is completely submerged. 

Draft. Usually refers to the maximum depth below water surface of any 
part of the float under given conditions. 

Forebody. That part of the float between the main step and the bow. 

Keel. The main longitudinal and continuous strength member located 
along the bottom and on the center line. 
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Keel, false. A protective member attached along the bottom center line 
on the outside, to prevent damage from handling or grounding. 
Metacenter. A point through which the resultant vertical buoyant force 
passes for all small displacements from the position of equilibrium. 

Metacentric height. The distance from the center of gravity to the 
metacenter. 

Porpoising. Any pronounced pitching oscillation in a moving float. 

Speed, getaway. The speed at which the entire weight of the seaplane is 
carried by the wings. 

Speed, hump. The speed or speeds at which the water resistance is a 
maximum. 

Sponsons. Lateral projections added to the sides of a float or hull to 
increase planing area or transverse stability. 

Spray strips. Thin longitudinal strips of triangular cross-section attached 
to the bottom along the chine to ‘“‘beat down”’ the spray. 

Squatting. A pronounced tendency to trim by the stern. 

Step. A line of discontinuity in a surface. In its usual form, a sudden 
change in transverse sections. 

Stern. The extreme rear point, or portion of a float. 

Trim. The angle of pitch, usually measured between the deck line and 
the water line. 

Trim by bow. An angle of trim produced by depressing the bow and 
raising the stern and measured in the same manner as trim by stern. 

Trim by stern. An angle of trim produced by raising the bow and depress- 
ing the stern and measured from a level position of some reference 
line, usually the deck line. 


Model tests—corresponding speeds. The only way to deter- 
mine the water resistance of a float or hull is by actual test, 
either full scale or on a model. Since most of the resistance 
is ‘‘wave-making,”’ comparisons must be made on the basis of 
+ or a That is, the “‘corresponding speed”’ is proportional 
to the square root of the length. 

Froudes’ Law of Comparison states that at corresponding 
speeds the full-scale wave-making resistance is equal to the 
model wave-making resistance multiplied by the cube of the 
linear scale ratio of full size to model. 


Letting V, 
Ry 


model speed and V = full-scale speed 
model resistance and R = full-scale resistance 
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and L = linear scale ratio £ ull size 
model 


the following relations hold 


Corresponding speeds V = V,VL 

Resistance at corresponding speeds R = R, L3 

If a model is 4 full scale, i.e., if L = 9, then the corresponding 
speeds on the model are one-third of the full-scale speeds, and 
at corresponding speeds, the full-scale resistance is 93 = 729 
times the model resistance. Since both the displacement A and 
the resistance R vary as the cube of the scale ratio, it follows 


that the ratio S is the same for model and full-scale at corre- 


sponding speeds. Model test data on seaplane floats and hulls 


are usually given in the form of curves of o against the ratio ae 
G 


where Vg is the “getaway” speed at which the float leaves 
the water. 

A given model test representing a definite full-scale seaplane, 
also represents the infinite number of full-scale getaway speeds 
and linear scale ratios giving the same value of model getaway 
speed. For example, the corresponding getaway speed for a 
linear scale ratio of 12 and a full-scale getaway speed of 52 
knots, is 15 knots. A model tested under these conditions 
would represent any geometrically similar full-scale getaway but 


to a different scale ratio, which is obviously ZL = (=) . That 


is, it would represent a 45 knot getaway to the scale ratio 9, 
or a 60 knot getaway to the scale ratio 16, etc. 

In addition to resistance at various trims, model tests should 
include, measurements of righting moments at rest, trimming 
moments underway, planing and spray characteristics. The 
righting moments at rest give the metacentric heights, while 
the trimming moments under way bring out any undesirable 
or uncontrollable tendency to dive or ‘‘squat.”’ The planing 
characteristics are rather general, and include items such as 
tendency to “porpoise” or “sticking at high speed.” The 
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spray characteristics are observed ordinarily to avoid an arrange- 
ment that allows spray to enter the propeller disc, or perhaps 
more exactly in some cases, to reduce this spray to a mini- 
mum. 


Model basin methods. Since the model basin equipment 
and general methods of testing seaplane floats at the Washington 
Navy Yard have been fully described elsewhere,* the following 
résumé has been prepared chiefly for the benefit of aeronautical 
engineers to whom the original papers are not readily avail- 
able. 

The model basin is simply a tank about 42 ft. in width by 
470 ft. in length and 15 ft. in depth at the center. A “carriage” 
extending across the basin, is driven on carefully aligned tracks 
by four electric motors using the Ward-Leonard speed control. 
The carriage is fitted with all of the speed control and dyna- 
mometer gear; it weighs approximately 75,000 Ibs. Its maxi- 
mum speed is approximately 18 knots, but testing is very rarely 
ever done at speeds greater than 15 knots, chiefly on account of 
the short run available and the arresting difficulties at the 
higher speeds. 

In testing seaplane models, certain assumptions are necessary 
in order to arrive at practical methods. First, the model is 
attached to the dynamometer by means of an arm about 8 ft. 
in length. This arm is hinged at the dynamometer attachment. 
and guided so that it is free to move in a vertical plane. The 
model attachment to the dynamometer arm is such as to allow 
the model to pitch freely about the full-scale center of gravity 
on the ‘‘free-to-trim” runs. On the ‘fixed trim” runs the 
attitude of the model relative to the arm is fixed. It is assumed 
that the angular changes due to vertical motion of the model 
are negligible and this is easily verified. Second, the model 
displacement at rest is adjusted by means of weights carried on 
an extra flexible wire cable passing over two pulleys to the 


tD. W. Taylor, ‘The United States Experimental Model Basin,” Tra i 
the Society of Naval Architects and Marine Engineers, 1900; H. C. Richawiane Tee 


ayeed of pivisg Boat Development,” Journal of the American Society of Naval Engineers, 
ay 1926. 
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model. These weights also counterbalance the dynamometer 
arm. Obviously, all of these weights in the system seriously 
affect any comparison of model oscillations with full scale. 
Third, it is assumed that the wing lift varies directly as the 
speed so as to lift the model entirely clear of the water at the 
“corresponding” model getaway speed. This lift is automatically 
applied directly to the counterbalance cable by a small vane 
which is calibrated by trial runs to determine the correct setting. 
When this setting is once determined, it is not changed during 
the test although the wing lift is proportional to the angle of 
attack as well as to the square of the speed. It may be shown, 
however, that the error introduced by this approximation is 
negligible, and some allowance is made by taking the getaway 
speed from 5% to 10% higher than the calculated actual 
value. 

That these assumptions do not seriously affect the validity 
of test data has been amply proved by a number of comparisons 
between model and full scale. The best method yet found for 
an approximate check is in the comparison of predicted and 
actual maximum loa ‘s that can be taken off by a given seaplane, 
in a calm or in a wind of known velocity. The agreement in 
every case investigated has been uniformly close, indicating 
less than 5% difference. 


Typical model basin data. Table 21 contains the average 
values of S against the ratio a for single floats, twin floats, 
and flying boat hulls as given by Captain H. C. Richardson? in 
his historical paper ‘‘ Naval Development of Floats for Aircraft.” 
These average values may be used for any general purpose, 
but too much dependence should not be placed on their use in 
a performance calculation that indicates a low margin of thrust, 
because the possible deviation from the average is consider- 
able. 


2 Transactions of the Society of Naval Architects and Marine Engineers, Vol. 34, 1926. 
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DISPLACEMENT 
TABLE 21. AVERAGE RATIO ———————— FROM MopeEt Basin TESTS 
RESISTANCE 


Average Ratio of z 


Speed 
Jo 
Getaway 
Speed Single Twin Boats 
Floats Floats 

20 8.05 8.10 10.30 
25 6.50 6.60 7.40 
30 5.10 5.30 6.10 
35 4.30 4.50 5-30 
40 4.10 4.10 4.90 
45 4.05 4.10 4.80 
50 4.30 4.25 5-95 
55 4.70 4.65 5-25 
60 5.15 5.15 5.50 
65 5-55 5-45 5-50 
7° 5-75 5-50 O50 
75 5-55 5-30 5-00 
80 5.00 4.70 4.50 
85 4-15 3.90 3.90 
90 3.20 3.05 3.05 
95 2.00 2.30 1.80 
100 fe) fo) 


Take-off in acalm. The resistance due to the floats or hull 
in a take-off is readily obtained by the following steps: 


1. Calculate getaway speed Vg 
2. Assume a series of speeds V 


3. Find ratio (4) for each speed 


Ze 
A V : 
4. Read (= ) at each{ —) from model basin curves 
sa Ve 
5. Calculate A at eacn speed: A = (w I— 7) 
A 


6. R= 


The total resistance is found by adding the air resistance of 
the airplane to each value of R. The air resistance is found by 


Ch. 15] SEAPLANE FLOATS 259 
assuming a constant angle of attack in take-off, say 8°, and 


calculating (or simply estimating) a value of 4 corresponding to 


this angle. The air resistance is D= Sed where Ly is the 


r 2 D 
wing lift= W oS ) . If the total resistance and maximum 


propeller thrust are plotted against speed on the same diagram, 
the margin of thrust available for acceleration is easily obtained 
and from this one can obtain the probable take-off performance. 


Take-off in a wind. The take-off may be made either into 
or down wind, but the effect of the wind is so great that few 
seaplanes can take off down wind when the wind velocity V, is 
large. If there is any wind at all, the air speeds and water speeds 
are no longer equal, and the wing lift no longer follows the 
curve assumed in making the Model Basin test. For take-off 
into wind, the wing lifts are greater and for take-off down wind 
they are less than the assumed values. This is simply equivalent 
to changing the load on the float at a given speed. Let us 
consider the effect of changes in load and changes in getaway 
speed on the values of 7 Figure 154 is a plot of a against =. 
for an NC type of boat hull with three loads and with constant 
getaway speed, as obtained in tests at the Washington Navy 


Yard. The values of S fall upon a single curve, indicating that 
with constant getaway speed the effect of a 25% change in load 
carried by the float is negligible. Figure 155 is a plot of a 
against for an Fs5 type of hull with three loads and with 


getaway speed varied according to the load, also from Wash- 
ington Navy Yard test data. As might be expected, there is a 


slight decrease in the values of a as the load and getaway 


H A= 440LBS,\:1665 K 
© A.= 5.50 LBS, Vs:I665 K 
x A .= 6-60 LBS,Vel6.65 K 


DISPLACEMENT — 
RESISTANCE 


fe) 20 40 60 80 100 
Figure 154. Effect of Variable Load on = with Constant Getaway Speed 


x A 633 LBS., Vo=13.86 K 
© A 7.56LBS., Ve=I5.13 K 
O A 889LBS, Ve=I6.41 K. 


i 

Z/ ul 
ui] O 
lz 
Ole 
<u 
Qi 
Y 
ray 


Le) 20 40 60 80 100 
A 


Figure 155. Effect of Variable Load and Getaway Speed on R 
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speed are increased, but the changes are negligible for engineer- 


ing purposes. Values of S may therefore be read directly from 


‘ paar fae : 
the standard curve, if the ratio y,38 determined from the water 
G 


speeds. Resistances are found from these values of a and the 


corresponding displacements. The procedure will be outlined 
in detail, for clarity. 


Take-Off into Wind: Denoting wind velocity by V, and 
the getaway air speed by Vg, the getaway water speed is Vew = 


Vco—Vo. If the water speed is Vy, then the values of are 


determined by the ratios of 7, and resistances are obtained 
as follows: ake 


. Calculate getaway air speed, Vg 

. Assume value of wind velocity Vo 

. Calculate Vew = Ve — Vo 

. Assume a series of water speeds, Vy 


mn PWN A 


. Find ratios Vw 
Vew 


6. Read = from standard curve, assuming i equiva- 


R 


V 
lent to Ve 
. Obtain air speed, V = Vw+ Vo 


. Find wing lift at each water speed, assuming Ly = 


Wye) 


9. Load carried by float A = W— Ly 


on 


10. Water resistance bree (2) 
R) 


11. Air resistance D= (4) 
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Take-Off Down Wind: In this case the water speeds are 
greater than the air speeds and it is necessary to consider Vo 
negative. 

The method just outlined must be considered as an approxi- 
mation for use in the absence of more exact data. It gives 
results which agree quite closely with observed full-scale maxi- 
mum loads which can be taken off in a wind of known velocity. 

Thrust may be assumed to vary linearly with speed, making 
it necessary to calculate the thrust at two speeds only, one of 
which may be the static thrust. (See Chapter VI.) 


Metacentric height. Metacentric height may be defined by 
considering a floating prism having its c.c. at a point G and its 
center of buoyancy at a point B. A line perpendicular to the 
water surface and passing through B will also pass through G. If 
the prism be inclined through a small angle 6 while retaining 
the same volume of displacement, the center of buoyancy will 
shift to a point B’. A vertical through B’ will intersect the 
original vertical BG at a point M, which is called the metacenter. 
The distance GM is called the metacentric height. 

It is easily shown that the metacentric height is a measure 
of static stability. Considering a slight inclination @ and 
taking moments about the original center of buoyancy B, it is 
seen that the disturbing moment is W -BG-sin® and the 
righting moment is W-BM-sin 9. The total resultant 
moment is 


M = W- (BM — BG) -sin80 


= W-GM-sin0@ 

and the slope of the resultant moment curve is 

d 

oe = W-GM-cos@ 
from which 

d I F 

GM = “add : se, (for 6 in radians) 

c Gye uae 


a a (for 6 in degrees) (180) 
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This relation is used to determine metacentric heights from 
inclination tests on models. 

The metacenter may be found by the use of the formula 
BM = I/V, where I is the moment of inertia of the waterplane 
about its center line and V is the total volume of displacement 
(see any book on Naval Architecture). J varies as L4, and V 
varies as L3, so that BM and the metacentric height vary directly 
as the length, or scale ratio. 

Approximate metacentric heights for seaplane floats and 
hulls may be obtained from the empirical formulas given below. 


Metacentric height required. Analysis of the performance 
of a great number of seaplanes indicates that satisfactory static 
stability is obtained when 


Transverse GM = longitudinal GM = 1.4 (A)! (181) 


where A is the gross weight of the seaplane. 


Transverse metacentric height of twin floats. It has been 
shown’ that with the design proportions in common use, the 
transverse metacentric height for twin floats is given closely 
by the empirical formula 
KL s*B 

a 
where L is the length and B the beam of each float in ft., s the 
spacing on center lines in ft., A the gross weight of the seaplane, 
and K, a constant varying from 17.7 to 20.8 with an average 
value 19.5. Equation (182) may be used to determine the spacing 
necessary for stability by substituting the average value of 
GM from equation (181). That is 


i era: LS a0) 0.28 Ai 18 
aes ¥ | /LB ( 3) 


GM = (182) 


Longitudinal metacentric height. It has been shown in 
N. A.C. A. Technical Note No. 183 that the longitudinal meta- 


3 W. S. Diehl, ‘Static Stability of Seaplane Floats and Hulls,”” N. A. C. A. Technical 
Note No. 183 (1924). 
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centric height for either single or twin floats is given with suff- 
cient accuracy by the empirical equation 


y 
GM = EE: "ee (184) 


where is the number of floats (i.e., one or two), B the beam 
of each float in ft., ZL the length in ft., A the gross weight of 
the seaplane, and K, a constant normally varying between 1.90 
and 2.40 with an average value of 2.10. 

Equation (184) may be used to determine the minimum 
length of a seaplane float for longitudinal stability by substituting 
the value of GM from equation (181) 


LAAT) O06 Le 


— = 
L 2.10nB nB (185) 


Transverse stability—single floats and flying boats. Single 
float seaplanes and flying boats require the use of auxiliary 
floats on the wing tips in order to secure transverse stability. 
The usual method is to install the tip floats so that at normal 
trim and draft they clear the water from 3 in. to 8 in. according 
to the span, and make the volume large enough to give a righting 
moment which is greater than twice the upsetting moment when 
one float is just submerged. 

The equation of moments about the c.B. is obviously 


A:lcos8 = C- Wh- sin 98 
from which, the wing tip float displacement is 
A; = C- WI- tan 0 (186) 
where W is the gross weight, / the height of the c.c. above the 
center of buoyancy (% = BG), J the distance from the c.B. of 
the tip float to the c.B. of the main float or hull, ® the angle 


of inclination required to submerge the tip float, and C a con- 
stant greater than 2.0 and preferably between 2.5 and 3.5. 


Section coefficients. The ratio of maximum cross-section to 
the product of beam by depth varies according to type as follows: 
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A 
BD 
Average single float.. SER Snitaro se amee BTS 
Average twin float......... Sey DOI OSE ERR Rana he 
Average boat hull, no ) sponsons.. Pe Camera 70 
Average boat hull, with sponsons............... 55 


These values are useful in making performance estimates. 
They should be regarded simply as averages, for this or similar 
purposes. 


Prismatic coefficients. The ratio of total (submerged) vol- 
ume to the product of length ZL, beam B, and depth D is the 
prismatic coefficient. This coefficient is variable with the type 
of float or hull, and has the following average values: 


Prismatic 
Type Coefficient, Cy 
Single oats tae accikenie meals ce ser vee eae .70 
wan Hoa tsiti «seminar iiaete oaeisiaid aome strc ee 72 
Boatenullsy MO! SPONSODSn qa. sirohinctacie aera cia ase .50 
Poatenullstwith SPOMSONS messes eae ee nes 5B 
Wine tip oatsees cries oe ara yt danes eeeeies 55 


The submerged displacement in sea water weighing 64 
Ibs./cu. ft. is given by 
As = 64 Cy LBD (187) 
L, B, and D being in ft. 


APPENDIX 
STANDARD ATMOSPHERE 


The “standard atmosphere” is an arbitrary variation of 
temperature, pressure, and density, with altitude, which is used 
for numerous aeronautical purposes, but chiefly as a basis for 
comparing performance. The standard atmosphere recom- 
mended by the National Advisory Committee for Aeronautics’ 
and adopted in 1925 by all interested Government Departments 
for official use in the United States is based on the following 
assumptions: 


Ground temperature, t, = 15°C. = 59° F. 

Isothermal temperature, f = —55°C. = —67°F. 
Temperature gradient, a = .0065° C./m. = .003566° F./ft. 
The air is a dry, perfect gas. 


The resulting equations are: 


T =T, — ah (188) 
p (2) T 
miReeh or =|). 18 
p Ses Bey aan (189) 
(Fim) 
h= — (=) log. (— 190 
| oh OREN (190) 
Tm = Harmonic mean temperature 
ey ah 
meats) (191) 
pee T, — ah 


Below the isothermal level (10,769 m. or 35,332 ft.) the 
following relations exist 


ISee ‘Standard Atmosphere—Tables and Data,” N.A.C.A. Technical Report No. 
218 (1925). 
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a i (2) Gre eS, 
C)- = es 
-Q"-O" 
= = (1 - = h) (195) 
) = (1 = ny (196) 
ae = (1 See h) (197) 
(= (4) (198) 


In the foregoing equations, the subscript o refers to the 
standard conditions at sea level. / is the altitude and R the 
gas constant for air. 


Standard atmospheric relations used in performance reduc- 
tion. The relation between pressure, temperature, and density 
ratio, at any altitude is 


e pb mm. Hg. 


ao 0.3789 (f°. + 273) (199) 
- p in. Hg. 
i 9.624 [t SO we 273°] (200) 

p in. Hg. 


= 17.32 (201) 


[¢°F. + 459.4°] 
The relation between an increment of pressure Ap= 


(p, — p2), the average density em between p, and p,, and the 
increment of altitude Az is, in metric units: 


— 13.59 Ap mm. Hg. 


Ah= (202) 
£ Om 
and in English units 
Akh= — 70.67 Ap in. Hg. (203) 


20m 
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Figures 156, 157, and 158 are plots of a - and p against 


altitude in standard atmosphere. Density ratios may be cal- 
culated by equation (199), and the corresponding density 
altitudes found from Figure 156. Pressure altitudes may be 
read directly from Figure 158. 


Figure 159 is a plot of Pe against altitude. 
e 


Velocity of sound in air. The velocity of sound in any gas is 


ae eo 


where # is the pressure, e the density, and K the ratio of specific 
heat at constant pressure to that at constant volume. The 
equation may be written in the form 


a= K be (2) (2 = - f/K@xt (205) 


for air K= 1.41, p= ws) ) ft., and eo = .002378 lbs. 
sec.?/ft.4 Substituting these values gives 


2116 T T 
a kite cag Cee Uae ae 206 
a TeAtex "002378 ( ) 1,120 388 (206) 


T being in °C. absolute. The variation with temperature is as 
follows: 


etc. —20 oO 10 15 20 30 
eat sec, -1,050 1,000 -Ij;11Q §; 1,120: 1,130. 1,150 


GENERAL CONVERSION FACTORS 


The following table of conversion factors departs from the 
conventional form in that the factors are given to as mary 
significant figures as possible. While in most cases four or five 
significant figures are sufficient, it has been the author’s experi- 
ence that greater accuracy is often required. It is not intended 
that the factors be used as given, unless such accuracy is required. 
For example, the exact conversion factor from cubic inches to 
cubic centimeters is 16.3871624, but 16.39 or even 16.4 is often 
close enough. 


Fundamental conversion factors. 


I meter = 39.37 inches (Act of U. S. Congress, 
28 July, 1866) 
3,.2808333 feet 


I 


I pound = 453.5924277 grams (International Bu- 
reau of Weights and Measures, July 


1893) 


Specific weight of dry air with normal CO, content 
at 760 mm. Hg. and 0° C. = .0012930 gr./cm. 


Specific weight of mercury at 
is OMwe SE Cta es MR eS WON Be i 


Standard gravity (Inter- 
national) g. 


. = 980.665 
= 32.174 ft./sec.? 
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MULTIPLY 


Atmospheres 


Bars 
Be toa. (mean) 


sé 
4c 


Centimeters (cm.) 
“ce 


cm. of mercury 
““ 


cm. per sec. 
cubic centimeters 
“cc 


cubic feet 
se 


cubic feet per min, 

be 4 
cubic feet of water 
cubic inches 

“ee 

“ce 


cubic meters 
oe 


cubic yards 
“ec 


Degrees (arc) 
ynes 
dynes per sq. cm. 


Ergs 


a aR a SE Aa se 


APPENDIX 


1,013,250. 


1.0 


777-98 
1,054.8 


107.560 


- 393700 

. 0328033 
5- 352391 

- 4460326 

. 193368 

27. 84507 
135-9510 

. 0328083 

- 000999973 

.06102338 

1,728.0 

I 


9 
7.480519 
28,317.017 

28.31625 
.028317017 
-471704 
. 028317 

62. 42833 

16. 3871624. 
.0163876 

1/231 
106 
61,023. 3753 
35-3144548 
I. 307943 
264.170 
2 


Tis 
- 78455945 
017453292 
. OO101972 
1.0 


1 U9) 
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CONVERSION FACTORS 


To OBTAIN 


cm. mercury 
inches mercury 

feet of water 

kilograms per sq. meter 
pounds per sq. in. 
pounds per sq. ft. 

bars 


dynes per sq. cm. 
foot-pounds 
joules 
kilogram-calories 
kilogram-meters 


inches 

feet 

inches of water 

feet of water 
pounds per sq. in. 
pounds per sq. ft. 
kilograms per sq. meter 
feet per sec. 

liters 

cubic inches 

cubic inches 

cubic yards 

gallons 

cubic centimeters 
liters 

cubic meters 

liters per sec. 

cubic meters per min. 
pounds 
cubic centimeters 
liters 

gallons 

cubic centimeters 
cubic inches 

cubic feet 

cubic yards 

gallons 

cubic feet 

cubic meters 


radians 
grams 
bars 


dyne-centimeters 
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CONVERSION FACTORS 


a ————— 
——— EEE 


MULTIPLY By To OBTAIN 
Fathoms 6.0 feet 
es 1.82880 meters 
feet 12.0 inches 
ue 1/3 yards 
Ke 30. 4800613 centimeters 
ie . 3048006 meters 
feet of water .029500 atmospheres 
ae - 433530 pounds per sq. in. 
us 62. 428327 pounds per sq. ft. 
se 304. 8006 kilograms per sq. meter 
36 . 882671 inches of mercury 
ut . 24199 centimeters of mercury 
feet per min. . 0113636 miles per hr. 
ss . 018288 kilometers per hr. 
Oe . 508001 centimeters per sec. 
feet per sec. .681818 miles per hr. 
a 1.09728220 kilometers per hr. 
ce 30. 48006 centimeters per sec. 
“s . 3048006 meters per sec. 
se - 5920858 knots 
foot-pounds . 138255 meter-kilograms 
foot-pounds per min. 1/33,000 horsepower 
foot-pounds per sec. 1/550 horsepower 
Gallons 231.0 cubic inches 
ue . 133680 cubic feet R 
4 3.785332 liters 
be . 832680 imperial gallons 
gallons-Imperial 1.20094 gallons 
ee ue 227.4176 cubic inches 
ue ve 4.5459631 liters 
grams . 0647988 grams 
cs 15.43236 grains 
ae - 0352739 ounces 
es . 0022046223 | pounds 
ue 1,000.0 milligrams 
ae .OOI kilograms 
ae 980.665 dynes 
gram-calories . 0039685 Bet. 
gram-centimeters 980. 665 ergs 
grams per cm. O.1 kilograms per meter 
as .06719702 pounds per foot 
ue .0055914 pounds per inch 
grams per cu. cm. 1,000.0 kilograms per cni. 
“ 62.42833 pounds per cu. ft. 
Horsepower 33,000.0 foot-pounds per min. 


550.0 foot-pounds per sec. 
76.04039 kilogram-meters per sec. 
“ . 986318 metric horsepower 


66 
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CONVERSION FACTORS 


MULTIPLY 


horsepower, metric 
“ce 


horsepower-hours 
“é 
cc 

Inches 

inches of mercury 
“cc 


6a 
sé 
“6 
6s 


inches of water 
“cé 


Joules 


Kilograms 
“é 
“6 
kilogram-calories 
oe 
“6 


kilogram-meters 
4é 
kilograms per cu. meter 
oe sé 


kilograms per meter 
kilograms per sq. meter 
oc ‘ 


kilometers 
4c 
sé 

kilometers per hr. 
oe 


‘ 66 


66 66 


knots 


By 


75.0 
1.013872 


2,545.06 


1,980,000. 


273,745-4 


2.54000508 
.0334211 
13.5951 
. 132925 
-4911570 
70. 72661 
345.3162 
-0735559 
. 1868324 
.0361275 
5.202360 
25. 400051 


107 


- 7375006 
. 1019716 


- 


2.20462234 
35- 273957 


7. 2329983 
g.80665 X 107 
. 06242833 

Oo! 


.6719702 

. 00142234 
. 2048169 

. 00289590 
. 003280833 
I 


O. 
3,280. 833 
. 6213700 
- 539553 
. 9113426 
.6213700 
+2777 
- 539553 
Ce) 


. 688944 
- 151553 
- 853249 
- 514791 


SS = eS 


kilogram-meters per sec. 


To OBTAIN 


horsepower 

etre Abie 
foot-pounds 
kilogram-meters 


centimeters 
atmospheres 

inches of water 

feet of water 

pounds per sq. in. 
pounds per sq. ft. 
kilograms per sq. meter 
inches of mercury 
centimeters of mercury 
pounds per sq. in. 
pounds per sq. ft. 
kilograms per sq. meter 


ergs 
foot-pounds 
kilogram-meters 


foot- pounds 
kilogram-meters 
foot-pounds 

ergs 

pounds per cu. ft. 
grams per cu. cm. 
pounds per ft. 
pounds per sq. in. 
pounds per sq. ft. 
inches of mercury 
feet of water 
grams per sq. cm. 
feet 

miles 

nautical miles 
feet per sec. 

miles per hr. 
meters per sec. 
knots 

nautical miles per hr. 
feet per sec. 

miles per hr. 
kilometers per hr. 
meters per sec. 


SS 
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CONVERSION FACTORS 


MULTIPLY By To OBTAIN 
Liters 1,000. 027 cubic centimeters 
4 61.02503 cubic inches 
es .035315411I | cubic feet 
+ . 264178 gallons 
“ .219975 imperial gallons 
Meters 39.37 inches 
“ 3.280833 feet 
“ 1.093611 yards 
meters per sec. 3.280833 feet per sec. 
. 2.2369317 miles per hr. 
+ 3.600 kilometers per hr. 
microns .OO00T centimeters 
miles 5,280.0 feet 
> 1.609347 kilometers 
e . 8683925 nautical miles 
miles per hr. 1.46666 feet per sec. 
we - 4470409 meters per sec. 
a 1.609347 kilometers per hr, 
‘ . 8683925 knots 
miles per hr. squared Pei seee! feet per sec. squared 
mills .OOI inches 
. 025400 millimeters 
Nautical miles 6,080. 20 feet 
I. 1515530 miles 
. 1,853.2486 meters 
Ounces 1/16 pounds 
28. 349527 grams 
ounces per sq. yd. 33.906096 grams per sq. meter 
Poundals .0310810 pounds 
¥ 13,825. 561 dynes 
pounds 453. 5924277 grams 
. -45359243 | kilograms 
16 ounces 
ss 32.174 poundals 
pounds-feet . 1382552 kilogram-meters 
pounds per ft. 1. 4881612 kilograms per meter 
pounds per cu. ft. 16. 018369 kilograms per cu. meter 
< ae .016018369 | grams per cu. cm. 
pounds per cu. in. 1,728 pounds per cu. ft. 
be a 27.6797424 grams per cu. cm. 
pounds per sq. ft. . 192220 inches of water 
if ns 4.8824088 kilograms per sq. meter 
pounds per sq. in. 2.036009 inches of mercury 
. rk 2.306645 feet of water 
: : ¥ .0680457 atmospheres 
703. 06687 kilograms per sq. meter 


ee ee eee 


MULTIPLY 


Quarts (dry) 
pea (Hatid) 


Radians 
radians per sec. 
6e 
oe 


revolutions 


revolutions per min. 


Slugs : 
square centimeters 


square feet 
“cc 
square inches 
“é 


square kilometers 
square meters 
ac 


square miles 
square yards 


Tons, long 
4é 
tons, short 
4c“ 


tons, metric 
6c 


Yards 


APPENDIX 


CONVERSION FACTORS 


By 


. 092903412 
645.162581 

6.45162581 

- 3861006 

10. 76386736 

1.1959853 
2.59000 

. 8361307 


1,016.047 
907. 1849 


2,204.622 


91.44018 
-9144018 


To OBTAIN 


cubic inches 
cubic inches 


degrees (arc) 
degrees per sec. 
revolutions per sec. 
revolutions per min. 
radians 

radians per sec. 


pounds 

square inch 
square feet 

square centimeters 
square meters 
square millimeters 
square centimeters 
square miles 
square feet 

square yards 
square kilometers 
square meters 


pounds 
kilograms 
pounds 
kilograms 
kilograms 
pounds 


centimeters 
meters 
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USEFUL FORMULAS 


Powers and Roots. 


(Peo Se am. gq" = q(m+n) 
a” 
a” = g(m—n) (a™)n me (an)m = qmn 
a 
a= Ve at = Vow 
i 
(a™ n — Gee 
a\" a” = 
ae = a" hn (5) Sgt a*b-* 
fhe. b oe 
(ab)" = a"b" = Vab 
Logarithms. 
Let 67=N (6 > 1.0) 
Then: logs N = x log,b = 1.0 
log,sI =o log,0 = —© 


log, MN = log,M + log,N 


logs 7 = log,M — logsN 


log, N’ = r- log,N 


log, V/ N’ = <- logs 


€ = 2.71828183 


loge2 = 0.69314718 


loge 10 = 2.30258509 
logio€ = 0.43429448 


log. N = log, 10- log,,. N 
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Derivatives. 


FUNCTION 


Sis 


u” . ° ° 


sin u 


cos u wilh ee 


tan u ees 


cot u 


c is any constant. 
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DERIVATIVE 


ep YE Sy ae du 
er eee es. ee OTe "Gg 1 Milos Te 


du 
PBT t Sut CUMS (COS Dy sitert 


dx 


CS CH 2 


uw and v are any functions of x, 
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Integrals. 


FUNCTION INTEGRAL 

The aig) eater Te. ae +¢ @un 
ou Soa logeu +e 

a“- du E aor +c 

et - dy : a + ¢ 

A) ern ae Se Se OE Ls . uw— fo-du 


u and v are any functions of x. 


INDEX 


A 


Absolute angles of attack, 6, 94-95, 
100-I0I 
Absolute ceiling, 155-157, 166, 168- 
170, 235-238 
definition of, 155-157 
estimation of, 235-238 
graphical solution for, 155-170 
Absolute coefficients, 4, 147 
Adjustable pitch propellers, 132-133 
Aileron area, 107-110 
balanced, 114-117 
efficiency, 107-109 
moments, 19-21, [10-114 
types, IIO-1II 
Air-cooled engines, drag of, 88 
Airfoil sections for tail surfaces, 102- 
103 
Airplane axes, 92 
model tests, 53-64 
parts, drag of, 87 
Airship, C-class offsets, 73 
Air-speed indicators, calibration of, 
196 
Air-speed, 
for best climb, 162-164, 183-186, 
197-200 
determination of, 199-200 
most economical, 201-202 
Altitude, 
climbed in given time, 159, I65- 
166, 168-170, 239-244 
effect on specific fuel consumption, 
124, 127 
variation of B.Hp. with, 120-123 
variation of T.Hp. with, 137-139 
Aneroid, 199-200 
Angle of attack, 
absolute, 6, 94-95, IO0O-I01 
blade setting, 132-133 
glide, 213-220 
induced, 30-32 
tail setting, 100 
Areas of control surfaces, 93-119 
ailerons, 107-110 
elevators, 102 
rudders, 106 
Aspect ratio, 29, 30, 42-45, 65, 67, 
95-99, IOI, 163, 172-186, 188- 
189, 231-238 


Aspect ratio—Continued 
effect on, 
absolute ceiling, 179-180, 237- 
238 
air speed for climb, 183-186 
downwash factor, 96 
horizontal tail area, 97-98 
maximum L/D, 180-182 
maximum speed, 173-175, 178 
power required, 173-175 
rate of climb curve, 95 
slope of lift curve, 95 
speed for minimum power, I8I, 
183, 
speed range constant, 182-184 
Axes, airplane, 92 


B 


Balanced controls, 114-119 
Bank, proper angle of, 217-220 
Bendemann dynamometer, 
with, 122-123 
Bernoulli’s theorem, 24-25 
Biplane, drag of, 32-43 
proportions of most efficient, 35, 43 
Blade setting, propeller, 132-133 
Body interference, effect on propeller, 
130-131 
Bombing range, 208-209 
Brake horsepower, 
general curves for, 120-122 
variation with altitude, 120-123 
variation with air speed, 139-141 
Breguet’s formulas for range and 
endurance, 204-208 


C 


Cable, drag of, 75-78 
Calibration of air-speed indicators, 
196 
Camber, mean, 14-15 
Ceiling, 
absolute, 137, 155-157, 166, 168- 
170, 235-238 
service, 137, 157-158, 166, 168-169, 
238-239 
Cellular radiators, drag of, 86 
Center of gravity location, 94, 98-102 
Center of pressure, 6-7, 14, 45-46, 
94-95, IOI 
Centrifugal force, 216, 220 
Characteristics, wing section, 5 


tests 


283 


284 


Circling flight, 216-220 
power required in, 217-219 
Circular discs, drag of, 66-68 
Circulation, 26-29 
Climb, 
air speed in, 162, 164, 183-186, 
197-200 
effect of assumed rates of, 168-170 
effect of decrease in weight, 170 
effect of temperature on, 189-192, 
198-199 
graphical solution for rate of, 145, 
. 154-155 
initial rate of, 155, 161, 165-166, 
168-171, 234-235 
in ten minutes, 242, 254 
reduction to standard, 
194-195 
sawtooth, 199-200 
time of, 158-159, 165-166, 168-170, 
239-244 
true rate of, 165-171 
Coefficients, 3, 4 
Compressibility, 6, 50-51, 53-54 
Compression ratio, effect on fuel 
economy, 123-125 
Control, 93-119 
Control horns, drag of, 81 
Controls, balanced, 114-119 
Conversion factors, 
for absolute coefficients, 4 
general, 273-279 
Corresponding speeds, 254-255 
Cut-out trailing edge, 22-23 
Cylinders, drag of, 68-70 


D 


189-192, 


Density, 
altitude, 187-188 
effect on B.Hp., 121-123 
effect on speed in a dive, 222 
ratio, 187-188 
standard, 4- 
Derivatives, table of, 281 
Diameter of propeller, 127-132 
Dihedral, 93 
longitudinal, 99 
Directional pee 93, 103-106 
Discs, drag of, 66-6 
Dive, speed in a, oe 221 
Diving start, effect on measured 
speed, 220-221 
Downwash, 28-29, 96, 98 
Drag coefficients, paca 3, 4, 65 
for wires, 75-76, 78 
Drag, 
biplane, 32-43. 
correction for, in model tests, 55-57 


INDEX 


Drag—Continued 
general equation for, 50-51, 53-54 
induced, 28-45, 148-149 
profile, 7, 11, 13, 29, 45, 71-72, 148 
Dynamic pressure, 4, 5 
stability, 92-93 


E 


Economical air speed, 201-202, 209- 
212 
effect of wind, on 209-212 
Efficiency, propeller, 129-133 
effect of body interference on, 
130-131 
effect of tip speed on, 130, 132 
general curve for, 128-130 
maximum, 130-133 
Elevators, 
area of, 102 
balanced, 114-117 
Ellipsoids, drag of, 68 
Endurance, 201-212, 245-250 
effect of mixture control on, 204- 
208 
estimation of, 245-250 
formulas for, 203-208 
maximum, 204, 249-250 
Engines, 
corrections for test stand data, 
124-127 
drag of air-cooled, 87 
general power curves, 
140-141 
specific fuel consumption, 123-127, 
201-211, 245 
variation of B.Hp. with altitude, 
120-123 
Equivalent altitude in performance 
reduction, 189 
Estimation of performance, 232-252 
gc power for climb, 146, 154, 155) 
161 


120-122, 


F 


Factors, conversion, 4, 273-279 
Fin area, vertical, 103-106 
Fineness ratio, 72, 74, 85-86 
Fittings, drag of, 75, 79-81 
Flaps, 15-21, 23 
Flat plates, apes of, 65-66 
Flat spin, 230-231 
Flight test data, 122-123, 138 
Flight testing, 196-200 
Floats, 

drag of, 83-85, 148 

excess buoyancy, 83 

seaplane, 253-265 


INDEX 


Floats—Continued 
take-off calculations, 258-262 
wing tip, 85, 148, 264 
Flutter from overbalance, 110, 116 
Flying boat hulls, drag of, 85, 148 
Fore-and-aft slot, 47-49 
Fore-and-aft spacing, effect on wire 
_ drag, 77 
Frictional drag, 70-71 
Frise balance, 116-117 
Fuel consumption, 123-124, 126, 201- 
211, 245 
effect of compression ratio on, 
123-125 
part throttle, 123-124 
variation with altitude, 124, 126 
Fuselage drag, 81-83, 148 


G 


Gap/chord ratio, 163 
General curves for 
B.Hp., 120-122, 140-141 
power coefficients, 135-137 
propeller efficiency, 128-130 
r.p.m., I4I-143 
T.Hp., 142-144 
General equation for drag, 50-51, 
53754 
Getaway speed, 255, 257-262 
Gliding flight, 213-220 
rate of descent in, 214 
with power, 214~216 
without power, 213-214 
Graphical calculation of performance, 
146, 154-159 
Ground run, 223-227 
effect of wind on, 224-227 
Ground speed, 211 
Guns, drag of, 98 


H 


Handley Page, 
balance, 116-117 
slots, 19, 21, 23 

Hemispheres, drag of, 69 

Hinge moments for flaps, 17, 19, 20, 
IIO-I12 


Horizontal flight, power required in, 


147-149 |. 
Horizontal tail area, 94-103 
equations for, 95, 100 
Horns, drag of control, 81 
Hulls, drag of flying boat, 85 
Hydrodynamical definitions, 24-28 


I 


Inclination, _ 
effect on wire drag, 77 
of flight path, 213-216 


285 


Indicated air speed, 196, 199~200 
Induced angle of attack, 30-32 
Induced drag, 28-45, 148-139 

in performance calculations, 148- 


149 
Induced power, 175 
Induction, 29 
Initial rate of climb, 155, 161, 165- 
166, 168-171, 176-177, 234-235 
effect of aspect ratio on, 176-177 
estimation of, 234-235 
Integrals, 282 
Interference, 8, 31-35, 49-50, 89-91, 
130-131 
wall, 49-50 
Trrotational motion in a fluid, 26-28 


L 


Lag in flight test instruments, 200 
Lamblin radiators, 86 
Landing gear, drag of, 88-89 
Landing run, 227-229 
Lateral stability, 92-93, 264 
Lean mixture, 123-125 
Lift coefficients, definition of, 3, 4 
Lift, increase due to flaps, 16, 17, 19, 
21, 23 
Limiting speed, 
airplane, 213, 221-223 
parachute, 69, 229, 230 
Location of c.G., 98-102 
Logarithms, relations for, 280 
Longitudinal dihedral, 79 
stability, 92-93, 94-102 


M 


Machine gun, drag of, 89 
Mass distribution, effect on spins, 231 


- Maximum endurance, 204, 249-250 


Maximum lift due to flaps, 19 
Maximum L/D, 180-181 
Maximum propeller efficiency, 130- 
132, 145 
Maximum range, 201-212, 245-249 
effect of mixture control on, 204- 
208 m 
effect of wind on, 209-212 
estimation of, 245-249 
formulas for, 203-208 
with a bomb load, 208-209 
Maximum speed, 
at altitude, 162, 164 
effect of aspect ratio on, 176, 178 
estimation of, 233-234 
flight tests, 197 
reduction to standard, 192-193 
Mean camber, 14-15 


286 


Metacenter, 254, 262 
Metacentric height, 254, 262-263 
Minimum power, 155, 181, 183 
Mixture control, 123-125 
Model basin, 
data from, 257-262 
description of, 256-257 
tests, 254-258 
Models, tests on airplane, 53-64 
performance from, 60-64 
Moment coefficient at zero lift, 98- 
101 
Moment curves, slope of, 58-60 
Moments, 
pitching, 57-59 
rolling, 59 
wing, 45-47 
yawing, 59-60 
Munk’s 
interference factor, 31-32 
span factor, 30, 31, 33-45 


N 


Nacelles, drag of, 81-83, 90-91, 148 
Negative thrust, 134-135 
Non-dimensional 
engine data, 120-122, 140-144 
propeller data, 128-130, 135-137 
Normal force coefficient, 46 


O 


Oil consumption, 
203-204, 208 


1B 


Paddle balance for control surfaces, 
116-117 
Parachutes, 69, 229-230 
Parallel middle body, 85 
Parasite drag, 
in performance calculations, 148- 
152, 172-185 
variation with. angle of attack, 
148-152 
Parasite drag coefficient, 172-185 
effect on absolute ceiling, 179-180 
air speed for climb, 183-186 
maximum L/D, 180-182 
maximum speed, 173-175, 178 
power required, 173-176, 178 
rate of climb, 176-177, 178 
speed for minimum power, 181, 
183 
speed range constant, 182-184 
Parasite drag data, 65-91, 172 
Performance, 
at altitude, 161 
from wind tunnel tests, 60-64 
reduction to standard, 187-194 


effect on range, 


INDEX 


Performance calculation, 146-164 
Performance estimation, 232-252 
absolute ceiling, 235-238 
initial rate of climb, 234-235 
maximum endurance, 250 
maximum speed, 232 
range and endurance, 245-250 
service ceiling, 238-239 
stalling speed, 232 
time of climb, 239-245 
Pitch, propeller,.129, 132-135 
Pitching moments, 57-59 
slope of curve of, 58-59, 94-95, 
100-102 
Polar curves, full scale, 172 
Porpoising of seaplane floats, 254-255 
Potential, velocity, 26 
Power absorbed by test club, 124-127 
Power available, 153-162 
at altitude, 158-159, 162 
example of calculation for, 160 
Power coefficient, curves of, 135-137 
Power curves, general, 120-122, 140- 


144 ; 
Power, effect of pressure and tem- 
perature on, 120-123, 187-189 
Power required, 63-64, 147-162 
at altitude, 159-160, 162 
example of calculation for, 160 
in horizontal flight, 63-64, 147-149, 
187 
in a turn, 217-219 
Powers and roots, 280 
Prandtl’s interference factors, 32-35 
Pressure altitude, 187-189, 192-193, 
199-200, 271 
Pressure, dynamic, 4 
effect on B.Hp., 120-123, ae 
Prismatic coefficient for seaplane 
floats, 265 
Profile res 7, ET, 13,29; 45) i= 725 
148 
Propeller, 
blade setting, 132-133 
diameter, 127-132 
drag coefficients for, 134-135 
efficiency, 128-132 
general curves for, 128-130, 138 
idling, 135 
maximum efficiency, 130-133, 145 
negative thrust, 134-135 
pitch, 129, 132-135 
static thrust, 133-134, 225-226 
R 
Radiators, drag of, 86 
Range, 
calculation of, 201-212 
estimation of, 245-249 


INDEX 


Rate of climb, 
by graphical methods, 145, 154- 
155, 161-163 
constants in equation for, 165-171 
effect of assumed linear, 168-170 
estimation of, 234-235 
exact equation for, 165 
general curve for, 239-241 
initial, 155, 161, 165-166, 168-171, 
234-235 a 
method of obtaining constants for, 
167-168 
slope of curve, 165, 167, 169-171, 
239-241 
true, 165-171 
Reduction of performance to stand- 
ard, 187-194 
Reynolds Number, 3, 5, 10, 11, 51-54, 
65-69, 71, 232 
determination of, 51-52 
effect on max. lift, 10, 232 
effect on min. drag, Io 
Rolling moments, 59, 112-114 
Rotation in a fluid, 27-28 
Round wire, drag of, 75-78 
r.p.m., 
of engine, 120, 122, 133, 139-141 
propeller, 127-128, 133, 139-143 
Rudder and fin area, 106 
Rudders, balanced, 114-119 
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Scale effect, 9-11, 54, 56-57, 71-74 
Sea level, 
power available at, 153-154 
power required at, 152 
Seaplane floats, 253-265 
corresponding speeds, 254-256 
drag of, 83-84, 148 
metacentric height, 262-264 
model basin data, 257-258, 260 
model basin tests, 256-257 
prismatic coefficients, 265 
section coefficients, 264-265 
stability longitudinal, 263-264 
stability transverse, 263-264 . 
take-off calculations, 258-259, 261- 
262 
Searchlight, drag of, 89 
Section coefficients for floats, 264-265 
Sections for 
tail surfaces, 102-103 
wings, 4-II 
Service ceiling, 157-158, 166, 168- 
169, 238-239 
definition of, 157 
estimation of, 238-239 
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Skin friction, 50, 70-72 
Slipstream corrections, 162-163 
Slope of lift curves, 94-96, 99-102 
table of, 101 
Slope of moment curves, 58-60, 94- 
_ 95, 100-106 - 
in pitch, 58-59, 94-95, 100-102 
in roll, 59 
in yaw, 59-60, 103-106 
Slot, 
fore-and-aft, 47-49 
Handley Page, 19, 21, 23 
Sound, velocity of, 50, 51, 53, 54, 273 
equation for, 273 
Spacing, effect on wire drag, 77 
Span Hoa 30, 31, 33-45, 148-149, 
I 


3 

at altitude, 162-164 

Specific fuel consumption, 123-126, 
201-211, 245 

Speed, maximum, 

effect of aspect ratio on, 176, 178 

estimation of, 233-234 

flight tests, 197 

reduction to standard, 192-193 
Speed, minimum, 61-62, 147, 232 
Speed for minimum power, 181-183 
Speed, stalling, 61-62, 147, 232 
Speed range, constant for, 182-184 

formula for, 182, 233-234 
Spheres, drag of, 68 
Spheroid, drag of, 68 
Spin, flat, 230-231 
Squashing in a turn, 219 
Stability airplane, 92-106 

dynamic, 92-93 

lateral, 92-93 

longitudinal, 92-102 

seaplane floats, 262-264 

static, 92-102 
Stabilizer area, 102 

setting, 99-102 
Stagger, 231 
Stalling speed, 61-62, 147, 232 
Standard atmosphere, 121, 138-139, 

187, 230, 267-272 

charts, 269-272 

definition of, 267 

formulas, 267-268 

reduction of performance to, 187- 


194 

Static stability, 92-102 

thrust, 133-134, 225-226 
Steady wind, effect on range, 211-212 
Stream function, 25 
Streamline, 25 

forms, 72-74, 81-82, 83-85 

wire, 77, 79-81 
Struts, 54, 72-74 
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Tail area, 
horizontal, 94-103 
vertical, 103-106 
Tail setting, equation for average, 100 
Tail skid, drag in landing, 228-229 
Tail surface, sections, 102-103 
plan form, 103-104 
Take-off run, 223-227 
approximation formulas for, 225- 
226 
effect of wind on, 224-227 
general curve for, 226-227 
seaplane, 258-262 
Temperature, effect on 
B.Hp., 120-123, 187-189 
maximum speed, 192-193 
rate of climb, 189-192 
Terminal velocity, 
of an airplane, 213, 221-223 
of a parachute, 229-230 
Tests, airplane model, 53-64 
Test-stand data, 124-127 
Three-dimensional flow, 26 
Throttling, effect on fuel consump- 
tion, 123-124 
Thrust, negative, 134-135 
Thrust power, 
calculation of, 147-160 
variation with altitude, 137-139 
variation with r.p.m. and V, 142- 


145 

Thrust, static, 133-134, 224-226 
Tip speed, 130, 132 
Trail angle, 12-15 
Trailing edge, cutting off, 22-23, 72 
Transverse stability, 263-264 
Trim, definition of, 254 
Turbulence, 50-51, 53-54 
Turn, 

power required in a, 217-219 

squashing in a, 219 
Turnbuckles, drag of, 75 
Two-dimensional flow, 25 


Vv 


Variable density wind tunnel, 3, 
11-13, 67-68 

Vector diagram, 57-58 

Velocity of sound, 50-54, 273 

Velocity potential, 26 
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Vertical tail area, 103-106 
Viscosity, formula for, 51 
Vortex theory, 27 
Vorticity, 27-28 
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Wall interference correction, 49-50 
Warping of wind tunnel models, 56 
Wave-making resistance, 50-54, 254 
Wheels, drag of, 88-89 
Wind, 
effect on economical speed, 209-212 
effect on range, 209-212 
Wind tunnel, 
variable density, 3 
wall interference, 49-50 
Wind tunnel tests, 
on airfoils, 5-6 
on airplane models, 53-64 
on float drag, 83-84 
Wing drag, 
induced, 28-45, 148-149 
in performance calculations, 148- 


149 

profile, 7, II, 13, 29, 45, 71-72, 148 
Wing section, 

analysis of data, 8 

drag, 71-72 

methods of plotting data, 6 

moment, 45-47 

selection of, 6-13 

stability, 94-95, 101 
Wing section characteristics, 

definition, 5 

methods of plotting, 6 
Wing tips, 23 
Wires, 

cable, 75-78 

streamline, 77, 79-81 
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Yaw, effect on drag of streamline 
wire, 81 
Yawing moments, 59, 60, 103-106 


Z 


Zero lift, 46-47, 98-99, IoI 
moment coefficient at, 98-101 
Zero moment, 46-47 
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